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

#define TOL 10.0E-6
#define MAXITE 128


int ifsym(gsl_matrix *A);
void vsub(gsl_vector *a,gsl_vector *b,gsl_vector *ans);
int QREvals(gsl_matrix *A,gsl_vector *e, int m, float Tol);