[Scilab-users] solve a system of equations in scilab

[Hermes]

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