integrate(expr1*f(expr2),x,minf,inf)
This is just one of many problems with this type of integral. I don't know how easy it is to fix it because you can get
same situation by choosing many different possible expressions for expr1 and expr2. Also I think in some cases the
final answer "may" be wrong even if you give consistent answers. I have encountered "und" sometimes which I think is
the case I tried that was wrong. I could try to reproduce that one, prove it numerically and post it if I can.
Here is another really bad case (one of many).
(%i1) assume(a>0, b>0,c>0,%sigma>0);
(eq1) [a>0,b>0,c>0,%sigma>0]
(%i2) f(x):=(%e^((-(x)^2)/(4*%sigma^2))/(2^(1/4)*%pi^(1/4)*sqrt(%sigma)))^2;
(eq2) f(x):=(%e^((-x^2)/(4*%sigma^2))/(2^(1/4)*%pi^(1/4)*sqrt(%sigma)))^2
(%i3) limit(integrate((x-5)^5*f(5*x-8), x, minf, inf), %sigma, 0);
(eq3) -1419857/15625
(%i4) integrate(f(x)*x^2, x, minf, inf)-integrate(f(x)*x, x, minf, inf)^2;
(eq4) %sigma^2
(%i5) limit(integrate((x-5)^5*f(5*x^4+x^3), x, minf, inf), %sigma, 0);
Is yx positive or negative? p;
Is log(yx) negative or zero? p;
Is log(yx) negative or zero? p;
Is log(yx) negative or zero? n;
Is 32000 * sqrt(2) * %sigma * sqrt(-log(yx)) - 27 positive, negative, or zero? p;
Is yx positive or negative? n;
Is sqrt(log(-yx)^2+%pi^2) + log(-yx) positive, negative, or zero? z;
Is sqrt(log(-yx)^2+%pi^2) - log(-yx) positive, negative, or zero? n;
Is 1728000 * %sigma * log(-yx) * sqrt(log(-yx)^2+%pi^2) + 3456000 * %sigma * log(-yx)^2 positive, negative, or zero?
p;
Is 2097152000000000000 * %sigma^4 * log(-yx)^6 + sqrt(log(-yx)^2+%pi^2) *
(-2097152000000000000*%sigma^4*log(-yx)^5-1492992000000*%sigma^2*log(-yx)^4+4194304000000000000*%pi^2*%sigma^4*log(-yx)^3+6291456000000000000*%pi^4*%sigma^4*log(-yx)+1492992000000*%pi^4*%sigma^2)
+ 1492992000000 * %sigma^2 * log(-yx)^5 +
(531441-3145728000000000000*%pi^2*%sigma^4) * log(-yx)^4 - 5971968000000 * %pi^2 * %sigma^2 * log(-yx)^3 +
(12582912000000000000*%pi^4*%sigma^4+1062882*%pi^2) * log(-yx)^2 + 4478976000000 * %pi^4 * %sigma^2 * log(-yx) +
1048576000000000000 * %pi^6 * %sigma^4 + 531441 * %pi^4
positive, negative, or zero?
Rich
----- Original Message -----
From: "Robert Dodier" <robert.dodier at gmail.com>
To: "Richard Hennessy" <rvh2007 at comcast.net>
Cc: "Maxima List" <maxima at math.utexas.edu>
Sent: Thursday, October 09, 2008 10:43 AM
Subject: Re: [Maxima] Gausian bug maybe?
On 10/9/08, Richard Hennessy <rvh2007 at comcast.net> wrote:
> I think this is right above but when the exponent is n you get an
> unending sequence of questions as long as you don't
> answer no to any of them.
Well, this is a pretty good example. Please make a bug report
if you have the time. I don't know when we'll fix it but let's
not lose track of it. For this particular example it may be relatively
easy to make progress. I'll try to remember to look at it.
Robert Dodier