function [output] = newton(error)
xold=2;xnew=1;
while abs(xnew-xold)> error
     xold=xnew;
     xnew=xnew/2+3/(2*xnew);
end
output=xold;