       Re: DSolve for a real function

• To: mathgroup at smc.vnet.net
• Subject: [mg127939] Re: DSolve for a real function
• From: Andreas Talmon l'Armée at smc.vnet.net
• Date: Tue, 4 Sep 2012 05:45:06 -0400 (EDT)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com
• Delivered-to: l-mathgroup@wolfram.com
• Delivered-to: mathgroup-newout@smc.vnet.net
• Delivered-to: mathgroup-newsend@smc.vnet.net
• References: <k1k8np\$90f\$1@smc.vnet.net> <20120901062746.B3D7F687B@smc.vnet.net>

```Hi,

My initial conditions are the following and I am pretty sure that my
solution consists only of a real part.
The Solution has four eigenvalues and they are complex conjugated. With
all variables and all parameters real numbers I must be able to retrieve
a real solution. But being a mathematica newbie,  I do not understand
how to do it with mathematca.

y''[-c] == ic0, y''[c] == ic0, y'''[-c] == ic1, y'''[c] == -ic1

http://dl.dropbox.com/u/4920002/DGL_4th_Order.nb

Clear[sol]

\$Assumptions = {a \[Element] Reals, ic0 \[Element] Reals,
ic1 \[Element] Reals, c \[Element] Reals};

sol[a_, x_] =
y[x] /. DSolve[{y''''[x] + a y[x] == 0, y''[-c] == ic0,
y''[c] == ic0, y'''[-c] == ic1, y'''[c] == -ic1}, y[x],
x][] // FullSimplify

Reduce[Element[sol[a, x], Reals], a, Reals]

Andreas

On 09/01/2012 08:27 AM, Bob Hanlon wrote:
> What are your initial conditions?
>
> Clear[sol]
>
> sol[a_, x_] = y[x] /. DSolve[
>       {y''''[x] + a y[x] == 0,
>        y == ic0, y' == ic1,
>        y'' == ic2, y''' == ic3},
>       y[x], x][] // FullSimplify
>
> (1/(2*a^(3/4)))*
>     (Cosh[(a^(1/4)*x)/Sqrt]*
>          (2*a^(3/4)*ic0*Cos[(a^(1/4)*x)/
>                   Sqrt] + Sqrt*
>               (Sqrt[a]*ic1 + ic3)*
>               Sin[(a^(1/4)*x)/Sqrt]) +
>        (Sqrt*(Sqrt[a]*ic1 - ic3)*
>               Cos[(a^(1/4)*x)/Sqrt] +
>             2*a^(1/4)*ic2*Sin[(a^(1/4)*x)/
>                   Sqrt])*
>          Sinh[(a^(1/4)*x)/Sqrt])
>
> Reduce[Element[sol[a, x], Reals], a, Reals]
>
> a > 0
>
>
> Bob Hanlon
>
>
> On Fri, Aug 31, 2012 at 3:59 AM,  <"Andreas Talmon l'Arm=E9e"@smc.vnet.net> wrote:
>> Hi All
>>
>> Is there a way to tell mathematica to solve only for real solutions. My
>> differential equation is of the kind
>>
>> y''''[x]+a y[x]==0
>>
>> a= constant coefficient
>>
>> I know that I get 4 komplex eigenvalues which are complex conjungated.
>> But y[x] is a real function.
>> Solving this equation with DSolve always gets a complex function y[x].
>>
>> Any Ideas.
>>
>> Thanks, Andreas
>>

```

• Prev by Date: Re: Creating animated 3d graphic in Mathematica for export
• Next by Date: Re: split the sublists into parts according to some rules
• Previous by thread: Re: DSolve for a real function
• Next by thread: Re: DSolve for a real function