舉個例子，希望有所幫助。代碼function main()

clc; clear all;

x=[0.1 0.2 0.15 0 -0.2 0.3];

y=[0.95 0.84 0.86 1.06 1.5 0.72];

x0=[-0.2:0.01:3];

y0=lglrcz(x,y,x0);

figure; hold on; box on;

plot(x, y, 'ro');

plot(x0, y0, 'k-');function s=lglrcz(x,y,x0)

n=length(x);

m=length(x0);

for i=1:m

t=0.0;

for j=1:n

u=1.0;

for k=1:n

if k~=j

u=u*(x0(i)-x(k))/(x(j)-x(k));

end

end

t=t+u*y(j);

end

s(i)=t;

end

結果