[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Supplying Jacobian to lsode

From: vHF
Subject: Re: Supplying Jacobian to lsode
Date: Sun, 6 Sep 2009 04:13:19 -0700 (PDT)

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!


Marek (vHF)

View this message in context:
Sent from the Octave - General mailing list archive at

reply via email to

[Prev in Thread] Current Thread [Next in Thread]