#include <gsl/gsl_eigen.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_vector.h>

int main(int argc, char** argv) {
   gsl_eigen_symmv_workspace* wsp;
   gsl_matrix* in, *out;
   gsl_vector* vec;

   in = gsl_matrix_alloc(2,2);
   out = gsl_matrix_alloc(2,2);
   vec = gsl_vector_alloc(2);
   
   wsp = gsl_eigen_symmv_alloc(2);
       
   gsl_matrix_set(in, 0, 0, 1);
   gsl_matrix_set(in, 0, 1, GSL_POSINF);
   gsl_matrix_set(in, 1, 0, GSL_POSINF);
   gsl_matrix_set(in, 1, 1, 1);
	
   gsl_eigen_symmv(in, vec, out, wsp);
	
   gsl_eigen_symmv_free(wsp);
}