MathGroup Archive 1996

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

Search the Archive

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