diff options
Diffstat (limited to 'buch/chapters/050-differential/uebungsaufgaben/airy.cpp')
| -rw-r--r-- | buch/chapters/050-differential/uebungsaufgaben/airy.cpp | 8 |
1 files changed, 2 insertions, 6 deletions
diff --git a/buch/chapters/050-differential/uebungsaufgaben/airy.cpp b/buch/chapters/050-differential/uebungsaufgaben/airy.cpp index 20cad38..bbd5166 100644 --- a/buch/chapters/050-differential/uebungsaufgaben/airy.cpp +++ b/buch/chapters/050-differential/uebungsaufgaben/airy.cpp @@ -34,8 +34,6 @@ double h0f1(double c, double x) { k++; term *= x / ((c+k-1) * k); s += term; - // if (debug) - // fprintf(stderr, "term = %.14f\n", term); counter++; } if (debug) @@ -43,11 +41,9 @@ double h0f1(double c, double x) { return s; } -double f0(double x) { - //return gsl_sf_hyperg_0F1(2/3, x*x*x/9.); - return h0f1(2./3., x*x*x/9.); -} double f1(double x) { + // unfortunately, gsl_sf_hyperg_0F1 does not work if c<1, because + // it uses a relation to the bessel functions //return gsl_sf_hyperg_0F1(2/3, x*x*x/9.); return h0f1(2./3., x*x*x/9.); } |