### Info

이 질문은 마감되었습니다. 편집하거나 답변을 올리려면 질문을 다시 여십시오.

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

조회 수: 1(최근 30일)
MINATI 2 May 2019
마감: Walter Roberson 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표시숨기기 이전 댓글 수: 1
MINATI 3 May 2019
ok
But the basic idea is to share the knowledge.
Pl. share @
minatipatra456@gmail.com

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by