求助一循环求解的问题，急

问题来源：http://www.ilovematlab.cn/thread-60022-1-1.html

问题描述：

http://hi.csdn.net/attachment/200912/19/37536_1261192846fTX3.png

Forcal求解代码：

//http://www.ilovematlab.cn/thread-60022-1-1.html
!SetRealStackMax(1000);
Rm(t,m_1:j,s,a,b,h:hVm)=
{
h=0.0001,
sys::call(2,hVm :t+h,m_1 : &a),  //计算s=Vm(t+h,m_1)
sys::call(2,hVm :t,m_1 : &b),    //计算s=Vm(t,m_1)
s=[a-b]/h,
j=0,(j<=m_1).while{
sys::call(2,hVm :t,j : &a), sys::call(2,hVm :t,m_1-j : &b),  //计算a=Vm(t,j),   b=Vm(t,m_1-j)
s=s+a*b,
j++
},
s-[1-which(m_1<1,0,1)] //最后返回这个值。看这句有没有问题：取m_1<1，而不是m_1<=1
};
//Fkappa将返回s的值。mm是模块变量，因调用Rm(t,mm-1)时mm将改变，故用a预先保存，函数返回前恢复mm的值
Fkappa(t:a,s,kappa:mm)= kappa=0, a=mm, s=t^(2*kappa)*Rm(t,mm-1), mm=a, s;
Vm(t,m:h,kappa:mm,hFkappa)=
{
h=-1, kappa=0,
if[m<=0,return(t)], if[m==1,return(1/3*h*t^3)],
which(m<=1,0,1)*Vm(t,m-1)+h*{mm=m, XSLSF::fpqg[hFkappa,0,t,1e-6]}  //XSLSF::fpqg即徐士良算法中的连分式积分法
};
main(t,m:k,s:hFkappa,hVm)=
{
hFkappa=HFor("Fkappa"), hVm=HFor("Vm"),  //预先获得函数Fkappa和Vm的句柄，函数调用时将用到
s=0,k=0,(k<=m).while{
s=s+Vm(t,k),
k++
},
s  //返回s
};
main[1,1];
main[1,2];
main[1,3];
main[3,1];
main[3,2];
main[3,3];

结果：

0.6666666666666667
0.8000500036352434
0.7460237331214614
-6.
26.4004500100092
-91.63284600184863