Метод Монте-Карло для решение систем линейных алгебраических уравнений

void write(double* array, int n)
{
    for (int i = 0; i < n; i++)
        cout << array[i] << endl;
    cout << endl;
}
void MonteKarlo(int n, double** a, double* b)
{
    srand(time(0));
    double s,x;
    int N = 100000;
    int** v = new int* [n];
    for (int i = 0; i < n; i++)
    {
        v[i] = new int[n];
        s = 0;
        for (int j = 0; j < n; j++)
            s += abs(a[i][j]);
        s+= abs(b[i]);
        b[i] /= -s;
        for (int j = 0; j < n; j++)
        {
            if (j == i)
                a[i][j]--;
            a[i][j] /= s;
            if (a[i][j] > 0)
                v[i][j] = -1;
            else
            {
                v[i][j] = 1;
                a[i][j] *= -1;
            }
        }
        for (int j = 1; j < n; j++)
            a[i][j] += a[i][j-1];
        cout << "S= " << s << endl;
    }
    int* m = new int[n];
    double* y = new double[n];
    for (int i = 0; i < n; i++)
    {
        y[i] = 0;
        s = 0;
        for (int j = 0; j < n; j++)
            m[j] = 0;
        for (int k = 0; k < N; k++)
        {
            x = 1.*rand() / RAND_MAX;
            for (int j = 0; j < n; j++)
            {
                if (x < a[i][j])
                {
                    m[j]++;
                    j = n;
                }
            }
        }
        for (int j = 0; j < n; j++)
            y[i] += v[j][i] * b[j] * m[j];
        y[i] /= N;
    }
    write(y, n);
    for (int i = 0; i < n; i++)
        delete[] v[i];
    delete[]v;
    delete[]m;
    delete[]y;
}