Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2008

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

Search the Archive

Re: NDSolve problem

  • To: mathgroup at smc.vnet.net
  • Subject: [mg86606] Re: NDSolve problem
  • From: "Kevin J. McCann" <kjm at KevinMcCann.com>
  • Date: Sat, 15 Mar 2008 03:07:56 -0500 (EST)
  • References: <frclra$rom$1@smc.vnet.net>

You have a lot of typos or bad Mathematica to mailer (at least to mine) 
translations, e.g. MaxSteps=AE 100000. There are others.

Kevin

Payal Sheth wrote:
> Hello everyone,
> 
> I am trying to simulate a chemical reaction to be solved using simulataneous differential equation shown below
> 
> ndsolutioncombined0=NDSolve[{dtProCat[t]=8A-kcat*ProCat[t],
> dtPro[t]==kcat*ProCat[t],
> dtCat[t]==kcat*ProCat[t],
> dtProCatCompd[t]==k1*(ProCat[t]*Compd[t])-k2*ProCatCompd[t],
> dtCompd==k2*ProCatCompd[t]-k1*(ProCat[t]*Compd[t]),
> ProCat[0]=8A10,Pro[0]==0, Cat[0]==0, Compd[0]==0, ProCatCompd[0=
> ]==0},{Pro,Cat,ProCatCompd,Compd, ProCat}, {t,0,10000}, MaxSteps=AE 100000]
> 
> To this, I would like to add a restriction Pro[t]+ProCat[t]+ProCatCompd[t]=10.
> 
> Basically, at t=0, the reactant is ProCat=10, which as the reaction progresses with time is broken down into different products. However, the total reactant+product concentration shouldn't exceed the initial reactant concentration. I am unable to keep that condition under the currect setup that I have. For example, at different Compd concentration, at certain time points I have >10 of Prot+ProCat+ProCatCompd
> 
> Can anyone tell me how I could add the restriction to my differential equations?
> 
> Thanks,
> PS
> 
> 


  • Prev by Date: Re: Update 6.0.2
  • Next by Date: Re: Re: Re: Should I uninstall the old version before upgrading?
  • Previous by thread: Re: NDSolve problem
  • Next by thread: Integrate[] confuses real and imaginary parts