aboutsummaryrefslogtreecommitdiffstats
path: root/buch/chapters/050-differential/uebungsaufgaben/airy.cpp
diff options
context:
space:
mode:
authorAndreas Müller <andreas.mueller@ost.ch>2021-12-08 20:22:01 +0100
committerAndreas Müller <andreas.mueller@ost.ch>2021-12-08 20:22:01 +0100
commit38bf01cfbcb836fc0f95a1349bd1756c1bc789af (patch)
tree76924e8f0f3e27e365d6862dc9441fc56d13b3c1 /buch/chapters/050-differential/uebungsaufgaben/airy.cpp
parentadd new stuff about airy and hypergeometric functions (diff)
downloadSeminarSpezielleFunktionen-38bf01cfbcb836fc0f95a1349bd1756c1bc789af.tar.gz
SeminarSpezielleFunktionen-38bf01cfbcb836fc0f95a1349bd1756c1bc789af.zip
complete diff problems
Diffstat (limited to 'buch/chapters/050-differential/uebungsaufgaben/airy.cpp')
-rw-r--r--buch/chapters/050-differential/uebungsaufgaben/airy.cpp8
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.);
}