Re: NDSolve
- To: mathgroup at smc.vnet.net
- Subject: [mg42899] Re: NDSolve
- From: "Bill Bertram" <wkb at ansto.gov.au>
- Date: Sat, 2 Aug 2003 04:12:24 -0400 (EDT)
- Organization: Australian Nuclear Science and Technology Organisation
- References: <bgctse$9l9$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
"Steffen" <nnnx at gmx.de> wrote in message news:bgctse$9l9$1 at smc.vnet.net... > Hi all, > unfortunately Mathematica is not able to solve the following: > > NDSolve[{ > > y1'[x] == -k1 y1[x] y2[x] + km1 y3[x], > > y2'[x] == -k1 y1[x] y2[x] + km1 y3[x], > > y3'[x] == k1 y1[x] y2[x] + km2 y4[x] M - (km1 + k2 M) y3[x], > > y4'[x] == -km2 M y4[x] + k2 M y3[x], > > y1[0] = 10^13, > > y2[0] = 10^16, > > y3[0] = 0, > > y4[0] = 0, > > } , {y1 , y2 , y3 , y4} , {x , 0 , 50 10^-6}] > Hi, Replace the single "=" signs by double "==" . Also remove the comma after y4[0]=0. Now, if you assign values to all of your constants you should be OK. Cheers, Bill