Re: [Maxima-discuss] quad_qagi call inside function fails

[Edwin W. <woo...@ch...>]

thanks for the correction. that was my error.

Ted
---------------------------------------------------------------------------
-----Original Message----- 
From: Kris Katterjohn 
Sent: Saturday, April 11, 2015 2:10 PM 
To: Edwin Woollett ; maxima-discuss 
Subject: Re: [Maxima-discuss] quad_qagi call inside function fails 

On 04/11/2015 01:57 PM, Edwin Woollett wrote:
> -----------------------------------------------------
> 1.)   quad_qagi inside function  atest1 with interval (minf,0)

> (%i1) fpprintprec:8$
> (%i2) integrate(exp(6*x),x,-inf,0);
> (%o2) 1/6
> (%i3) quad_qagi(exp(6*x),x,minf,0);
> (%o3) [0.166667,2.48306918E-11,75,0]
> (%i4) atest1(kv) := 
> block([ i1L, i1],
>    i1L : quad_qagi(kv*x, x, minf, 0),
>   print(" i1L = ",i1L),  
>   if i1L[4] = 0 then i1 : i1L[1]
>     else ( print(" error i1 "),
>              return() ),
>  i1)$
> (%i5) atest1(6);

It looks like in %i3 you're doing the integral of exp(6*x) from minf to
0, but in %i5 you're doing the integral of 6*x from minf to 0.  I think
you're missing an exp() in the quad_qagi in atest1 (and also in your
atest2).

Kris Katterjohn