function y = Newton_solve(f,a,x0)
n = 1000;
delta = 0.001;
y = x0;
for i = 1 : n
    d = (f(y+delta) - f(y))/delta;
    y = y - (f(y)-a)/ d;
end
end