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] ====