Problema 1 ME

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

//Problema 1 
clc
function y=F(x)
f1=6.*x(1)-2.*cos(x(2).*x(3))-1
f2=9.*x(2)+sqrt(x(1).^(2)+ sin(x(3))+1.06)+0.9
f3=60.*x(3)+ 3.*(exp(-x(1).*x(2)))+10.*%pi-3
y=[f1;f2;f3]
endfunction
function y=J(x)
J11=6
J12=2.*x(3).*sin(x(2).*x(3))
J13=2.*x(2).*sin(x(2).*x(3))
J21=x(1)./(sqrt(x(1).^(2) + sin(x(3))+1.96))
J22=9
J23=(0.5.*cos(x(3)))./(sqrt(x(1).^(2)+sin(x(3)) + 1.96))
J31=-3.*exp(-x(1).*x(2)).*x(2)
J32=-3.*exp(-x(1).*x(2)).*x(1)
J33=60
y=[J11,J12,J13;J21,J22,J23;J31,J32,J33]
endfunction
X0 =[ 0; 0; 0 ]
w = naolinear(F, J, X0, 0.0001)
disp(w)
v = w(:,$)
disp("VETOR SOLUÇÃO OBTIDO: ")
disp(v)
disp("Substituindo a solução em F")
disp( F(v) )