#include <stdio.h>
#include <stdlib.h>
#include <gsl/gsl_blas.h>

#define MAXFILELEN 80
#define TOL 10.0E-6
#define MAXITE 100

//pre-defined functions

//GC function
//Solve Ax= b 
int fGC(gsl_matrix_float *A,gsl_vector_float *x,gsl_vector_float *b,float Tol,int MaxIte);
//sub functions for FGC
void VectorSub(gsl_vector_float *a,gsl_vector_float *b,gsl_vector_float *ans);
void VectorScale(gsl_vector_float *a,double x,gsl_vector_float *ans);
void VectorAdd(gsl_vector_float *a,gsl_vector_float *b,gsl_vector_float *ans);

int myin(gsl_matrix_float **A,gsl_vector_float **b);
int myout(gsl_vector_float *x);