c=newton_interp(y,x);
    x3=newton_eval(c,y,0);
    x=[x(2); x(3) ; x3];
    y=[y(2) ; y(3) ; x3^2-3];