A = [0    -6    -1;
     6     2   -16;
    -5    20   -10];                 %矩陣A
x0=[1;1;1];                           %初值x(0)
t=0:0.01:1;                           %時間t
xt=[];                               %計算結果x(t)
for i=1:length(t),                      %計算每個時間點上的x(t)
    xt(i,:)=expm(t(i)*A)*x0;
end;
plot3(xt(:,1),xt(:,2),xt(:,3),'-o')           %畫三維圖顯示x(t)
grid on;                             %給圖添加輔助網格
xlabel('xt(:,1)');                       %給圖的坐標軸x命名
ylabel('xt(:,2)');                       %給圖的坐標軸y命名
zlabel('xt(:,3)');                       %給圖的坐標軸z命名
title('矩陣指數函數結果');             %給圖添加標題