for (int i = 0; i <= N; i++) //zerowanie wektora aktualnego i wektora poprzedniego
		Wek_Aktualny[i] = 0, Wek_Poprzedni[i] = 0;
	for (int i = 0; i <= N; i++) //zapis wektora beta do wektora aktualnego
		Wek_Aktualny[i] = Wek_Beta[i];
	do
	{
		Liczba_iteracji++;
		for (int i = 0; i <= N; i++) //zapis wektora aktualnego do wektora poprzedniego
			Wek_Poprzedni[i] = Wek_Aktualny[i];
		for (int i = 0; i <= N; i++) //zerowanie wektora aktualnego
			Wek_Aktualny[i] = 0;
		for (int i = 1; i <= N; i++) //obliczanie nowego wektora aktualnego
		{
			for (int j = 1; j <= i - 1; j++)
				Wek_Aktualny[i] += Mac_Alfa[i][j] * Wek_Aktualny[j];
			for (int j = i + 1; j <= N; j++)
				Wek_Aktualny[i] += Mac_Alfa[i][j] * Wek_Poprzedni[j];
			Wek_Aktualny[i] += Wek_Beta[i];
		}
		double tmp = 0.0;
		for (int i = 1; i <= N; i++) //sprawdzanie stopu obliczeń iteracyjnych
		if (tmp < abs(Wek_Aktualny[i] - Wek_Poprzedni[i]))
			tmp = abs(Wek_Aktualny[i] - Wek_Poprzedni[i]);
		if (tmp <= E)
			break;
	} while (Liczba_iteracji < MLI);