Implementacja metody gaussa-seidla

#include <iostream>

using namespace std;
#define R 4

void gaussSeidel( double ** macierz )
{
    double m = 0;
    double * X = new double[ R ];
    double * Xprev = new double[ R ];
    int iteracje = 5;
   
    for( int i = 0; i < R; i++ ) X[ i ] = 0, Xprev[ i ] = 0;
   
    while( iteracje-- ) {
       
        for( int i = 0; i < R; i++ ) {
            m = 0;
            for( int j = 0; j < i; j++ )
                 m += macierz[ i ][ j ] * X[ j ];
           
            for( int j = i + 1; j < R; j++ )
                 m -= macierz[ i ][ j ] * Xprev[ j ];
           
            X[ i ] =( macierz[ i ][ R ] - m ) / macierz[ i ][ i ];
           
        }
       
        for( int i = 0; i < R; i++ ) {
            Xprev[ i ] = X[ i ];
        }
       
        for( int k = 0; k < R; k++ ) {
            cout << "x" << k << " " << X[ k ] << endl;
        }
        cout << endl;
       
       
    }
   
    delete[] X;
    delete[] Xprev;
}

int _tmain( int argc, _TCHAR * argv[] )
{
   
    double ** macierz = NULL;
    int k = 0;
   
    macierz = new double *[ R ];
    for( int i = 0; i < R; i++ )
         macierz[ i ] = new double[ R + 1 ];
   
    double tab[] = { - 1.0, 1 / 10, - 1 / 5, 0, - 3 / 5,
        1 / 11, - 1.0, 1 / 11, - 3 / 11, - 25 / 11,
        - 1 / 5, 1 / 10, - 1.0, 1 / 10, 11 / 10,
        0, - 3 / 8, 1 / 8, - 1.0, - 15 / 8 };
   
    for( int i = 0; i < R; i++ ) {
        for( int j = 0; j < R + 1; j++ ) {
            macierz[ i ][ j ] = tab[ k ];
            k++;
        }
    }
   
    gaussSeidel( macierz );
   
    for( int i = 0; i < R; i++ )
         delete[] macierz[ i ];
   
    delete[] macierz;
    return 0;
}