I am trying to reproduce in Scilab in method of Draghilev. The document is developed from Maple. For the previous system develop the following script:
/function Sys=Q(x,z) Sys(1)=(x(1)-2)^2+(x(2)+2)^2+z^2-9; Sys(2)=x(1)^6+x(2)^6+z^6-12; endfunction z=0; [j,v,info]=fsolve([2;-0.1],list(Q,z)) disp(j,v,info) / 1. 0. 1.776D-15 1.496447 *0.9574371* The results do not coincide 100%. The result that differs is highlighted in black. The expected result is as follows: <http://mailinglists.scilab.org/file/t497622/InitialPonit_Dragilev_sMethod.jpg> It will be possible to graph the function to guide me to select the initial values. something like this: function y = fct1(x) y = sin(x) - .25*x; endfunction fplot2d([-5:.1:5], fct1) show_window() x0 = [-5:5]; [x,v,info]=fsolve(x0, fct1); fplot2d(x', fct1,-1) disp(x,v) Gracias Hermes -- Sent from: http://mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html _______________________________________________ users mailing list [email protected] http://lists.scilab.org/mailman/listinfo/users
