Index: tol_tests/gsl/ticket_001398/Jamroot
===================================================================
--- tol_tests/gsl/ticket_001398/Jamroot	(revision 84140e009ab1777c937a3976bdaa1b958a45deb6)
+++ tol_tests/gsl/ticket_001398/Jamroot	(revision 84140e009ab1777c937a3976bdaa1b958a45deb6)
@@ -0,0 +1,4 @@
+lib gsl ;
+lib gslcblas ;
+
+exe test_cdf_beta : test_cdf_beta.c gsl gslcblas ;
Index: tol_tests/gsl/ticket_001398/test_cdf_beta.c
===================================================================
--- tol_tests/gsl/ticket_001398/test_cdf_beta.c	(revision 84140e009ab1777c937a3976bdaa1b958a45deb6)
+++ tol_tests/gsl/ticket_001398/test_cdf_beta.c	(revision 84140e009ab1777c937a3976bdaa1b958a45deb6)
@@ -0,0 +1,36 @@
+#include <gsl/gsl_cdf.h>
+
+int main()
+{
+  int i;
+  struct cdf_beta_args {
+    double x;
+    double a;
+    double b;
+  } test_args[] =
+    {
+      {0.9060444950591773,0.2041502528011575,3.577470746807153},
+      {0.9106045018415898,0.1294490217232419,4.295225903094223},
+      {0.9306478824000806,0.159713134506609,4.357303934038536},
+      {0.9015749655663967,0.1769167491596307,3.709593420822396},
+      {0.9056045790202916,0.1844906478517543,3.606277883345159},
+      {0.9067085161805153,0.1903702202487794,3.697913210813661},
+      {0.9155602918472141,0.2029333494515557,3.52015751714325},
+      {0.9124320398550481,0.148724728854337,4.097669113723164},
+      {0.9031198702286929,0.07202908789468146,7.504563898196396},
+      {0.9078737252857536,0.1940092241212835,3.449984239856465},
+      {0.9086117066908628,0.1761127284729697,3.611606788307608},
+      {0.9105308051221073,0.2000373082584609,3.467828750535305},
+      {0.9149131434969604,0.1729368532663472,4.128703338450808},
+      {0.9363559170160443,0.1799456963706892,4.207224837486437},
+      {0.9446900135371834,0.06873371518266619,6.27892521836533},
+      {0.9999999868417623,24.95643274099755,918.6566331776214}
+    };
+
+  for ( i = 0; sizeof( test_args ) / sizeof( test_args[ 0 ] ); ++i ) {
+    double p = gsl_cdf_beta_Pinv( test_args[i].x,
+                                  test_args[i].a,
+                                  test_args[i].b );
+  }
+  return 0;
+}