From 38bf01cfbcb836fc0f95a1349bd1756c1bc789af Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Wed, 8 Dec 2021 20:22:01 +0100 Subject: complete diff problems --- buch/chapters/050-differential/uebungsaufgaben/airy.cpp | 8 ++------ 1 file changed, 2 insertions(+), 6 deletions(-) (limited to 'buch/chapters/050-differential/uebungsaufgaben/airy.cpp') 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.); } -- cgit v1.2.1