Re: Supplying Jacobian to lsode

[vHF]

Dear John,


John W. Eaton-3 wrote:
> 
> Try
> 
>   x = lsode (address@hidden, @jac}, init, u);
> 
> instead.
> 

Brilliant -- this works.


John W. Eaton wrote:
> 
> I think there is an error in your Jacobian.  
> ...
> BTW, this system appears to be unstable, and all the interesting
> behavior seems to happen before u = 1.
> 

I should have mentioned that I submitted the smallest example 
I could think of, not the real thing. Now I am going back to 60+ 
equations, so your performance tips are also going to be very useful!

Thanks,

Marek (vHF)

-- 
View this message in context: 
http://www.nabble.com/Supplying-Jacobian-to-lsode-tp25311485p25316691.html
Sent from the Octave - General mailing list archive at Nabble.com.