function [ ev_y ] = lagrange_interpol( ev_x, x, y )

ev_y=zeros(size(ev_x));
for j=1:size(ev_x,2)
    for i=1:size(x,2)
        ev_y(1,j)=ev_y(1,j)+y(1,i)*prod(ev_x(1,j)-x(1,[1:i-1,(i+1):end]))./prod(x(1,i)-x(1,[1:i-1,(i+1):end]));
    end
end
