MathAction #198 integrate(sin(x**2),x) is not handled

Topics

FrontPage
- IssueTracker
  - #198 integrate(sin(x**2),x) is not handled <-- You are here.

Axiom doesn't seem to do the integral of sin(x^2), but both Maxima and Mathematica (per http://integrals.wolfram.com/ anyway) produce answers. Confirmed by Martin Rubey, and uploaded at his request to IssueTracker.

fricas

integrate(sin(x^2),x)

$\label{eq1}{fresnelS \left({x \ {\sqrt{2 \over \pi}}}\right)}\over{\sqrt{2 \over \pi}}$

(1)

Type: Union(Expression(Integer),...)

However:

 Maxima 5.9.1.1cvs http://maxima.sourceforge.net
 Using Lisp CMU Common Lisp 19b (19B)
 Distributed under the GNU Public License. See the file COPYING.
 Dedicated to the memory of William Schelter.
 This is a development version of Maxima. The function bug_report()
 provides bug reporting information.
 (%i1) integrate(sin(x**2),x);
                                             (sqrt(2) %i + sqrt(2)) x
 (%o1) sqrt(%pi) ((sqrt(2) %i + sqrt(2)) erf(------------------------)
                                                        2
                                     (sqrt(2) %i - sqrt(2)) x
        + (sqrt(2) %i - sqrt(2)) erf(------------------------))/8
                                                2
 (%i2)

(new) --Bob McElrath?, Fri, 19 Aug 2005 11:13:26 -0500 reply

It's worth noting that both Maple and Mathematica produce the FresnelS? function, which is defined in terms of this integral. However, the Maxima answer appears to be correct as well.

... --kratt6, Sat, 20 Aug 2005 09:40:27 -0500 reply

Note that

is a D-finite function:

fricas

f := sin(x^2)

$\label{eq2}\sin \left({{x}^{2}}\right)$

(2)

Type: Expression(Integer)

fricas

4*x^3*f - D(f,x) + x*D(f,x,2)

$\label{eq3}0$

(3)

Type: Expression(Integer)

Thus integration should be "easy"...

... --kratt6, Sat, 20 Aug 2005 11:05:33 -0500 reply

I browsed the web a little more and came to the conclusion that the Risch algorithm only deals with elementary functions whose integral is elementary, too. "Clearly" (looking at maxima's output or browsing the web),

does not have an elementary antiderivative. Hence, I suspect that it should be treated just as $e^{-x^2}$ , by the PatternMatchIntegration package. This would not be too difficult, probably.

Martin

... --kratt6, Sat, 20 Aug 2005 12:30:16 -0500 reply

Although this does not really solve the original problem, I think I found a bug in INTPM. Currently, there is an operation:

           pmComplexintegrate(f, x) ==
             (rc := splitConstant(f, x)).const ^= 1 =>
               (u := pmintegrate(rc.nconst, x)) case "failed" => "failed"
               rec := u::ANS
               [rc.const * rec.special, rc.const * rec.integrand]
             cse := (rec := matcherfei(f, x, true)).which
             cse = ERF  => [rec.coeff * erf(rec.exponent), 0]
             "failed"

It is pretty obvious that the third line should read:

               (u := pmComplexintegrate(rc.nconst, x)) case "failed" => "failed"

instead. If we perform this change, we get instead of

fricas

complexIntegrate(-%i/2*e^(%i*x^2),x)

$\label{eq4}{{\sqrt{\pi}}\ {\sqrt{-{i \ {\log \left({e}\right)}}}}\ {\erf \left({x \ {\sqrt{-{i \ {\log \left({e}\right)}}}}}\right)}}\over{4 \ {\log \left({e}\right)}}$

(4)

Type: Expression(Complex(Integer))

the correct answer, same for

fricas

complexIntegrate(%i/2*e^(-%i*x^2),x)

$\label{eq5}{{\sqrt{\pi}}\ {\sqrt{i \ {\log \left({e}\right)}}}\ {\erf \left({x \ {\sqrt{i \ {\log \left({e}\right)}}}}\right)}}\over{4 \ {\log \left({e}\right)}}$

(5)

Type: Expression(Complex(Integer))

For some reason, it still won't do

fricas

complexIntegrate(sin(x^2),x)

$\label{eq6}\int^{ \displaystyle x}{{{-{{\sqrt{- 1}}\ {{{e}^{{{\%A}^{2}}\ {\sqrt{- 1}}}}^{2}}}+{\sqrt{- 1}}}\over{2 \ {{e}^{{{\%A}^{2}}\ {\sqrt{- 1}}}}}}\ {d \%A}}$

(6)

Type: Expression(Integer)

Curiously, the pattern matcher is not even invoked in this case... Even if we enter the integral as

fricas

complexIntegrate(-%i/2*e^(%i*x^2)+%i/2*e^(-%i*x^2),x)

$\label{eq7}\int^{ \displaystyle x}{{{-{i \ {{{e}^{i \ {{\%L}^{2}}\ {\log \left({e}\right)}}}^{2}}}+ i}\over{2 \ {{e}^{i \ {{\%L}^{2}}\ {\log \left({e}\right)}}}}}\ {d \%L}}$

(7)

Type: Expression(Complex(Integer))

it fails, although in this case the pattern matcher is invoked. It would need to be invoked on each summand seperately, though.

Martin

... --kratt6, Sat, 20 Aug 2005 15:56:51 -0500 reply

There is another issue I don't quite understand. Currently axiom returns the whole integral unevaluated if it does not manage to evaluate it completely. If it were not for the bug #199, the following were an example:

fricas

integrate(exp(x)/x+exp(-x^2)+1/log(x),x)

$\label{eq8}{{{\erf \left({x}\right)}\ {\sqrt{\pi}}}+{2 \ {li \left({x}\right)}}+{2 \ {Ei \left({x}\right)}}}\over 2$

(8)

Type: Union(Expression(Integer),...)

Although axiom produces the intermediate result integrate(exp(x)/x+exp(-x^2),x)+li(x), it discards it. The responsible line of code is in rinteg$FSINT :

    rinteg(i, f, x, h, comp) ==
      not elem? i => integral(f, x)$F
      empty? rest(l := [mkPrimh(f, x, h, comp) for f in expand i]) => first l
      l

Does this make sense?

Martin

... --kratt6, Mon, 24 Dec 2007 02:42:09 -0800 reply

partially fixed in FriCAS

... --test1, Wed, 23 Apr 2014 18:45:42 +0000 reply

Status: open => closed

Subject: (replying) Be Bold !!

	( 14 subscribers )