Can anyone help me to run this coding?

RekhaAl2O3CuSiO2TiO2Water_My14PaperCoding()

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 NaN 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -4.95261114166381e-14

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -1.70783380211543e-14

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -1.03182089464865e-14

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -7.39216407886729e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -5.75901780523335e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -4.71691309070662e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -3.99416219840075e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -3.46347054135655e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -3.05726185938896e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -2.73633440296135e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -2.47638306282591e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -2.26153736683265e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -2.08099471768041e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -1.92714705636603e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * -1.49011611938477e-14 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 1.49011611938477e-14 Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * 0 0 Inf 1.49011611938477e-14 -1.79448129415503e-15

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 0 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

0 0 -Inf 0 0

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 1.49011611938477e-14 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

1.0e+00 * NaN 0 NaN -1.49011611938477e-14 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

ff = 5×1

NaN 0 NaN 0 NaN

<mw-icon class=""></mw-icon>

<mw-icon class=""></mw-icon>

Error using bvp4c (line 196)
Unable to solve the collocation equations -- a singular Jacobian encountered.

Error in solution>RekhaAl2O3CuSiO2TiO2Water_My14PaperCoding (line 118)
sol = bvp4c (@OdeBVP, @OdeBC, solinit);

function RekhaAl2O3CuSiO2TiO2Water_My14PaperCoding
clear
close all
format long g

%Define all Parameters
global N r s B G M K Pr Q p Nr V t a b c d A1 A2 A3 A4 A5 P1 P2 P3 P4 P5 C1 C2 C3 C4 C5 
global D1 D2 D3 D4 D5 B1 B2 B3 B4 B5 B6 B7 B8 K1 K2 K3 K4 K5 K6 K7 K8 H1 H2 H3 H4 H5 H6 R1 R2 R3 R4 R5 E1 E2 E3 E4 E5


%Boundary layer thickness & stepsize
etaMin=0;
etaMax1=15;
stepsize1=etaMax1;

%Input for the Parameters
N=3;                % shape factor
r=1;                      
s=1;           
B=1;                     
G=1;                     
M=2;                     %Schmidt parameter 
K=1;                    %Porosity parameter 
Pr=1;                    %Activation Energy parameter 
Q=1;                   %Prandtl parameter 
p=0.001;                   %Brownian parameter 
for Nr=[0.1 0.2 0.3]                    %Chemical reaction parameter 
    V=1;
    t=pi/6;
    
    
    a=0.05;                  %phi1-1st nanopaeticle concentration (MgO)
    b=0.15;                   %phi2-2nd nanoparticle concentration (TiO2)
    c=0.25;                   %phi2-2nd nanoparticle concentration (CoFe2O4)
    d=0.05;
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    %%%%%%%%%%%%%%%%%%%1st nanopaeticle properties (Al2O3)%%%%%%%%%%%%%%%%%%%
    P1=3970;                     %P=rho=density 
    C1=765;                      %C=specific heat 
    K1=40;                       %K=thermal conductivity 
    B1=36.9*10^6;             %electrical condutivity 
    R1=0.85*(10)^(-5);          % thermal exapnsion coeffcient
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    %%%%%%%%%%%%%%%%%%%2nd nanopaeticle properties (Cu)%%%%%%%%%%%%%%%%%%%
    P2=8933;                     %P=rho=density 
    C2=385;                    %C=specific heat  
    K2=400;                   %K=thermal conductivity  
    B2=5.96*(10)^7;              %electrical condutivity 
    R2=1.67*(10)^(-5);          % thermal exapnsion coeffcient
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%3rd nanoparticles (SiO2) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    P3=2200;                     %P=rho=density 
    C3=754;                      %C=specific heat 
    K3=1.4013;                      %K=thermal conductivity  
    B3=3.5*(10)^6;              %electrical condutivity 
    R3=5.5*(10)^(-7);          % thermal exapnsion coeffcient
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    %%%%%%%%%%%%%%%%%%% 4th nanoparticles (TiO2)%%%%%%%%%%%%%%%%%%%
    P4=4250;                    % rho density 
    C4=686.2;                     %C=specific heat 
    K4=8.9538;                    %K=thermal conductivity 
    B4=2.4*(10)^6;                     %electrical condutivity 
    R4=0.90*(10)^(-5);          % thermal exapnsion coeffcient
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    %%%%%%%%%%%%%%%%%%% Base fluid(Water)%%%%%%%%%%%%%%%%%%%
    P5=997.1;                    %P=rho=density 
    C5=4179;                     %C=specific heat 
    K8=0.613;                    %K=thermal conductivity 
    B8=0.05;                     %electrical condutivity 
    R5=21*(10)^(-5);             % thermal exapnsion coeffcient
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%multiplier%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    H1=((1-d)*((1-c)*(((1-a)*(1-b)*P5+a*P1)+b*P2)+c*P3)+d*P4)/P5;                    %(Rho)mnf/(Rho)f
    
    H2=((((1-a)^(2.5))*((1-b)^(2.5))*((1-c)^(2.5))*((1-d)^(2.5))))^(-1);                            %(mu)mnf/(mu)f
    
    D1=P1*C1;                                                                  %rho*cp nanopartocale1
    D2=P2*C2;                                                                  %rho*cp nanopartocale2
    D3=P3*C3;                                                                  %rho*cp nanopartocale3
    D4=P4*C4;
    D5=P5*C5;                                                                  %rho*cp base fluid
    H3=((1-d)*((1-c)*(((1-a)*(1-b)*D5+a*D1)+b*D2)+c*D3)+d*D4)/D5;                   %(Rho*Cp)khnf/(Rho*Cp)bf
    
    E1=P1*R1;                                                                  %rho*cp nanopartocale1
    E2=P2*R2;                                                                  %rho*cp nanopartocale2
    E3=P3*R3;                                                                  %rho*cp nanopartocale3
    E4=P4*R4;
    E5=P5*R5;                                                                  %rho*cp base fluid
    H5=((1-d)*((1-c)*(((1-a)*(1-b)*E5+a*E1)+b*E2)+c*E3)+d*E4)/E5;
    
    B7=(B1+(N-1)*B8-(N-1)*a*(B8-B1))*B8/(B1+(N-1)*B8+a*(B8-B1));          %(k)mnf/(k)f
    B5=(B2+(N-1)*B7-(N-1)*b*(B7-B2))*B7/(B2+(N-1)*B7+b*(B7-B2));              %(k)hnf
    B6=(B3+(N-1)*B5-(N-1)*c*(B5-B3))*B5/(B3+(N-1)*B5+c*(B5-B3)); 
    H4=(B4+(N-1)*B6-(N-1)*d*(B6-B4))*B6/(B4+(N-1)*B6+d*(B6-B4))*B8;
    
    K7=(K1+(N-1)*K8-(N-1)*a*(K8-K1))*K8/(K1+(N-1)*K8+a*(K8-K1));          %(k)mnf/(k)f
    K5=(K2+(N-1)*K7-(N-1)*b*(K7-K2))*K7/(K2+(N-1)*K7+b*(K7-K2));              %(k)hnf
    K6=(K3+(N-1)*K5-(N-1)*c*(K5-K3))*K5/(K3+(N-1)*K5+c*(K5-K3)); 
    H6=(K4+(N-1)*K6-(N-1)*d*(K6-K4))*K6/(K4+(N-1)*K6+d*(K6-K4))*K8;
    
    A1=H2/H1;
    A2=H5/H1;
    A3=H4/H1;
    A4=H6/H3;
    A5=(H3)^(-1);
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%first solution%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    solinit = bvpinit (linspace(etaMin, etaMax1, stepsize1), @OdeInit1);
    sol = bvp4c (@OdeBVP, @OdeBC, solinit);
    eta = linspace (etaMin, etaMax1, stepsize1);
    y = deval (sol, eta);
    
    figure(1)   %velocity przfile
    plot(sol.x,sol.y(2,:),'linewidth',1.5)
    xlabel('\eta')
    ylabel('f(\eta)')
    %xlim([0 10])
    % ylim([0 1])
    hold on
    
    figure(2)   %Velocity profile
    plot(sol.x,sol.y(5,:),'linewidth',1.5)
    xlabel('\eta')
    ylabel('g(\eta)')
    %xlim([0 10])
    %ylim([0 1])
    hold on
    
    
    %Displaying the output for first solution
    fprintf('\nFirst solution:\n');
    fprintf('f''(0) = %7.9f\n',y(3));     %reduced skin friction
    fprintf('-theta(0)= %7.9f\n',-y(6));  %reduced local Nusselt number
    fprintf('CfxRe(1/2) = %7.9f\n',H2*y(3));   %skin friction
    fprintf('NuxRe(-1/2) = %7.9f\n',(-1)*H6*y(6));     %local Nusselt number
    fprintf('\n');
    % fprintf('%3.2f %7.6f %7.6f\n', B,A1*((2*O*(r^2)/v)^(-1/2))*(((y(3))^2+(y(5))^2)^(1/2)), (-1)*r*((2*O/v)^(1/2))*(A4+(4/3)*Nr)*y(7), (-1)*y(9))
end


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Define the ODE function
    function ff = OdeBVP (x, y)
        
        ff = [y(3)
            y(4)
            ((A1*(1+2*r*x)*y(4)+2*A1*r*y(3)+(1+s)*(y(1)*y(3)-((y(2))^2))+A1*r*B*(y(2)*y(3)-3*y(1)*y(4))+A1*B*(1+2*r*x)*((y(2))^2)+A2*G*cos(t)*y(5)-(A3*M+A1*K)*y(2)))/(A1*B*(1+2*r*x)*y(1))
            y(5)
            ((Pr*(y(2)*y(5)-y(1)*y(5)+p*(((y(2))^2)*y(5)-y(1)*y(2)*y(5)-y(1)*y(3)*y(5)))-A2*Q*y(5)-2*A1*r*y(5))/(A4*(1+2*r*x)-Pr*p*((y(1))^2)+A5*Nr))
            ]
    end

%Define the boundary condition 
    function res = OdeBC (ya, yb)
        
        res = [ya(1)-0
            ya(2)-1 
            ya(5)-(1-V)
            yb(2)-0
            yb(5)-0];
    end


%Setting the initial guess for first solution
    function v = OdeInit1 (~)
        
        v = [0
            0
            0
            0
            0];
    end

end

Answer 1

[James Blanchard]

I would suggest that you start by not using an initial guess of zero for all of your variables. I think this is leading to some divide by zero issues. A better initial guess will give you something to work with. You should have a general idea of what values you expect for all of your variables, so use that knowledge to select an initial guess.

Answer 2

[Torsten]

OdeBVP returns NaN for ff(3) at the start (see above). You will have to find out the reason (initial values for y, x, parameters used ... ?).