Eu estou tentando escrever uma função que implementa o método de Newton em duas dimensões e enquanto eu tenho feito isso, eu tenho que ajustar meu script para que os parâmetros de entrada da minha função sejam f (x) em um vetor coluna, a matriz Jacobiana dof(x), o palpite inicialx0 e a tolerância onde a funçãof(x) e sua matriz jacobiana está em arquivos .m separados.

Como exemplo de um script que escrevi que implementa o método de Newton, tenho:

n=0;            %initialize iteration counter     
eps=1;          %initialize error     
x=[1;1];        %set starting value

%Computation loop     
while eps>1e-10&n<100 
    g=[x(1)^2+x(2)^3-1;x(1)^4-x(2)^4+x(1)*x(2)];         %g(x)      
    eps=abs(g(1))+abs(g(2));                             %error     
    Jg=[2*x(1),3*x(2)^2;4*x(1)^3+x(2),-4*x(2)^3+x(1)];   %Jacobian     
    y=x-Jg\g;                                            %iterate     
    x=y;                                                 %update x     
    n=n+1;                                               %counter+1     
end 

n,x,eps       %display end values

Então, com esse script, eu tinha implementado a função e a matriz Jacobiana no script atual e estou lutando para descobrir como criar um script com os parâmetros de entrada necessários.

Obrigado!