//这里是代码窗口,请将Forcal代码写在下面
i: OutVector(p:k,i)= k=FCDLen(p),printff{"\r\n"},i=0,(i
!using["XSLSF"]; //使用命名空间XSLSF
//函数定义,用于计算微分方程组中各方程右端函数值,连分式法对微分方程组积分一步函数pbs1将调用该函数f。
//t为自变量,x为函数值,dx为右端函数值(即微分方程的值)。
//该函数被调用时将用到参数A,k,这2个参数正好是拟合参数,用模块变量传递这2个参数。
f(t,x,dx::A,k)={
dx=A*x/(k+x)
};
//计算相邻两组实验数据的理论值与实验值的误差:第一组数据作为理论值的起点,用以计算出终点处的理论值,与第二组数据代表的实验值做比较。
//hf为函数f的句柄,预先获得了该句柄;a为工作数值;step为积分步数;eps为积分精度;t1,t2为积分的起点和终点;
//x1是对应t1的实验数据,即积分起始点处的函数值;x_1是对应t2的实验数据,与积分终点处的函数值做比较。
//h,i为自动变量。
t_i_j(hf,a,step,eps,t1,t2,x1,x_1:h,i)=
{
a.setra(0,x1), //工作数组赋初值
h=(t2-t1)/step,