Re: Beginner--Help on using FindRoot to solve the system of equations

• To: mathgroup at smc.vnet.net
• Subject: [mg66221] Re: Beginner--Help on using FindRoot to solve the system of equations
• From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
• Date: Thu, 4 May 2006 05:21:12 -0400 (EDT)
• References: <e39kql\$cud\$1@smc.vnet.net>
• Sender: owner-wri-mathgroup at wolfram.com

Hi,

a) FindRoot[] needs the start values as
FindRoot[eqns,{{var1,startVar1},{var2,startVar2},..}]
and not as you type
b) what is Dr ? you have not set it and you don't ask FindRoot[]
to search for it.

Regards
Jens

xiaochu at che.utexas.edu wrote:
> This is a code to solve vapor-liquid equilibrium by van der Waals Equation of State.
>
>
>
> (* 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]]
>
>

• Prev by Date: Re: Beginner--Help on using FindRoot to solve the system of equations
• Next by Date: Re: Faster Random Walk Simulation ?!?
• Previous by thread: Re: Beginner--Help on using FindRoot to solve the system of equations
• Next by thread: Re: Beginner--Help on using FindRoot to solve the system of equations