MathGroup Archive 1998

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

Search the Archive

Using NDSolve for matrix equations

  • To: mathgroup at
  • Subject: [mg13538] Using NDSolve for matrix equations
  • From: "Eric R. Bittner" <bittner at UH.EDU>
  • Date: Fri, 31 Jul 1998 04:33:46 -0400
  • Organization: university of houston
  • Sender: owner-wri-mathgroup at

I need to solve a set of coupled differential equations  using NDSolve.
The basic form is:

X''[t] == -F[X[t]] + Constraint[X[t],X'[t]]; (basically a mechanical
problem with constraints) No problem if X is a simple scaler.  But,
what if X[t] is a table?   
Still, no problem so long as I can 
expand everything and write it out explicily as eom = {
X1''[t]== -F[X1[t],X2[t],.... ] + ..., X2''[t]== -F[X1[t],X2[t], ...]  +
and so on.

The problem is that I have about 200 to 2000 coupled  equations
(depending upon the problem). 

The obvious attempt:

eom = {X''[t]==-0.1 X[t],
X'[t]=={{0,0},{0,0}} };
soln = NDSolve[eom,X[t],{t,0,10}];

NDSolve complains that these are not well formed  equations. 

What is the proper way to input these equations in matrix form into

-- Thanks, 

Eric R. Bittner                 bittner at Dept. of Chemistry       
ph: 713-743-2775 University of Houston          fax: 713-743-2709
Houston TX 77204    
 "confutatis maledictis, flammis acribus addictis..."

  • Prev by Date: Solve for DEQ
  • Next by Date: JAva & Mathematica
  • Previous by thread: Solve for DEQ
  • Next by thread: JAva & Mathematica