matlab搞得我头大了！！！！

[中少]

项目已经结束了。虽然自己因为各种原因没有做到底。自己还是学到了很多的东西。对于matlab的热爱也增加了 不少。记得有个师兄说过，软件的熟悉需要有人带着你，这样的话一两个月就可以熟悉了。因为我是是深刻体会了各自为战的苦楚啊！！大家戒之~~~

[yangmp]

function f=language(x,y,x0)
syms t;
if length(x)==length(y)
n=length(x);
else
disp('x和y的维数不相等');
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);

if i==n
if nargin==3
f=subs(f,'t',x0);
else
f=collect(f);
f=vpa(f,6);
end
end
end


function f=newton(x,y,x0)
syms t;
if length(x)==length(y)
n=length(x);
c(1:n)=0.0;
else
disp('x和y的维数不相等！');
return;
end

f=y(1);
y1=0;
l=1;

for i=1:n-1
for j=i+1:n
y1(j)=(y(j)-y(i))/(x(j)-x(i));
end
c(i)=y1(i+1);
l=l*(t-x(i));
f=f+c(i)*l;
simplify(f);
y=y1;

if i==n-1
if nargin==3
f=subs(f,'t',x0);
else
f=collect(f);
f=vpa(f,6);
end
end
end