Метод Гаусса. Не работает проверка

#include <locale>
#include <cmath>
#include <malloc.h>
#include <iostream>
using namespace std;
Значения х находит верно, значения q[i] тоже.
Проверка заключается в том, чтобы q[i]*x-b было равно нулю. А получаются какие-то сумасшедшие числа
 
const int n=4;
int i, j, m, k;
double aa, bb;
 
double a[n][n+1]={{  3,   4, -9,   5, -14},   //матрица коэффициентов, где последний столбец - правая часть системы
{-15, -12,  50, -16, 44},
{-27, -36,  73,   8, 142},
{  9,  12, -10, -16, -76}};
 
int main()
{
    setlocale(LC_ALL, "Russian");
    double* x;
    x = (double *)malloc(n* sizeof(double));
 
    for (k = 0; k<n; k++) //Поиск максимального элемента в первом столбце
    {
        aa = abs(a[k][k]);
        i = k;
        for(m = k+1; m<n; m++)
            if(abs(a[m][k])>aa)
            {
                i = m;
                aa = abs(a[m][k]);
            }
 
            if (aa == 0)   //проверка на нулевой элемент
            {
                cout<<"Система не имеет решений"<<endl;
            }
 
            if (i != k)  //  перестановка i-ой строки, содержащей главный элемент k-ой строки
            {
                for (j=k; j < n+1; j++)
                {
                    bb = a[k][j];
                    a[k][j] = a[i][j];
                    a[i][j] = bb;
                }
            }
            aa = a[k][k];//преобразование k-ой строки (Вычисление масштабирующих множителей)
            a[k][k] = 1;
            for (j=k+1;j<n+1;j++)
                a[k][j] = a[k][j]/aa;
            for (i = k+1; i < n; i++)//преобразование строк с помощью k-ой строки
            {
                bb = a[i][k];
                a[i][k] = 0;
                if (bb!=0)
                    for (j=k+1; j< n+1; j++)
                        a[i][j]=a[i][j]-bb*a[k][j];
            }
    }
 
    for (i=n-1; i>=0; i--)   //Нахождение решений СЛАУ
    {
        x[i] = 0;
        aa = a[i][n];
        for (j = n; j>i; j--)
            aa = aa-a[i][j]*x[j];
        x[i] = aa;
    }
 
    cout<<"Решение системы:"<<endl;  //вывод решений
    for (i = 0; i< n; i++)
    {
        cout<<"x["<<i<<"]="<<x[i];
        cout<<endl;
    }
 double b_1[n] = {-14,44,142,-76};
 double a_1[n][n]={{  3,   4, -9,   5},
                   {-15, -12,  50, -16},
                   {-27, -36,  73,   8},
                   { 9,  12, -10, -16}};
double q[n];
    for (int i = 0; i < n; i++)            //проверка
      for (int j = 0; j < n; j++)
     q[i]+= a_1[i][j]*x[j];     
    for (int i=0;i<n;i++)
   cout<<q[i]<<endl;
   cout<<endl;
   for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
        {
   cout<<q[i]-b_1[i]<<endl;
    }
 
  return 0;
}