Метод простой итерации в общем виде

#include <iostream>
#include <cmath>
 
using namespace std;
int main()
{
    int N,M, MAX,count;;// åñëè ñêàçàòü, ÷òî ðåøåíèå ìàòðèöû, ãäå êîë-âî óðàâíåíèé íå ñîâïàäàåò ñ êîë-âîì íåèçâåñòíûõ íå èìååò ñìûñëà, òî M=N
    int MASSIV[M];//ìàññèâ, ñîäåðæ-èé êîýôôèöèåíòû ìàòðèöû "â ñòðî÷êó"
    int diag[N];
    double column[N];      // ñòîëáåö èç N ñâîáîäíûõ ÷ëåíîâ
    double x1 = 0, x2 = 0, x3 = 0, x4 = 0;          // íà÷àëüíîå ïðèáëèæåíèå
    double x1n, x2n, x3n, x4n;
    double eps = 0.0000000000001;
    
//******************************************************************************************
 
    // äàëåå öèêë, ïðîâåðÿþùèé ñõîäèìîñòü ìåòîäà
    int numb=0;
    numb=M/N;//êîë-âî ñòðîê
    for(int i=0; i<M; i+N)
    if ()
         cout << "Diagonal'noe preobladaniye ne narushaetsya" << endl;
    else 
    {
        cout << "Diagonal'noe preobladaniye narushaetsya" << endl;
    }
    count = 0;
    do 
    {
        x1n = (b1 - (a12*x2 + a13*x3 + a14*x4))/a11; // êàê ýòî çàïèñàòü â îáùåì âèäå?
        x2n = (b2 - (a21*x1 + a23*x3 + a24*x4))/a22; // êàê çàèñàòü öèêë?
        x3n = (b3 - (a31*x1 + a32*x2 + a34*x4))/a33;
        x4n = (b4 - (a41*x1 + a42*x2 + a43*x3))/a44;
        count = count + 1;
        if(fabs(x1n-x1) < eps && 
            fabs(x2n-x2) < eps && 
            fabs(x3n-x3) < eps &&
            fabs(x4n-x4) < eps) break;
        x1 = x1n;
        x2 = x2n;
        x3 = x3n;
        x4 = x4n;
           
    }while(1);  
    x1 = x1n;
    x2 = x2n;
    x3 = x3n;
    x4 = x4n;
    cout << "Kolichestvo iter. = " << count << endl << "x1 = " << x1 << endl 
         << "x2 = " << x2 << endl << "x3 = " << x3 << endl << "x4 = " << x4 << endl;       
}