ken.m

来自「有兴趣学系matlab的朋友」· M 代码 · 共 10 行

A=[1/3,1/2;2/3,1/2];
a=[];
z(1:2,1)=[1/3;2/3];
for i=1:10
    z(1:2,i+1)=A*z(1:2,i);
    a=[a z(1:2,i+1)];
end
x=a(1,1:10)
y=a(2,1:10)