#include <stdio.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_poly.h>

int
main()
{
  double PolyParm[6];
  double z[10];
  int it;

  gsl_poly_complex_workspace *w = gsl_poly_complex_workspace_alloc(6);

  PolyParm[0] = 0;
  PolyParm[1] = 1792548; /* 1792549 */
  PolyParm[2] = 324810;
  PolyParm[3] = -288348;
  PolyParm[4] = -43875;
  PolyParm[5] = 20715;

  gsl_poly_complex_solve(PolyParm,6,w,z);

  for (it=0;it<5;it++)
    printf("z%d = %.18f + %.18f I\n",it,z[2*it],z[2*it+1]);

  gsl_poly_complex_workspace_free(w);

  return 0;
}