function FuncyNewt(nIter,Xini)
 
Xold = Xini;
 
for i = 1:nIter
    Xnew = 100*(1/i -Xold);
    Xold = Xnew;
    X(i) = Xnew;
end
Xfinal = Xnew;
X'
% print using disp
mystr = ['after ', num2str(nIter),' iterations, approximate square root of 5 is ', num2str(Xfinal)];
disp(mystr);


function PS1func(nIter,Xini)
% nIter = 3;
Xold = Xini;
 
for i=1:nIter
    Xnew = 100*(1/i - Xold);
    Xold = Xnew;
    [i Xnew]
end
Xfinal =Xnew;
% print using disp
str = ['after ', num2str(nIter),' iterations, approximate solution is ', num2str(Xfinal)];
disp(str);

%This is version 1 script before we made it into a function
Xold = 1.5;
nIter = 2;
% X0 is initial guess
% does 4 iterations, returns X4
% Xnp1 = Xn - f(Xn)/f'(Xn)
% F(X) = X^2 -5
for i = 1:nIter
    Xnew = Xold - (Xold^2 -5)/(2*Xold);
    Xold = Xnew;
end
Xfinal = Xnew;
% print using disp
mystr = ['after ', num2str(nIter),' iterations, approximate square root of 5 is ', num2str(Xfinal)];
disp(mystr);