Graph is not coming with this MATHEMATICA code (new to this)

TypeEQ = 2;

NumEQ = 3;

f[1, z_, {f_, g_, s_}, sigma_] := D[f, {z, 2}] + Gr*Pr*g - Nr*s + sigma;

f[2, z_, {f_, g_, s_}, sigma_] := D[g, {z, 2}] + Nb*D[g, z]*D[s, z] + Nt*(D[g, z])^2 - f;

f[3, z_, {f_, g_, s_}, sigma_] := D[s, {z, 2}] + Nt/Nb*D[f, {z, 2}] - Le*f;

NumBC = 7;

BC[1, z_, {f_, g_, s_}] := f /. z -> -1;

BC[2, z_, {f_, g_, s_}] := f /. z -> 1;

BC[3, z_, {f_, g_, s_}] := g /. z -> -1;

BC[4, z_, {f_, g_, s_}] := g /. z -> 1;

BC[5, z_, {f_, g_, s_}] := s /. z -> -1;

BC[6, z_, {f_, g_, s_}] := s /. z -> 1;

BC[7, z_, {f_, g_, s_}] := Integrate[f, {z, 0, 1}] - Ra*Pr;

zL[1] = -1; zR[1] = 1;

zL[2] = -1; zR[2] = 1;

zL[3] = -1; zR[3] = 1;

U[1, 0] = eps1 - 3/2*(-25 + 4eps1)z^2 + 5/2*(-15 + 2eps1)*z^4;

U[2, 0] = eps2*(1 - z^2);

U[3, 0] = eps3*(1 - z^2);

L[1, u_] := D[u, {z, 2}]; L[2, u_] := D[u, {z, 2}]; L[3, u_] := D[u, {z, 2}];

Nr = 3/20; Nt = 1/20; Nb = 1/20; Le = 10; Gr = 5; Pr = 1; Ra = 5;

c0[1] = c0[2] = c0[3] = h;

GetOptiVar[3, {}, {eps1, eps2, eps3, h}];

BVPh[1, 10]

[KSSV]

This is a MATLAB forum...how yu expect us to help for mathematica?

[MINATI]

ok

But the basic idea is to share the knowledge.

Pl. share @

[e-mail]