用MATLAB进行龚伯兹曲线和罗吉斯蒂曲线预测

电脑安装不了MATLAB，所以运行不了程序，希望有安装的好心人帮我运行一下。
分别预测t=24和26的时候的数据。麻烦把输出的结果全都复制给我，谢谢了。

龚伯兹曲线的程序：
这个是可以运行的了，最好可以画个图，谢谢。
function f=fungmtz(y)
n=length(y);
r=n/3;
L1=0;
L2=0;
L3=0;
for i=1:r
L1=L1+log(y(i));
L2=L2+log(y(r+i));
L3=L3+log(y(2*r+i));
end
b=((L3-L2)/(L2-L1))^(1/r)
a=exp((L2-L1)*((b-1)/(b*(b^r-1)^2)))
k=exp((L1-(b*(b^r-1)/(b-1))*log(a))/r)
t=input('t=')
k*a^(b^t)
end
y=[1.59 1.21 1.38 1.46 1.34 1.78 2.28 3.47 5.94 7.85 9.26 11.19 14.78 17.55 21.86 20.34 21.82 20.21 19.9 18.32 17.44];
fungmtz(y)

罗吉斯蒂曲线的程序：
还是求t=24和26，k=60，最好有图，谢谢了。
Function f=funglogl(y) %输入时间序列y，按列输入。
K=input(‘k’)%输入增长极限，即饱和值。
Y1=log((k./y)-1);
X=[[1:length(y)]’,ones(length(y),1)];
[b,bint,r,rint,stats]=regress(y1,x);%调用线性回归程序。
Stats
If stats(3)<=0.10 %判断利用线性回归方程估计的参数是否可用。
A=-b(1),m=exp(b(2)),%显示估计的模型参数。
Fprintf(‘input time’)
T=input(‘t=’)%输入预测时间点。
k/(1+m*exp(-a*t))%输入预测时间点。
Else
Fprintf(‘the model is bad’)
End
End
y=[1.59 1.21 1.38 1.46 1.34 1.78 2.28 3.47 5.94 7.85 9.26 11.19 14.78 17.55 21.86 20.34 21.82 20.21 19.9 18.32 17.44];

都没法运行。。。表示遗憾。。。

我不会~~~但还是要微笑~~~：）