METODO INTERPOLACION POR LAGRANGE

Un código sencillo para hacerlo es:

function y0 = lagrange_interp(x, y, x0)
y0 = 0;
n = length(x);
for j = 1 : n
t = 1;
for i = 1 : n
if i~=j
t = t * (x0-x(i))/(x(j)-x(i));
end
end
y0 = y0 + t*y(j);
end

En donde (x,y) en la entrada son los valores conocidos. x0 es el valor a interpolar.

la fuente es tomado  de :

http://www.lawebdelprogramador.com/foros/Matlab/1201500-S.O.S_Interpolacion_Lagrange_en_matlab.html