关于差分方程的一个问题，求助

yxy

第一次编matlab程序，出现很多问题，还是让高手来看看吧！
matlab程序如下：

function x=yxy21(x0,n,a1,a2,b,c)
x=x0;x1=a1*b*c*x0;x2=a1*b*c*x1+a2*b*(1-a1)*b*c*x0;
for k=1:n
  x(k+2)=a1*b*c*x(k+1)+a2*b*(1-a1)*b*c*x(k);
end
k=(0:20)';
y1=yxy21(100,20,0.5,0.25,0.18,10);y2=yxy21(100,20,0.5,0.25,0.19,10);
y3=yxy21(100,20,0.5,0.25,0.20,10);round([k,y1',y2',y3']);
plot(k,y1,k,y2,':',k,y3,'--');

不知为何，那个function部分总是出现各种错误，改了很久也不明所以。

jimin

function y=yxy21(a,b)
x0=1;
y=[];
for k=1:10
y=[y;a*x0+b];
x0=y(end);
end
>> y=yxy21(2,3)

y =

           5
          13
          29
          61
         125
         253
         509
        1021
        2045
        4093
我随便写了个函数,和你的差不多，你可以参考一下，然后写出你自己需要的

yxy

谢谢了