newtonRaphson

newtonRaphson.m

% Third stab at matlab
function out = newtonRaphson(funct,old)
    
    delta = 0.001;
    y = funct(old);
    m = (funct(old+delta)-y)/delta;
    new = (m*old-y)/m;

    apprx = funct(new);
    
    % Let's check to see if it hit the needed level
    if abs((apprx - y)/2) < 0.01 || funct(apprx) == 0,
        out = new;
        return;
    end
    
    out = newtonRaphson(funct,new);