#include <stdio.h>
#include <stdlib.h>
#include <assert.h>
#include <limits.h>
#include <math.h>
#include <complex.h>
#include <string.h>
#include <gsl/gsl_poly.h>

#include "../numeric/SU21-numeric.h"
//#include "../numeric/SU21-numeric-base.c"


int main(int argc,char *(argv[])) {
  FLOAT_t coeff[7]={1,-21,175,-735,1624,-1764,720};
  FLOAT_t sol[12],*solR,*solI;
  int i;

  gsl_poly_complex_workspace *ws
    = gsl_poly_complex_workspace_alloc(7);

  gsl_poly_complex_solve(coeff,7,ws,sol);

  gsl_poly_complex_workspace_free(ws);

  solR=&(sol[0]);
  solI=&(sol[1]);
  for(i=0;i<6;i++) {
    printf(FOSpCnv FOPlCnv"*i\n",*solR,*solI);
    solR+=2; solI+=2;
  }
}