diff options
author | Andreas Müller <andreas.mueller@ost.ch> | 2021-12-08 20:22:01 +0100 |
---|---|---|
committer | Andreas Müller <andreas.mueller@ost.ch> | 2021-12-08 20:22:01 +0100 |
commit | 38bf01cfbcb836fc0f95a1349bd1756c1bc789af (patch) | |
tree | 76924e8f0f3e27e365d6862dc9441fc56d13b3c1 /buch/chapters/050-differential/uebungsaufgaben/airy.cpp | |
parent | add new stuff about airy and hypergeometric functions (diff) | |
download | SeminarSpezielleFunktionen-38bf01cfbcb836fc0f95a1349bd1756c1bc789af.tar.gz SeminarSpezielleFunktionen-38bf01cfbcb836fc0f95a1349bd1756c1bc789af.zip |
complete diff problems
Diffstat (limited to '')
-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.); } |