function xc = newton_tdelillo(f,fprime,x0,tol)
% Newton's method
n=4;
x(1)=x0;
for j=1:n
   x(j+1)=x(j) - f(x(j))/fprime(x(j));
end
%xc=x(n+1);
xc = x;