function [p0,err,k,y]=newton2(fx,dfx,p0,delta,epsilon,max1)
%Input	- fx is the function input as a symbolic string 'fx'
%	- dfx is the derivative of f input as a symbolic string 'dfx'
%	- p0 is the initial approximation to the zero of f
%	- delta is the tolerance for p0
%	- epsilon is the tolerance for the function
%	- max1 is the maximum number of iterations
%Output	- p0 is the Newton-Raphson approximation to the zero
%	- err is the error estimate for p0
%	- k is the number of iterations
%	- y is the function value f(p0)
for k=1:max1
   x = p0;
   p1=p0-eval(fx)/eval(dfx);
   err=abs(p1-p0);
   relerr=2*err/(abs(p1)+delta);
   p0=p1;
   x=p0;
   y=eval(fx);
   X = [k,x,y,relerr];
   disp(X)
   if (err<delta)|(relerr<delta)|(abs(y)<epsilon), break, end
end