differential equations

• To: mathgroup at smc.vnet.net
• Subject: [mg4522] differential equations
• From: N Kinoshita <nk204 at cam.ac.uk>
• Date: Wed, 7 Aug 1996 04:17:32 -0400
• Sender: owner-wri-mathgroup at wolfram.com

```Hi, I am a research student of economics. I want to solve
a system of differential equations with Mathematica.

I wrote the follwoing command based on the example given
in p.698 of 'Mathematica: A system for doing mathematics
by computer'

------------------------------------------------------------
ks=0.380884;cs=0.0282254;ws=0.846308;
k0=ks/10;
cprimess=0.00778013/0.257126;
wprimess=0.966347/0.257126;

NDSolve[
{c'[k]==
If[k==ks,cprimess,
(-0.012+0.952757k/w[k]^.7+0.11938/w[k]^.6
+2.6603c[k]/w[k]^.1-0.635171w[k])
/(0.875c[k]+0.952757k^2/w[k]^.7+1.04457k/w[k]^.6
+2.6603c[k]k/w[k]^.1-0.635171k*w[k]^.3-0.447674w[k]^.4)
],
w'[k]==
If[k==ks,wprimess,
(1.1938w[k]^.4+6.65074c[k]w[k]^.9-1.58793w[k]^1.3)
/(0.875c[k]+0.952757k^2/w[k]^.7+1.04457k/w[k]^.6
+2.6603c[k]k/w[k]^.1-0.635171k*w[k]^.3-0.447674w[k]^.4)
],
c[ks]==cs, w[ks]==ws},
{c,w},{k,ks,k0}]
-----------------------------------------------------------

Then I got the error message saying:
NDSolve::deql:
True is not a valid equation or list of valid equations.

I was able to solve a single differential equation before,
but do not know how to cope with my problem now.

I would like to get some comments.

==== [MESSAGE SEPARATOR] ====

```

• Prev by Date: Looking for experiences integrating Mathematica with other apps
• Next by Date: Re: Pairings
• Previous by thread: Looking for experiences integrating Mathematica with other apps
• Next by thread: Re: differential equations