function ne=Newton(fx, w, n) % initial guess w, do Newton method n times
g=diff(fx);
x=w;
for i=1:n   
    ag=eval(g);              % ag=g(x)
    af=eval(fx);             % af=f(x)
    if ag~=0 & af~=0      
       b=x-af/ag;
        x=b;
     X=[i x eval(fx)];
    % Y(i)=x-1;
     disp(X)
   end
end
   % for j=1:n-2
    %    if abs(Y(j))>0.000001
     %   Z(j)=Y(j+1)/(Y(j)^2);
     %  end 
     %end
    % Z

%function y=f(x)
%y=3*x+sin(x)-exp(x);

%function z=g(x)
%z=diff(f(x));