Mathematica 9 is now available
Student Support Forum
-----
Student Support Forum: 'equation err' topicStudent Support Forum > General > "equation err"

Next Comment >Help | Reply To Topic
Author Comment/Response
zy
10/28/10 8:54pm

u[\[Rho]_] = \[Rho]*2;
K[T_] = T + 1;
n[T_] = T + 2;
ReNum[ \[Rho]_, K_, n_, u_] = ( \[Rho] + K + n + u);
\[Lambda]Value[ReNum_, n_] = (ReNum + n);
p\[Rho]p[p_, T_] = (p + T);
p\[Rho]T[p_, T_] = (p - T);
z1 = 0;
z2 = 2.0;
\[Rho]0 = 1.05;
p0 = 5.5;
Solequpre[\[Rho]_, T_, p_, z1_, z2_] :=
Module[{u0, K0, n0, \[Lambda]Value0, ReNum0, p\[Rho]p0, p\[Rho]T0},
(*---微分方程计算参数准备----*)
u0 = u[\[Rho]];
K0 = K[T];
n0 = n[T];
ReNum0 = ReNum[ \[Rho], K0, n0, u0];
\[Lambda]Value0 = \[Lambda]Value[ReNum0, n0];
p\[Rho]p0 = p\[Rho]p[p, T];
p\[Rho]T0 = p\[Rho]T[p, T];
(*----------*)

NDSolve[
{p'[z] ==
9.81*\[Rho][z] + (u0)^2*\[Rho]'[z] - (\[Lambda]Value0*\[Rho][
z]*(u0)^2)/(2*De),(*质量方程*)
\[Rho]'[z] == (p\[Rho]p0)*p'[z] + (p\[Rho]T0)*0.032,
p[z1] == p0, \[Rho][z1] == J0 },(*气体方程*)
{p[z], \[Rho][z]},(*压力、密度函数*)
{z, z1, z2}] ;

];

Solequpre[1.05, 283, 5, 0, 2.5]

the result:
NDSolve::deqn:
Equation or list of equations expected instead of False in the first argument {0==-33148.41.05[z],False,5[0]==5.5,1.05[0]==1.05}.

URL: ,

Subject (listing for 'equation err')
Author Date Posted
equation err zy 10/28/10 8:54pm
Re: equation err Peter Pein 10/30/10 5:36pm
Re: Re: equation err zy 11/01/10 08:59am
Next Comment >Help | Reply To Topic