clc, clear all ,close all;
format short e
%*** definiere Wertebereich & Anzahl der Matrizen
M = 10;
N = 5000;

%*** setze seed und initialisiere Zufallsmatrizen
rand('seed',121086);

%*** Initialisiere Quotient
q=1; Amin=zeros(4);

tic
for k=1:N
  A = round(2*M*rand(4,4)-M);
  
  %*** falls Gauss-Seidel konvergiert
  if ...%*** to DO !
    continue;
  end
  
  %*** to Do
  rhoGS(k)=...
  
  %*** Falls Gauss-Seidel konvergent
  if rhoGS(k)<1
    %** to DO
    rhoJac(k) = ...
    %*** berechne neuen Quotienten
    %*** to DO
    
    if qtmp < q
        %*** to DO
    end
    
  end
end
toc

%*** to DO



