1. for (int i = 0; i <= N; i++) //zerowanie wektora aktualnego i wektora poprzedniego
2.                 Wek_Aktualny[i] = 0, Wek_Poprzedni[i] = 0;
3.         for (int i = 0; i <= N; i++) //zapis wektora beta do wektora aktualnego
4.                 Wek_Aktualny[i] = Wek_Beta[i];
5.         do
6.         {
7.                 Liczba_iteracji++;
8.                 for (int i = 0; i <= N; i++) //zapis wektora aktualnego do wektora poprzedniego
9.                         Wek_Poprzedni[i] = Wek_Aktualny[i];
10.                 for (int i = 0; i <= N; i++) //zerowanie wektora aktualnego
11.                         Wek_Aktualny[i] = 0;
12.                 for (int i = 1; i <= N; i++) //obliczanie nowego wektora aktualnego
13.                 {
14.                         for (int j = 1; j <= i - 1; j++)
15.                                 Wek_Aktualny[i] += Mac_Alfa[i][j] * Wek_Aktualny[j];
16.                         for (int j = i + 1; j <= N; j++)
17.                                 Wek_Aktualny[i] += Mac_Alfa[i][j] * Wek_Poprzedni[j];
18.                         Wek_Aktualny[i] += Wek_Beta[i];
19.                 }
20.                 double tmp = 0.0;
21.                 for (int i = 1; i <= N; i++) //sprawdzanie stopu obliczeń iteracyjnych
22.                 if (tmp < abs(Wek_Aktualny[i] - Wek_Poprzedni[i]))
23.                         tmp = abs(Wek_Aktualny[i] - Wek_Poprzedni[i]);
24.                 if (tmp <= E)
25.                         break;
26.         } while (Liczba_iteracji < MLI);