function Xs = NewtonSol_Schultz(Fun,FunDer,Xest)
for i = 1:100
    Xi = Xest - Fun(Xest)/FunDer(Xest);
    plot(Xi,Fun(Xi),'rx')
    error =abs((Xi- Xest)/Xest);
    if error < 10E-6
        Xs = Xi;
        break;
    end
    Xest = Xi;
end
if i == 100
    fprintf('Solution was not obtained in %i iterations.',100)
    Xs = ('No answer');
else
    disp('Newton Method-')
    disp('Tolerance Reached')
    disp('Error:')
    disp(error);
    disp('Solution:')
    disp(Xs);
end
end