function [c] = my_interpol(x,y)

%This Function computes the Newton-Interpolation for
%Nodes (x_0,...,x_m) and values (y_0,...,y_m)
%x and y must be vectors of the same size
%c is the vector of coefficients of the Newton basis

m=max(size(y));

%Initialisation
c=y;

for k=1:m-1
    for i=m:-1:k+1
       c(i)=(c(i)-c(i-1))/(x(i)-x(i-k));
    end
end


end


