Matlab - necesitosoluciondeecuacionesnolinealespornewton

clear all, clc

syms x

f  = input('Función              :','s');

x0 = input('Aproximación inicial :');

tol= input('Tolerancia           :');

disp('Datos de salida')

df = diff(sym(f));

f  = inline(char(f));

df = inline(df);

x  = x0;

e  = abs(f(x));

i  = 0;

fprintf('%s\t\t\t%s\t%s\n','it.','x','f(x)');

while (e>=tol)

    fprintf('%d\t%f\t%f\n',i,x,f(x))

    nwr = x - (f(x)/df(x));

    e  = abs(f(x));

    x  = nwr;

    i  = i+1;

end