#include "header.h"
int main(void)
{
	gsl_vector_float *b;
	gsl_vector_float *x;
	gsl_matrix_float *A;
	float tol;        
	int maxite;       //maximal iterate
	int n;            // length of vector

	in(&A,&b);      
	n=(*b).size;
        x=gsl_vector_float_alloc(n);  
	gsl_vector_float_set_zero(x); //x=0;
	tol=TOL;                      //control the error
	maxite=MAXITE;                //control the iterate
	cg(A,x,b,tol,maxite);

	out(x);

	gsl_vector_float_free(x);
	gsl_vector_float_free(b);
	gsl_matrix_float_free(A);
	return 0;
}