This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Integrationproblem with qags
- From: Axel Hutt <hutt at wias-berlin dot de>
- To: gsl-discuss at sources dot redhat dot com
- Date: Thu, 03 Jul 2003 18:00:23 +0200
- Subject: Integrationproblem with qags
- Organization: Weierstrass-Institute for Applied Analysis and Stochastics Berlin
Dear all,
I have to solve
dV(x,t)/dt=\int_a^b r*e^(-r*|x|)*p(x,t)
for large r by spatial and temporal integration and I tried
to do the integration by gsl_integration_qags (due to the singularity
at x=0).
It turns out, that everything works fine in case of a smooth
function p(x,t). However, if
p(x,t)=V(x,t) ,
I get error messages in most of the time. My suspicion is,
that one reason is the discreteness of V(x,t), which is
unfortunate for the extrapolation steps in qags. Is this
the reason?
Is there a routine in GSL that works better?
Axel