Beginner--Help on using FindRoot to solve the system of equations
- To: mathgroup at smc.vnet.net
- Subject: [mg66190] Beginner--Help on using FindRoot to solve the system of equations
- From: xiaochu at che.utexas.edu
- Date: Wed, 3 May 2006 02:45:50 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
This is a code to solve vapor-liquid equilibrium by van der Waals Equation of State.
I don't why this code is not working, please help, thanks very much!
(* Name of EOS *)
EOSName = "Van der Waals";
Ttilde = .;
Dtilde = .;
Z = 3/(3-Dtilde)-9*Dtilde/8/Ttilde;
(* Related variables *)
Psi1 = (-9*Dtilde)/(8*Ttilde) - Log[3-Dtilde];
lnB = (1-Z)-Psi1;
mu = -lnB+Log[Dtilde];
P = Dr Ttilde Z;
DtildeG=10^-14;
DtildeL=2.91;
Dtilde1=.;
Dtilde2=.;
mu1=mu/.Dtilde->Dtilde1;
mu2=mu/.Dtilde->Dtilde2;
P1=P/.Dtilde->Dtilde1;
P2=P/.Dtilde->Dtilde2;
Ttilde = 0.1;
Result=
FindRoot[
{P1==P2,mu1==mu2},
{Dtilde1,DtildeG},
{Dtilde2,DtildeL},
MaxIterations->1000,
WorkingPrecision->16];
Link to the forum page for this post:
http://www.mathematica-users.org/webMathematica/wiki/wiki.jsp?pageName=Special:Forum_ViewTopic&pid=10101#p10101
Posted through http://www.mathematica-users.org [[postId=10101]]