clcclear	
x0=[110?148.5?208.5?323.5?441.1]‘;%注意這里為列向量
n=length（x0）;
lamda=x0（1:n-1）./x0（2:n）?%計(jì)算級比
range=minmax（lamda‘）?%計(jì)算級比的范圍
sita=[exp（-2/（n+1））exp（2/（n+1））]%級比可容范圍
x1=cumsum（x0）;%累加運(yùn)算
B=[-0.5*（x1（1:n-1）+x1（2:n））ones（n-11）];
Y=x0（2:n）;
u=B\Y
x=dsolve（‘Dx+a*x=b‘‘x（0）=x0‘）;
x=subs（x{‘a(chǎn)‘‘b‘‘x0‘}{u（1）u（2）x1（1）}）;
yuce1=subs（x‘t‘[0:n-1]）;
%為提高預(yù)測精度，先計(jì)算預(yù)測值，再顯示微分方程的解