[Octave-bug-tracker] [bug #32959] lsode use "1/f" instead "f"

[anonymous]

URL:
  <http://savannah.gnu.org/bugs/?32959>

                 Summary: lsode use "1/f" instead "f"
                 Project: GNU Octave
            Submitted by: None
            Submitted on: Срд 30 Мар 2011 21:28:49
                Category: Libraries
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Incorrect Result
                  Status: None
             Assigned to: None
         Originator Name: Valber
        Originator Email: address@hidden
             Open/Closed: Open
         Discussion Lock: Any
                 Release: 3.3.50
        Operating System: GNU/Linux

    _______________________________________________________

Details:

Solve Ordinary Differential Equations
     dy              
     -- = f (y, x)
     dx
I set inline function 
f=inline("y/x")

and lsode make eror

if set f=inline("x/y")

output is correct

TEST:
>>> f=inline("-x/y");
>>> lsode(f,5,x);
error: lsode: repeated convergence failures (t = -3.72945e-154perhaps bad
jacob
ian supplied or wrong choice of integration method or tolerances)
>>> f=inline("x/y");
>>> lsode(f,3,x);
>>> plot(x,ans)
>>> f=inline("-y/x")
>>> lsode(f,3,x);
>>> plot(x,ans)




    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?32959>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/