Im solving a system equation within a for loop, but it never ends But its ok to evaluate f1(n) outside the loop.... why?
where eqa is
in the other hand this works
so.. whats the problem with ode_solver?
It seems like this is just ridiculously slow, but I'm not sure why. I tried tweaking your code in a variety of ways and nothing helped. To try to figure out what was going on I wrote the following Runge-Kutta 4th order code in Cython; its pretty hackish but it makes it clear that the solution isn't doing anything very strange. The version below is hard-coded for your ll=2.
and then you can get a plot with:
Anyway for me it seems easier to implement such things directly than to wrestle with the GSL wrapping.
posted Aug 27 '10mhampton
359 ● 2 ● 6 ● 15
Asked: Aug 22 '10
Seen: 123 times
Last updated: Aug 27 '10
powered by ASKBOT version 0.7.22
Copyright Sage, 2010. Some rights reserved under creative commons license.