       RE: Oh, those complex values!!!

• To: mathgroup at smc.vnet.net
• Subject: [mg72127] RE: [mg72068] Oh, those complex values!!!
• From: "David Park" <djmp at earthlink.net>
• Date: Wed, 13 Dec 2006 06:39:16 -0500 (EST)

```Aaron,

I would much prefer to write your f definition as follows...

f[a_, b_, c_, d_][x_] = x Log[a + b x + Sqrt[c + 2d x + x^2]];

If you do an indefinite integration you must always remember that the answer
is the returned value plus a constant. That constant could be a complex
number that undoes the imaginary number you get in your result. Also the
returned result contains multivalued functions and this complicates the
picture more.

A general definite integral seems to take too long, but if you integrate
with specific values of a, b, c, d you can obtain exact real results in
reasonable time.

Integrate[f[2, 2, 10, 1][x], {x, 0, 1}]
% // N

NIntegrate[f[2, 2, 10, 1][x], {x, 0, 1}]

David Park

From: aaronfude at gmail.com [mailto:aaronfude at gmail.com]

Integrating a real function, getting complex values again. How to get
the correct answer here? How to systematically avoid this pitfall?

f=x Log[a+b x+Sqrt[c+2d x+x^2]];
t = Assuming[c>d^2&&a>0&&b>0,Integrate[f, x]];
u=t/.a\[Rule]2/.b\[Rule]2/.c\[Rule]10/.d\[Rule]1;

N[u/.x\[Rule]1-u/.x\[Rule]0]
9.50705\[InvisibleSpace]-29.7626 \[ImaginaryI]

NIntegrate[f/.a\[Rule]2/.b\[Rule]2/.c\[Rule]10/.d\[Rule]1, {x, 0, 1}]
0.954442