Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2001
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2001

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: differential equation with buondary conditions

  • To: mathgroup at smc.vnet.net
  • Subject: [mg30323] Re: differential equation with buondary conditions
  • From: "Allan Hayes" <hay at haystack.demon.co.uk>
  • Date: Sat, 11 Aug 2001 03:39:46 -0400 (EDT)
  • References: <9kqk0f$4g6$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Gustavo,

We have,

DSolve[{y''[x]+k y[x]==0, y[0]==0,y[a]==0},y,x]

{{y -> Function[{x}, 0]}}

This is the generic solution, the only other possibility is for special
values of k and/or a (see below)

Leave out the boundary conditions

sol=DSolve[{y''[x]+k y[x]==0},y,x]

{{y -> Function[{x}, C[1]*Cos[Sqrt[k]*x] +
      C[2]*Sin[Sqrt[k]*x]]}}

Get the function

ys = y/.sol[[1]]

Function[{x}, C[1]*Cos[Sqrt[k]*x] + C[2]*Sin[Sqrt[k]*x]]

The first boundary condition gives

ys[0]==0

        C[1] == 0

Using this on the second boundary condition we get

ys[a]==0/.C[1]->0

C[2]*Sin[a*Sqrt[k]] == 0

This entails C[2] =0 unless Sin[a*Sqrt[k]] =0.
--
Allan
---------------------
Allan Hayes
Mathematica Training and Consulting
Leicester UK
www.haystack.demon.co.uk
hay at haystack.demon.co.uk
Voice: +44 (0)116 271 4198
Fax: +44 (0)870 164 0565

"Gustavo Seabra" <gseabra at swbell.net> wrote in message
news:9kqk0f$4g6$1 at smc.vnet.net...
> Hello,
>
>     I'm trying to make Mathematica solve the following:
>
>                             y''[x] + k y[x] == 0
>
> subject to the boundary conditions:
>         y[x<0] = 0
>         y[x>a] = 0
> so that y[x] != 0 only if 0 < x < a.
> (yes, it's the "particle in a 1-d box problem.)
>
> If I just do: DSolve[{y''[x] + k y[x]  ==  0}, y[x], x]
> it works fine, giving:
> {{y[x] -> C[2] Cos[Sqrt[k] x] + C[1] Sin[Sqrt[k] x]}},
> which is perfectly ok.
>
> But if I include the boundary conditions y[0] == y[a] == 0,
> it doesn't work.
>
>     Any ideas?
> --
> -----------------------------------------------------------------
> Gustavo Seabra - Graduate Student
> Chemistry Department
> Kansas State University
> -----------------------------------------------------------------
>
>




  • Prev by Date: Re: Problem with FourierTransform
  • Next by Date: Re: Trick for getting the comples conjugate in symbolic calculations?
  • Previous by thread: Re: differential equation with buondary conditions
  • Next by thread: Compile[]