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

1 view (last 30 days)
MINATI on 2 May 2019
Closed: Walter Roberson on 3 May 2019
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]
##### 2 CommentsShowHide 1 older comment
MINATI on 3 May 2019
ok
But the basic idea is to share the knowledge.
Pl. share @
minatipatra456@gmail.com