home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #19 / NN_1992_19.iso / spool / sci / math / symbolic / 2312 < prev    next >
Encoding:
Internet Message Format  |  1992-09-03  |  1.6 KB
  1. Path: sparky!uunet!dove!rosentha
  2. From: rosentha@bldrdoc.gov (Peter Rosenthal 303-497-5844)
  3. Newsgroups: sci.math.symbolic
  4. Subject: NDSolve question
  5. Message-ID: <5363@dove.nist.gov>
  6. Date: 3 Sep 92 16:48:16 GMT
  7. Sender: news@dove.nist.gov
  8. Organization: National Institute of Standards and Technology
  9. Lines: 36
  10.  
  11. I would like to solve a set of coupled non-linear, non-stiff
  12. second order ODE's.  One of the terms in the equation will be
  13. a gaussian noise source.  The standard deviation of the 
  14. noise term is fixed by the time increment used by the solver.
  15.  
  16. I would like to use NDSolve in mathematica, but I need to fix
  17. the stepsize of the integrator.  There doesn't seem to be anyway
  18. to do this, based on what I've read in the manual.
  19.  
  20. One method I tried was to generate a noise function before calling NDSolve
  21. that consisted of an InterpolatingFunction generated from Gaussian random 
  22. data.  Then NDSolve worked fine.  The problem with this approach is
  23. that it requires a lot of extra memory to store a long random time series
  24. necessary to span the entire time span.  I would much rather generate the
  25. noise at each time step, and not have to store an entire sequence.
  26.  
  27. One approach would be to use a fixed step 4th order Runge Kutta
  28. routine that spits out an InterpolatingFunction, but I don't know
  29. if I have the time to write one.                     
  30.  
  31. Does anyone know how to set up a fixed time step version
  32. of NDSolve?
  33.  
  34. Thanks in advance
  35.  
  36. Peter Rosenthal
  37.  
  38.  
  39.  
  40.  
  41. if I have the time to code one.N Interpolating function
  42.  
  43. from data generated by a gaussian 
  44. -- 
  45. Peter Rosenthal       Email rosenthal@cmg.eeel.nist.gov
  46.  
  47.