Menu

#47 Problem with numeric solve with erf

None
closed
5
2015-04-17
2013-04-03
No

num_solve does not work with erf for me, eg num_solve(erf(x)=0.9,x=1)
gives error during function evaluation (e.g. singularity)

Discussion

  • Rainer Schöpf

    Rainer Schöpf - 2013-04-03

    The following works:

    1: load specfn;
    
    2: load numeric;
    
    3: on rounded;
    
    4: num_solve(erf(x)=0.9,x=1);
    
    {x=1.16308715273}
    

    It should work without "on rounded", though.

    I will investigate.

     
  • Rainer Schöpf

    Rainer Schöpf - 2013-04-03
    • assigned_to: Rainer Schöpf
    • milestone: -->
     
  • Rainer Schöpf

    Rainer Schöpf - 2013-04-03
    • status: open --> accepted
     
  • Rainer Schöpf

    Rainer Schöpf - 2013-04-06

    Your example works now if you load the specfn package, ie.

    load numeric,specfn;
    
    num_solve(erf(x)=0.9,x=1);
    
     
  • Rainer Schöpf

    Rainer Schöpf - 2013-04-06
    • labels: num_solve, erf --> num_solve, erf, specfn
     
  • Rainer Schöpf

    Rainer Schöpf - 2013-04-09

    I've just checked in a fix for your problem: specfn is now loaded automatically when erf, ei, si, or ci are evaluated numerically.

    Please try your original example again.

     
  • Rainer Schöpf

    Rainer Schöpf - 2013-04-10
    • status: accepted --> closed
     
  • Rainer Schöpf

    Rainer Schöpf - 2015-04-17
    • labels: num_solve, erf, specfn --> num_solve, erf, specfn, fixed
     

Log in to post a comment.

MongoDB Logo MongoDB