Metoda Seidla

int k, i;
double sum, sum1;
for( k = 0; k < MLI; k++ )
for( i = 0; i < n; i++ )
{
    if( k == 0 ) x[ i ] = wektorbeta[ i ];
   
    if( i == 0 && k != 0 )
    {
        for( int j = 1; j < n; j++ )
        {
            x_nowy[ 0 ] =( macierzalfa[ 0 ][ j ] * x[ j ] ) + wektorbeta[ 0 ];
        }
    }
    if( k != 0 && i != 0 )
    {
        for( int j = 0; j < i - 1; j++ ) sum = sum + macierzalfa[ i ][ j ] * x_nowy[ j ];
       
        for( int j = i + 1; j < n; j++ ) sum1 = sum1 + macierzalfa[ i ][ j ] * x[ j ];
       
        x_nowy[ i ] = sum + sum1 + wektorbeta[ i ];
    }
   
}

}