matlab如何求非线性差分方程？求解答

y1(1)=0;
y2(1)=2;
y3(1)=9;
n=100;
for k=1:n-1
d=@(y)[0.9*y1(k)+0.1*y2(k);0.3*y1(k)-0.01*y1(k)*y3(k)-1.01*y2(k);-0.0267*y3(k)+0.01*y1(k)*y2(k)+y3(k)];
p=fsolve(d,[0,0,0],optimset('Display','off'));
y(k+1)=subs(p);
end
y
plot(y);
******************************************************************
求的是下面的这个，求指导，到底怎么求？
y1(k+1)=0.9*y1(k)+0.1*y2(k)
y2(k+1)=0.3*y1(k)-0.01*y1(k)*y3(k)-1.01*y2(k)
y3(k+1)=-0.0267*y3(k)+0.01*y1(k)*y2(k)+y3(k)
---贴吧极速版 For UWP