function [f,f0]=Language(x,y,x0)%%你要輸出f和f0,那你的函數(shù)的輸出參數(shù)個(gè)數(shù)必須是2個(gè)
syms t;
if(length(x)==length(y))%%這里是判斷,不能用=,要用==
n=length(x);
else
disp('x和y的維數(shù)不相等');
return;
end
f=0.0;
for(i=1:n)
l=y(i);
for(j=1:i-1)
l=l*(t-x(j))/(x(i)-x(j));
end;
for(j=i+1:n)
l=l*(t-x(j))/(x(i)-x(j));
end;
f=f+l;
simplify(f);
end
f0=subs(f,'t',x0);
%%其它的不用改
>>x=0:0.5:3;
y=[ 00.47940.84150.99750.90930.59850.1411];
[f,f0]=Language(x,y,1.6)
f =
-799/3125*t*(t-1)*(t-3/2)*(t-2)*(t-5/2)*(t-3)+561/500*t*(t-1/2)*(t-3/2)*(t-2)*(t-5/2)*(t-3)-133/75*t*(t-1/2)*(t-1)*(t-2)*(t-5/2)*(t-3)+3031/2500*t*(t-1/2)*(t-1)*(t-3/2)*(t-5/2)*(t-3)-399/1250*t*(t-1/2)*(t-1)*(t-3/2)*(t-2)*(t-3)+1411/112500*t*(t-1/2)*(t-1)*(t-3/2)*(t-2)*(t-5/2)
f0 = 0.9996
%%%直接調(diào)用運(yùn)行,上面是運(yùn)行結(jié)果.