% Fluss von \dot(x)=Ax für verschiedene EW-Konstellationen
% by Dr. Ute Feldmann
% Das Programm ist 'gepaust' --> immer mit 'Enter' weiter

for i=[-2,2]
    close all;T=randn(2,2);A=T*[1 0; 0 i]*inv(T);PhasePortraitU(A);
    title('Zwei reelle Eigenwerte'),pause
end

close all;
%T=randn(2,2);A=-T*[1 1; 0 1]*inv(T);PhasePortraitU(A);pause
A=-[1 1; 0 1];PhasePortraitU(A);
title('Doppelter Eigenwert, geometr. Vielf.=1'),pause

for i=-1:1
    close all;T=randn(2,2);A=T*[i 2;-2 i]*inv(T);PhasePortraitU(A);
    title('Komplexe Eigenwerte'),pause
end

for i=1:4
    close all;A=randn(2,2);PhasePortraitU(A);
    title('Überraschung!'),pause
end
