[DSolve] Re: newbie question DSolve (revisited again)
- To: mathgroup at smc.vnet.net
- Subject: [mg52226] [DSolve] Re: newbie question DSolve (revisited again)
- From: Paul Abbott <paul at physics.uwa.edu.au>
- Date: Mon, 15 Nov 2004 03:17:35 -0500 (EST)
- Organization: The University of Western Australia
- References: <cn1p31$evf$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
In article <cn1p31$evf$1 at smc.vnet.net>, pdesai1 at umbc.edu wrote:
> The Dirac Delta forces a discontinuity in the system, so we actually have
> two parts of the string
> 1. between 0 and d y1[x]
> 2. between d and l y2[x]
>
> so the b.c are
>
> y1[0]=0
> y2[l]=0
>
> at the constraint
>
> y1[d]=y2[d] Shape is continuous
So, actually,
y1[d] == y2[d] == y[d]
> y1'[d]-y2'[d]=c*lamda*y[d] the jump discontinuity due to the delta function
>
> I am thinking of an alternate construction without Dirac Delta function
>
> y1''[x]-lamda^2*y1[x]=0 0<x<d
>
> with B.C
>
> y1[0]=0
> y1'[d]=y2'[d]-c*lamda*y[d]
>
> and
>
> y2''[x]-lamda^2*y2[x]=0 0<x<d
you mean d < x < l
> y2'[d]=y1'[d]+c*lamda*y[d]
>
> Please let me know if this is applicable and how I can implement in
> Mathematica. Also can mathematica give you Eigenfunctions-- non-trivial
> solutions to a homogenous BVP
Of course it can -- if you state the problem fully. Clearly, however,
the eigenvalue problem requires, in general, the solution to a
transcendental equation.
Here is the first equation with boundary condition and continuity.
sol1 = First[DSolve[{y1''[x] - lambda^2 y1[x] == 0,
y1[0] == 0, y1[d] == y[d]}, y1, x]]
Here is the second equation with boundary condition and continuity.
sol2 = First[DSolve[{y2''[x] - lambda^2 y2[x] == 0,
y2[l] == 0, y2[d] == y[d]}, y2, x]]
Now apply the derivative jump discontinuity condition.
FullSimplify[{
y1'[d] == y2'[d] - c lambda y[d],
y2'[d] == y1'[d] + c lambda y[d]} /. sol2 /. sol1]
As stated, however, c is not determined by the differential equations or
boundary conditions.
Cheers,
Paul
--
Paul Abbott Phone: +61 8 6488 2734
School of Physics, M013 Fax: +61 8 6488 1014
The University of Western Australia (CRICOS Provider No 00126G)
35 Stirling Highway
Crawley WA 6009 mailto:paul at physics.uwa.edu.au
AUSTRALIA http://physics.uwa.edu.au/~paul