#include <stdio.h>
#include <math.h>
#define MAXSIZE 50
#define EPS 1.0e-6
int solve(double A[], double b[], int nmax, int n, double x[]);
int decomp(double A[], int nmax, int n);
int residual(double A[][], double b[], double x[], int n);
int main()
{
int i, j, nmax = MAXSIZE;
int n = 3, chkerr = 0;
double A[MAXSIZE][MAXSIZE], b[MAXSIZE], x[MAXSIZE];
A[0][0] = 1.; A[0][1] = -3.; A[0][2] = 1.; A[0][3] = 7.;
A[1][0] = 1.; A[1][1] = 2.; A[1][2] = 5.; A[1][3] = 1.;
A[2][0] = 4.; A[2][1] = 1.; A[2][2] = -2.; A[2][3] = 0.;
A[3][0] = 1.; A[3][1] = 6.; A[3][2] = 2.; A[3][3] = 3.;
