Mejor método de resolver sistema no lineal

syms x y z t

%Punto inicial para iterar

X01=1 ; X02=1; X03=1; X04=1;

X=[x y z t];

%Ahora tengo que poner el sistema

f(1)=;

f(2)=;

f(3)=;

f(4)=;

b=[0 0 0 0];

TOLERANCIA=0.0000000000001;

X0=[X01; X02; X03; X04]; % Es el punto inicial que se escoge para iterar;

Jf=jacobian(f,X);

t=1;

j=1;

        while (t>=TOLERANCIA) && (j<100)

            A=double(subs(Jf,X,X0));

            B=b'-subs(f,X,X0)'+A*X0;

            S=A\B;

            X0=S;

            t=norma(b-subs(f,X,X0));

            j=j+1;

        end

%X tiene que estar en columna

function J=norma(X)

J=sqrt(X*transpose(X));