function kraftsystem()
   % Konstante Variablen belegen
   a1 = 70; a2 = 40; b = 70; F1 = 80; 
   A = [ -cosd(a1), -cosd(a2); sind(a1), -sind(a2) ];

   n = 1;  % Startindex für die Felder
   % Beginn 120 N, jeweils um 5 N verkleinern
   for( F2 = 120:-5:30 )
      % aktuelle äußere Kraft berechnen
      r = [ -F2 * sind(b); F2 * cosd(b) + F1 ];
      % aktuelle Aufgabe lösen 
      Fs = A \ r;
      % Felder Fs1, Fs2, Fe und y für den Plot belegen
      Fs1(n) = Fs(1);
      Fs2(n) = Fs(2);
      Fe(n)  = F2;
      y(n)   = 0;
      % Laufindex der Felder um 1 erhöhen
      n = n + 1;
   end

   % Grafikausgabe aller Daten
   plot( Fe, Fs2, 'b', Fe, Fs1, 'k:', Fe, y, 'k' );
   legend( 'Seilkraft Fs2', 'Seilkraft Fs1' );
   xlabel( 'Kraft F2 / N' );
   ylabel( 'Seilkräfte / N' );
end