function [root, ea, iter] = newtraph(f,dx,xr,es,maxit)
    if nargin < 3
        error('atleast 3 input arguments are required');
    end
    if nargin < 4 || isempty(es)
        es = 0.0001;
    end
    if nargin < 5 || isempty(maxit)
        maxit = 20;
    end
    iter = 0; ea = 100;
    while iter < maxit
        xrold = xr;
        xr = xr - f(xr)/dx(xr);
        iter = iter + 1;
        if xr ~= 0
            ea = abs((xr - xrold)/xr) * 100;
        end
        if ea <= es
            break;
        end
    end
    root = xr;
end