       Getting number out of NDSolve

• To: mathgroup at smc.vnet.net
• Subject: [mg24404] Getting number out of NDSolve
• From: colton at socrates.berkeley.edu (John S. Colton)
• Date: Tue, 18 Jul 2000 00:58:04 -0400 (EDT)
• Organization: University of California at Berkeley
• Sender: owner-wri-mathgroup at wolfram.com

```Hi, all,

I used Mathematica's function NDSolve to solve some coupled
differential equations.

Following the example in the manual, the command was something like
this:
sol = NDSolve[{ x1''[t] + eta x1'[t] + 2 Pi Sin[2 Pi x1[t]] - k
(x2[t]-x1[t]-a12[t]) == 0,
x2''[t] + eta x2'[t] + 2 Pi Sin[2 Pi x2[t]] - k (x3[t]-x2[t]-a23[t])
+ k (x2[t]-x1[t]-a12[t]) == 0,
x3''[t] + eta x3'[t] + 2 Pi Sin[2 Pi x3[t]] + k (x3[t]-x2[t]-a23[t])
== 0,
x1 == 0, x2 == 1, x3 == 2, x1' == x2' == x3' == 0},
{x1,x2,x3},{t,0,tmax}]

(a12 and a23 functions were defined earlier, as well as the eta and k
constants)

Then, to plot the solutions, I used something like this:
Plot[{Evaluate[x1[t] /. sol],Evaluate[x2[t] /. sol],Evaluate[x3[t] /.
sol]}, {t,0,tmax}]

Worked just great! Solutions looked like what I expected.

However, when attempting to get a number for the functions x1,x2, or
x3, for a particular time, all my attempts fail short.

I've tried things like N[x1[tmax]], or Evaluate[x1[tmax] /. sol, or
Evaluate[x1[t] /. sol /. t -> tmax, etc,etc. I can't get anything to
spit out a number.

I appreciate any help.

John

```

• Prev by Date: Re: How to use BestFitParameters?
• Next by Date: Re: Problem with antiderivtive of ArcSec
• Previous by thread: two variable recurrence relations
• Next by thread: Re: Need to reduce 2 lists so that only {x,y} pairs with same x remain