Newton Raphson como función en Matlab

function x=newtonraphson(fun,dfun,x0,tol,maxit)


k=0;
corr=tol+1;


while abs(corr)>tol & k<maxit
    
    k=k+1;
    fx0=feval(fun,x0);
    dfx0=feval(dfun,x0);
    if dfx0==0
        error('la derivada se anula')
    end
    corr=fx0/dfx0;
    x=x0-corr;
    x0=x;
end