Estoy tratando de escribir una función que implemente el método de Newton en dos dimensiones y aunque lo he hecho, ahora tengo que ajustar mi script para que los parámetros de entrada de mi función sean f (x) en un vector de columna, la matriz jacobiana. def(x), la conjetura inicialx0 y la tolerancia donde la funciónf(x) y su matriz jacobiana están en archivos .m separados.

Como ejemplo de un script que escribí que implementa el método de Newton, tengo:

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

Así que con esta secuencia de comandos, implementé la función y la matriz jacobiana en la secuencia de comandos real y estoy luchando para descubrir cómo puedo crear una secuencia de comandos con los parámetros de entrada necesarios.

¡Gracias!