%Hard code Newton's Method

x=1; e1=1; r=0;
numIters=10;

T=zeros(numIters+1,3);  %Table for storage
T(1,:)=[x e1 0];

for j=1:numIters
    y=x^2;
    dy=2*x;
    x=x-y/dy;
    e2=abs(x-r);
    T(j+1,:)=[x e2 e2/e1];  %Test for linear convergence?
    e1=e2;
end

T