#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "matrix9.h"

main( int argc, char * argv[] )
{
  MATRIX *a, *u;
  VECTOR *b, *y, *x;
  FILE * in, * out;

	/*
	 * Open in, out
	 */

  read_aug_matrix ( in, &a, &b );

  fprintf ( out, "Input matrix\n");
  print( out, a, b );

  reduce ( a, b, &u, &y );

  fprintf (out,"Reduced matrix\n");
  print( out, u, y);
  x = back_substitute( u, y );
  fprintf (out, "Solution\n");
  print_vector (out,x);

	/*
	 * etcetera
	 */
}