04-22.txt

来自「MATLABR2006a基础教程」· 文本 代码 · 共 10 行

>> A =[0,-6,-1;6,2,-16;-5,20,-10];

>> x0=[1;1;1];

>> for t = 0:.01:1
   X = [X expm(t*A)*x0]; 
end

>> plot3(X(1,:),X(2,:),X(3,:),'-o')