NDSolve

• To: mathgroup at smc.vnet.net
• Subject: [mg3688] NDSolve
• From: FERRUCCIO Renzoni <renzoni at fexphds04.tu-graz.ac.at>
• Date: Fri, 5 Apr 1996 02:54:11 -0500
• Sender: owner-wri-mathgroup at wolfram.com

```Dear Mma Users,

in a program I use several times the following:

o = -4.;

ls  = Table[ 0,{j,1,81}];

Do[
diff = NDSolve[ eq,var,{t,0.,tf},ms];
ls[[i]] = {o,NIntegrate[(p/. diff)[[1]],{t,ti,tf},
MinRecursion -> 5,MaxRecursion -> 12]};
qu[o,.5]  = ( p /. t -> .5 /. diff )[[1]] ;
qu[o,3.5] = ( p /. t -> 3.5 /. diff )[[1]];
qu[o,6.5] = ( p /. t -> 6.5 /. diff )[[1]];
qu[o,9.5] = ( p /. t -> 9.5 /. diff )[[1]];
qu[o,12.5] = ( p /. t -> 12.5 /. diff )[[1]];
qu[o,15.5] = ( p /. t -> 15.5 /. diff )[[1]];
ClearAll[diff];
o = o + 0.1,
{i,1,81} ];

So I would like write it as a routine at the beginning of the program (I
don't like packages because I have already lot of files!).

How can I write in compact way such a routine that has as input
"eq,var,tf,p" and as output ls and qu ?
Thanks a lot.

F. Renzoni

renzoni at fexphds04.tu-graz.ac.at

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

```

• Prev by Date: HELP!
• Next by Date: FastBinaryFiles for OS/2
• Previous by thread: Re: HELP!
• Next by thread: FastBinaryFiles for OS/2