求用MATLAB求解分数阶累加灰色模型（已有参考代码，但执行报错）

我是数学系本科生，所以计算机编程不是很好，我的毕业论文需要如题，已有代码如下

unction [MAPE,X0F] = CFEGM(r,X0,Time_X0,Time_X0F,nf)
n=numel(X0);
X1(1,1)=X0(1);
for i=2:n
X1(i,1)=X1(i-1,1)+X0(i,1)*(Time_X0(i)-Time_X0(i-1))/(Time_X0(i)^(ceil(r)-r));
end
%最小二乘法计算系数
for i=1:n-1
Z1(i)=0.5*(X1(i)+X1(i+1));
B(i,1)=-Z1(i);
B(i,2)=1;
Y(i,1)=(X1(i+1)-X1(i))/(Time_X0(i+1)-Time_X0(i));
end
P=inv(B'*B)*B'*Y;
a=P(1);
b=P(2);
%计算预测值
X1F(1,1)=X0(1);
for i=2:n+nf
X1F(i,1)=(X1F(1)-b/a)*exp(-a*(Time_X0F(i)-Time_X0F(1)))+b/a;
end
%还原
X0F(1,1)=X1F(1);
for i=2:n+nf
X0F(i,1)=Time_X0F(i)^(ceil(r)-r)*(X1F(i)-X1F(i-1))/(Time_X0F(i)-Time_X0F(i-1));
end
%计算误差
for i=2:n
e(i)=abs(X0F(i)-X0(i))/X0(i);
end
MAPE=mean(e(2:n));

我运行时会显示unction未定义，没办法运算，可以有人帮忙讲一下代码，顺便把代码补全
附需要解决的数学问题：

将unction改为function，因为函数的m文件夹的开头是function，有用望采纳