m=0;n=0;tol=1e-12;
while (m~=3)&(n~=3)
  disp(' ')
  A=input('Bitte geben Sie eine 3x3-Matrix ein: A=');
  [m,n]=size(A);
end
disp(' ')
disp('Charakteristisches Polynom von A berechnen: Bitte Taste drücken ...')
disp(' ')
pause
chi=poly(A);
str={'' '' ''};
if chi(2)>tol
  str{1}=[' + ' num2str(chi(2)) ' * lambda^2'];
elseif chi(2)<-tol
  str{1}=[' - ' num2str(abs(chi(2))) ' * lambda^2'];
end
if chi(3)>tol
  str{2}=[' + ' num2str(chi(3)) ' * lambda'];
elseif chi(3)<-tol
  str{2}=[' - ' num2str(abs(chi(3))) ' * lambda'];
end
if chi(4)>tol
  str{3}=[' + ' num2str(chi(4))];
elseif chi(4)<-tol
  str{3}=[' - ' num2str(abs(chi(4)))];
end
disp(['chi(lambda) = lambda^3' str{1} str{2} str{3}])
disp(' ')
disp('Vektor der Eigenwerte von A berechnen: Bitte Taste drücken ...')
pause
[V,Lambda]=eig(A);
Lambda=transpose(diag(Lambda))
disp('Matrix der Eigenvektoren von A berechnen: Bitte Taste drücken ...')
pause
V
disp('Inverse von A berechnen: Bitte Taste drücken ...')
pause
Ainvers=inv(A)
disp('Testprogramm beendet')
disp(' ')