Challenging Integral: Evaluating \int 3e^{2x}\sqrt{x+1}\ dx for ODEs

[twoflower]

Hi,

could you please help me in evaluating this integral?

<br /> \int 3e^{2x}\sqrt{x+1}\ dx<br />

I got to this point during solving one ODE and have no idea what to do with this one.

Thank you.

 

[benorin]

not pretty, Maple says it involves an erfc.

 

[twoflower]

not pretty, Maple says it involves an erfc.

Bad, very bad...

 

[fleon]

Hi twoflower

I solve your integral, and found the following answer

$$\frac{3(4e^{2+2x}\sqrt{x+1}-\sqrt{2\pi}Erfi(\sqrt{2(1+x)}))}{8e^{2}}$$

when

Erfi[z] gives the imaginary error function erf(iz)/i

 

[twoflower]

Hi twoflower
I solve your integral, and found the following answer
$$\frac{3(4e^{2+2x}\sqrt{x+1}-\sqrt{2\pi}Erfi(\sqrt{2(1+x)}))}{8e^{2}}$$
when
Erfi[z] gives the imaginary error function erf(iz)/i

Thank you, but that's definitely not what I was supposed to get..Kind of ugly for ODE solution, isn't it? :)

 

[fargoth]

try y=\sqrt{x+1}
and then pick 4ye^{2y^2} as dv and y for u... (http://mathworld.wolfram.com/IntegrationbyParts.html"

 

[benorin]

What was the ODE?

By the way, you can use www.integrals.com

 

[fleon]

Sorry, I'm learning to use this editor, the answer is as follow

\frac{3(4e^{2+2x}\sqrt{x+1}-\sqrt{2\pi}Erfi(\sqrt{2(1+x)}))}{8e^{2}}.

 

[fargoth]

oh yeah, and i forgot, after the integration in parts i suggested above.. you know the trick for e^{x^2} right?

 

[twoflower]

Thank you fargoth, going to try that.

What was the ODE?

By the way, you can use www.integrals.com

The ODE is

<br /> y&#039;&#039; + 2y&#039; + y = 3e^{x}\sqrt{x+1}<br />

 

[benorin]

is it e^{2x}\mbox{ or }e^{x}?

 

[twoflower]

is it e^{2x}\mbox{ or }e^{x}?

It's e^{x}

 

[benorin]

Well, if it is either of the above, then your answer is just nasty (erfc, and all that). But if, perhaps, you had, say y^{\prime\prime}+2y^{\prime} + y= 3e^{-x}\sqrt{x+1}, then your answer is rather sweet, namely y(x)=c_{1}e^{-x}+c_{2}xe^{-x}+\frac{4}{5}e^{-x}(x+1)^{\frac{5}{2}}. In fact, if the exponent of e in the nonhomogeneous part is anything other than -1, erfc (computer aided educated guess).

 

[twoflower]

Well, if it is either of the above, then your answer is just nasty (erfc, and all that). But if, perhaps, you had, say y^{\prime\prime}+2y^{\prime} + y= 3e^{-x}\sqrt{x+1}, then your answer is rather sweet, namely y(x)=c_{1}e^{-x}+c_{2}xe^{-x}+\frac{4}{5}e^{-x}(x+1)^{\frac{5}{2}}. In fact, if the exponent of e in the nonhomogeneous part is anything other than -1, erfc (computer aided educated guess).

You're an oracle, benorin! Professor wrote in wrongly, there has to be e^{-x} in the original ODE.

 

[saltydog]

Why should the mildly messy format of the first solution discourage you? The ODE as first written is perfectly solvable even for the second solution (assuming you're using variation of parameter and B(x) is defined as above) with the second one:

A(x)=-3\int xe^{2x}\sqrt{x+1}

That's just me though.

Edit: Oh yea, I haven't (yet) checked this.

Edit2: Oh yea: equal rights for special functions. You know I have to send a dollar to Lurflurf each time I use that saying don't you?

 

[twoflower]

Well, if it is either of the above, then your answer is just nasty (erfc, and all that). But if, perhaps, you had, say y^{\prime\prime}+2y^{\prime} + y= 3e^{-x}\sqrt{x+1}, then your answer is rather sweet, namely y(x)=c_{1}e^{-x}+c_{2}xe^{-x}+\frac{4}{5}e^{-x}(x+1)^{\frac{5}{2}}. In fact, if the exponent of e in the nonhomogeneous part is anything other than -1, erfc (computer aided educated guess).

Now I computed it and got another solution..

And the official one our proffesor has on his webpage is

<br /> y = \frac{3}{2} \log \left(x^2 + \sqrt{x^4 + 1}\right)e^{-x} -2(x+1)^{\frac{3}{2}}xe^{-x} + \alpha e^{-x} + \beta xe^{-x}<br />

 

[twoflower]

Why should the mildly messy format of the first solution discourage you? The ODE as first written is perfectly solvable even for the second solution (assuming you're using variation of parameter and B(x) is defined as above) with the second one:

A(x)=-3\int xe^{2x}\sqrt{x+1}

That's just me though.

Edit: Oh yea, I haven't (yet) checked this.

Edit2: Oh yea: equal rights for special functions. You know I have to send a dollar to Lurflurf each time I use that saying don't you?

What??

 

[saltydog]

Hey Twoflower, what's going on with this? Your initial ODE was:

y^{&#039;&#039;}+2y^{&#039;}+y=3e^x\sqrt{1+x}

Using variation of parameter, I get:

y(x)=A(x)e^{-x}+B(x)xe^{-x}+c_1e^{-x}+c_2xe^{-x}

with:

B(x)=\frac{3\left(4e^{2+2x}\sqrt{1+x}-\sqrt{2\pi}\text{Erfi}(\sqrt{2(1+x)})\right)}{8e^2}

A(x)=-\frac{3\left(4e^{2+2x}\sqrt{1+x}(4x-3)+7\sqrt{2\pi}\text{Erfi}(\sqrt{2(1+x)})\right)}{32e^2}

Back-substituting y(x) (via Mathematica) into the LHS of the ODE, I get the RHS.

When I back-substitute the "official solution" you posted above into the LHS of the ODE, I get an expression not even close to the RHS, minus x or not.

 

[saltydog]

Really, just what is anyway:

\frac{d}{dx}\text{Erfi}[x]

and:

\frac{d^2}{dx^2}\text{Erfi}[x]

And also, I think we should ask Twoflower to solve this for us:

y^{&#039;&#039;}+2y^{&#039;}+y=3e^x\sqrt{1+x},\quad y(0)=0,\;y&#039;(0)=1

not numerically neither: figure out what c1 and c2 are.

 

[twoflower]

Really, just what is anyway:

\frac{d}{dx}\text{Erfi}[x]

and:

\frac{d^2}{dx^2}\text{Erfi}[x]

And also, I think we should ask Twoflower to solve this for us:

y^{&#039;&#039;}+2y^{&#039;}+y=3e^x\sqrt{1+x},\quad y(0)=0,\;y&#039;(0)=1

not numerically neither: figure out what c1 and c2 are.

Well, this will definitely not appear in the test and I would probably not manage it.

 

[saltydog]

Well, this will definitely not appear in the test and I would probably not manage it.

As always, just a suggestion. Do as you wish. And at the bottom of your test, if you're sure you're going to make an A, write "PF powered!".

 

[twoflower]

As always, just a suggestion. Do as you wish. And at the bottom of your test, if you're sure you're going to make an A, write "PF powered!".

Well, a large portion of my analysis success at the faculty so far is PF powered :) Thank you.

 

[benorin]

I'm no oracle, and I didn't use Oracle [the database]: just Maple.

You're an oracle, benorin! Professor wrote in wrongly, there has to be e^{-x} in the original ODE.