fsolve — résout un système d'équations non-linéaires
[x [,v [,info]]]=fsolve(x0,fct [,fjac] [,tol])
vecteur réel (vecteur initial).
fonction externe (fonction Scilab ou chaîne de caractères ou liste).
fonction externe (fonction Scilab ou chaîne de caractères ou liste).
scalaire. Tolérance : l'algorithme s'arrête lorsque l'erreur
relative estimée entre x et la solution est inférieure à tol
(tol=1.d-10 par défaut).
vecteur réel (solution estimée).
vecteur réel (valeur de la fonction en x).
indicateur d'arrêt de l'algorithme
paramètres d'entrée incorrects
l'erreur relative estimée entre x et la solution est inférieure à tol
nombre d'appels à fct
tol est trop petit. L'algorithme stationne.
les itérations n'améliorent pas assez la solution.
Résout un système de n équations non linéaires à n inconnues avec un algorithme du type méthode hybride de Powell. La Jacobienne peut être fournie (c'est fortement conseillé).
0 = fct(x)
fct est une fonction "externe" Cette fonction
renvoie v=fct(x) pour x
donné.
La séquence d'appel de fct est :
[v]=fct(x).
Si fct est une chaîne de caractères, elle désigne
le nom d'une subroutine Fortran ou C liée dynamiquement à Scilab, avec
comme liste d'appel Fortran :
fct(n,x,v,iflag) integer n,iflag double precision x(n),v(n)
ou alors pour une routine C :
fct(int *n, double x[],double v[],int *iflag)
(voir link).
jac est une fonction "externe" aussi. Cette
fonction renvoie v=d(fct)/dx (x) pour
x donné.
La séquence d'appel de jac est :
[v]=jac(x).
Si jac est une chaîne de caractères, elle désigne
le nom d'une subroutine Fortran ou C liée dynamiquement à Scilab. Les
séquences d'appel sont les mêmes que pour fct (attention v doit être un
tableau n x n).
// un exemple simple
a=[1,7;2,8];b=[10;11];
deff('[y]=fsol1(x)','y=a*x+b');
deff('[y]=fsolj1(x)','y=a');
[xres]=fsolve([100;100],fsol1);
a*xres+b
[xres]=fsolve([100;100],fsol1,fsolj1);
a*xres+b
// voir SCI/modules/optimization/sci_gateway/fortran/Ex-fsolve.f
[xres]=fsolve([100;100],'fsol1','fsolj1',1.e-7);
a*xres+b