[SciPy-user] Runge-Kutta ODE vs. built-in tools

[Zane Selvans]

Anne Archibald wrote:
> On 17/03/2008, Zane Selvans <zane at ideotrope.org> wrote:
>
>>  Does anyone have a little piece of example code that uses the built in
>>  ODE solvers successfully on their system so I can test my system, and
>>  get a better idea of how to use the tools?  Maybe running on a set of
>>  equations with a known analytic solution for comparison?
> 
> There is (now) an entry in the scipy tutorial:
> http://www.scipy.org/SciPy_Tutorial#head-3daf2af2101f73b650604a8e5553301ba5d5cfc6
> 
> It's very simple, but it should be enough to get you started.

Thanks, I also found another couple of examples of varying complexity:

http://www.scipy.org/LoktaVolterraTutorial
http://www.scipy.org/SciPyPackages/Integrate

> You can also verify that your scipy is installed correctly with scipy.test().

Unfortunately, I can't actually.  The version I finally got to work came 
from the current SVN repository, and the test() framework is undergoing 
some changes.  It's currently not functional.  But I did run through the 
above examples in my interpreter, and it gives the right answers, so it 
seems I have the integrate package working.  Will play with it a bit now.

>>  I'm a grad student at the University of Colorado/Caltech working on the
>>  tectonics of icy satellites - I'm modeling the tidally induced stresses
>>  on the surface of an icy satellite that has a subsurface ocean (like
>>  Europa or Ganymede or Titan), so that I can compare those stresses to
> Fascinating! I'd love to see the results when you get it working...

Well, glad to know there's at least one person on Earth who thinks it's 
interesting :)

-- 
Zane Selvans
Amateur Human
zane at ideotrope.org
303/815-6866
PGP Key: 55E0815F
-------------- next part --------------
A non-text attachment was scrubbed...
Name: zane.vcf
Type: text/x-vcard
Size: 254 bytes
Desc: not available
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20080317/011f60ab/attachment.vcf>