%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第五章例5-5-9程序
% 對稱實矩陣矩陣指數的求法
% 西安電子科技大學出版社出版  陳懷琛編著 2007年7月
%
clear,clf,format compact
A = [ -2, 4, 1; 4, -2, -1; 1, -1, -3 ],
[p,D] = eig(A)
y=p*[exp(D(1,1)),0,0;0,exp(D(2,2)),0;0,0,exp(D(3,3))]*inv(p)