function s=Newton(fx,dfx,x0,eps)
if nargin==3
    eps=1.0e-6;
elseif nargin<3
    return
end
x=x0;
f0=eval(fx);
f1=eval(dfx);
if (f0~=0 && f1~=0)

    while abs(x1-x0)>eps    
        x1=x0-f0/f1;
        x=x1;
        f2=eval(fx);
        f3=eval(dfx);
        x1=x0;
        x0=x0-
        
end
i=1;
f2=1;
while abs(f2)>eps
    x2=x1-f1*(x0-x1)/(f0-f1);
    x=x2;
    f2=eval(fx);
    x=[i,x0,x1,x2,f2];
    disp(x);
    x0=x1;
    f0=f1;
    x1=x2;
    f1=f2;
    i=i+1;
end