%
% RT-Vorlesung
% Beispiel für numerische Probleme bei der Berechnung der Matrix-
%          exponentialfunktion exp(A)
%
% F. Svaricek, 21.05.2002
%

%
% Definition der Matrix A = T diag(-1 -37 ) inv(T)
%
T=[1 3; 2 4]
%
A = T * diag([-1 -37]) * inv(T)

%
disp([' Berechnung mit Hilfe der Matlab-Standardfunktion EXPM'])
%
EA = expm(A)

%
disp([' Berechnung mit Hilfe der Diagonalmatrix'])
%
EA1 = T * diag([exp(-1) exp(-37)]) * inv(T)

%
disp([' Berechnung mittels des Reihenansatzes (Matlab-Funktion EXPM2)'])
%
EA2 = expmdemo2(A)

%
disp([' Berechnung mit Skalierung und Quadrierung : (exp(A/2))^2'])
%
A2 = A/2

EA3 = expmdemo2(A2)

EA3 = EA3*EA3