function p=Newton(poly,X,Y)

tol=10^(-10);
N=10000;
p=0;
p0=inf;
i=0;
while abs(p-p0)>tol & i<=N 
    p0=p;
    p=p0-fun(poly,X,Y,p0)/funprime(poly,X,Y,p0);
    i=i+1;
end