各位会使MATLAB的XDJM，帮小弟这个忙，毕设一个关键问题，v1=(x(i+1)-t).^2/(2*h(i))我想求v1的表达式，t是变量，怎么搞，其中x=1:1:998，i=1:1:996,h(i)=x(i+1)-x(i)，如果维数不行可以改，其实就是得出t为变量的v1的表达式，当然后边t我还要单用个矩阵赋值。
这是我编的程序，
u(i)=1/2;
l(i)=1/2;
fii(i+1)=(fi(i)-fi(i+1))/2;
i=1:1:997;
d(i)=3*(fi(i)-fi(i+1));
n=998;
e=ones(n,1);
s2=spdiags([e/2,2*e,e/2],[-1,0,1],n-1,n+1);
sf=full(s2);
di=transpose(d);
m=sf\di;

i=1:1:996
syms t
v1=(x(i+1)-t).^2/(2*h(i))
v1=v1*(-m(i))
v2=(t-x(i)).^2/(2*h(i))
v2=v2*m(i+1)
v3=fi(i)
v4=h(i)/6
v4=v4*(m(i+1)-m(i))
v(i)=v1+v2+v3+v4
a1=(x(i+1)-t(i))/h(i)
a1=a1*m(i)
a2=(t(i)-x(i))/h(i)
a2=a2*m(i+1)
a(i)=a1+a2
plot(t,v,'r',t,a,'g')
就卡壳到这了，今天就要交论文了，急啊