k=input('请输入K');Y=(X-b/a)*exp(-a*k)+b/a%%%%%%%%%%%%%%%%%%%%%%灰色预测模型GM(1,1)
y(1)=(x-b/a)+b/a;for k=1:10 y(k+1)=(x-b/a)*exp^(-a*k)+b/a;endy
function Y=my_fun(k)X=1;%这里是初值,我随便假设的a=2;b=1;if k=0Y=1;elseY=(X-b/a)*exp(-a*(k-1))+b/a;end