[Scilab-users] numderivative Heart function

[Hermes]

I wish in this case that x passes as a parameter:

function val = Hxy(x,y)
    val=(1.25*y-sqrt(abs(x))).^2+x.^2-1; // switched .^ to ^ to handle
vectors
endfunction;

x=0.0;
yint1=fsolve(1,list(Hxy,x))//// ---> (25/16)*y^2-1
disp(yint1)
yint2=fsolve(-1,list(Hxy,x))
disp(yint2)
the expected value is 4/5 and all the time I receive 0.68034

What incorrect logic am I assuming?
Please, I would like to put an end to this ("Heart") !!!!
Gracias



--
Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html