#include <stdio.h>
/*#include <stdlib.h>*/
#include<gsl/gsl_blas.h>
#define maxfile 80
#define qloop 10.0E-10

int CG(gsl_matrix_float *A,gsl_vector_float *b,gsl_vector_float *x,int n); //prototype for CG function