Hello all,
I have a function, f2d, defined in the following steps:
f( x, l ) := quad_qagi( t^x * exp(-t) * log(t), t, l, inf )[1];
f1d(x, l) := f(x, l )/ gamma(x+1) - f(x, 0) * gamma_incomplete(x+1, l) / gamma(x+1)^2;
f2d( s, d, l1, l2) := f1d( (s+d)/2, l1 ) * f1d( (s-d)/2, l2 ) / 4;
Now, up to this point I can get Maxima to perform the numerical evaluation, eg
f2d( 4, -2, 2, 3), numer;
gives:
.01481272382896876
If I try to integrate f2d, though, maxima balks. For example:
quad_qagi( f2d(4, d, 2, 3), d, minf, inf );
refuses to evaluate. How do I get Maxima to perform this evaluation?
Thanks,
Sean