Bonjour mon ami !
Je te suggère d'utiliser l'
Optimization Toolbox, en particulier la fonction
fsolve qui permet de résoudre des systèmes d'équations non linéaires ;-)
Tape
help fsolve ou
doc fslove pour avoir plus d'informations sur cette fonction. Tu peux aussi consulter l'
aide disponible sur le site de Mathworks, c'est la même.
f=@(x)([0.4-cos(x(1))+cos(x(2))-cos(x(3))+cos(x(4)),...
0.5-cos(5*x(1))+cos(5*x(2))-cos(5*x(3))+cos(5*x(4)),...
0.5-cos(7*x(1))+cos(7*x(2))-cos(7*x(3))+cos(7*x(4)),...
0.5-cos(11*x(1))+cos(11*x(2))-cos(11*x(3))+cos(11*x(4))]);
x0=pi./[32,6,5,3];
xroot=fsolve(f,x0);
Voici le résultat obtenu.
>> xroot
xroot =
0.0969 0.5272 0.6211 0.9952
>> 180/pi*xroot
ans =
5.5491 30.2080 35.5886 57.0214
>> f(xroot)
ans =
1.0e-012 *
0.0151 0.3448 0.6006 0.8690 Je te souhaite une bien bonne journée !
-2x(1)² + 9x(2) + x(3) = 4
x(1) + 3 x(1) x(2) - 15 x(3) = 10
25 x(1) + 3 x(2)² - 2 x(2) x(3) = 12
-2x(1)² + 9x(2) + x(3) = 4
x(1) + 3 x(1) x(2) - 15 x(3) = 10
25 x(1) + 3 x(2)² - 2 x(2) x(3) = 12