MathGroup Archive 2006

[Date Index] [Thread Index] [Author Index]

Search the Archive

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]]



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