1 files changed, 26 insertions, 0 deletions
diff --git a/gsl-1.9/doc/examples/polyroots.c b/gsl-1.9/doc/examples/polyroots.c
new file mode 100644
index 0000000..6debbd5
--- /dev/null
+++ b/gsl-1.9/doc/examples/polyroots.c
@@ -0,0 +1,26 @@
+#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;
+}