#include <stdio.h>

#include <gsl/gsl_poly.h>

int

main (void)

{

  int i;

  /* coefficients of P(x) =  -1 + x^5  */

  double a[6] = { -1, 0, 0, 0, 0, 1 };  

  double z[10];

  gsl_poly_complex_workspace * w 

      = gsl_poly_complex_workspace_alloc (6);

  gsl_poly_complex_solve (a, 6, w, z);

  gsl_poly_complex_workspace_free (w);

  for (i = 0; i < 5; i++)

    {

      printf ("z%d = %+.18f %+.18f\n", 

              i, z[2*i], z[2*i+1]);

    }

  return 0;

}