mynewton

eps=0.000 1; %précision de la solution
nrMax =100; % Le nombre max des itéraions à ne pas dépasser
iter=0;
X0=[0;6;1;4];    % Valeurs initiales, arbitraire : tentez desvaleurs différentes
x=X0
TabF=[@Fonc1,@Fonc2,@Fonc3,@Fonc4];
ImgF=[Fonc1(x);Fonc2(x);Fonc3(x);Fonc4(x)];
vnorm = norm(inv(MaJacob2(TabF,x))*ImgF);
while(vnorm>eps & iter<NrMax)
    x=x-inv(MaJacob2(TabF,x))*ImgF;
    ImgF=[Fonc1(x),Fonc2(x),Fonc3(x),Fonc4(x)];
    vnorm = norm();