Jncik said:
Homework Statement
solve uxx = utt
for 0<x<π when u(0,t) = u(π,t) = 0 and u(x,0) = sin(2x) and ut(π/4,0) = 0
The Attempt at a Solution
i used the method of separation of variables and have found out that
[URL]http://latex.codecogs.com/gif.latex?u%28x,t%29%20=%20%28c_{1}e^{x\sqrt{\lambda}}%20+%20c_{2}e^{-x\sqrt{\lambda}}%29%28C_{1}e^{t\sqrt{\lambda}}%20+%20C_{2}e^{-t\sqrt{\lambda}}%29
[/URL]
This is wrong. I presume that, separating variables so that u(x,t)= X(x)T(t), you found that [itex]X''/X= \lambda= T''/T[/itex] so that you got the two ordinary differential equations [itex]X''= \lambda X[/itex] and [itex]T''= \lambda T[/itex]. What you have here is the solution
if [itex]\lambda> 0[/itex] and you have no reason to believe that.
But you can also say that [itex]u(0,t)= X(0)T(t)= 0[/itex] means either X(0)= 0 or T(t)= 0. Since the latter would make the solution 0 for all t, it cannot satisfy the initial conditions so we must have X(0)= 0. Similarly, [itex]u(\pi, t)= X(\pi)T(t)= 0[/itex] gives [itex]X(\pi)= 0[/itex].
So we are solving [itex]X''= \lambda X[/itex] with boundary conditions [itex]X(0)= X(\pi)= 0[/itex].
Consider three cases: [itex]\lambda= 0[/itex], [itex]\lambda> 0[/itex], [itex]\lambda< 0[/itex]
If [itex]\lambda= 0[/itex], the equation is just X''= 0 and, integrating twice, [itex]X(x)= C_1x+ C_2[/itex]. Then X(0)= C_2= 0 so [itex]X(\pi)= C_1\pi= 0[/itex] which implies that X is identically 0. But then we cannot satisfy the intial condition that u(x, 0)= sin(2x) so that is impossible,
If [itex]\lambda> 0[/itex], we can let [itex]\lambda= \alpha^2[/itex] for [itex]\alpha[/itex] non-zero. The differential equation is [itex]X''= \alpha^2X[/itex] and the general solution is [itex]X(x)= C_1e^{\alpha x}+ C_2e^{-\alpha x}[/itex] (this is what you had). Now we must have [itex]X(0)= C_1+ C_2= 0[/itex] and [itex]X(\pi)= C_1e^{\alpha\pi}+ C_2e^{-\alpha\pi}[/itex]. From the first equation, [itex]C_2= -C_1[/itex] and, putting that into the second, [itex]C_1e^{\alpha\pi}- C_1e^{-\alpha\pi}= C_1(e^{\alpha\pi}- e^{-\alpha\pi})= 0[/itex].
Now, since [itex]\alpha \ne 0[/itex], [itex]e^{\alpha\pi}[/itex] and [itex]e^{-\alpha\p}[/itex] are not equal (one is larger than 1 and the other less than 1) and their difference cannot be 0, so we must have [itex]C_1= 0[/itex] and then [itex]C_2= 0[/itex]. That is, X(x) is identically 0 so it cannot give a u that satisifies the initial condition.
Finally, if [itex]\lambda< 0[/itex], we can let [itex]\lambda= -\alpha^2[/itex] for [itex]\alpha[/itex] non-zer. The differential equation is [itex]X''= -\alpha^2X[/itex] and the general solution to that is [itex]X(x)= C_1cos(\alpha x)+ C_2 sin(\alpha x)[/itex]. Now we must have [itex]X(0)= C_1= 0[/itex] and [itex]X(\pi)= C_2 sin(\alpha\pi)= 0[/itex]. This would give [itex]C_2= 0[/itex] and so again a "trivial" solution
unless [itex]sin(\alpha\pi)= 0[/itex] which
does happen if [itex]\alpha\pi[/itex] is a multiple of [itex]\pi[/itex]- that is, if [itex]\alpha= n[/itex] and [itex]\lambda= -n^2[/itex] with n and integer.
That is, in order to have a non-trivial solution, we must have [itex]\lambda= -n^2[/itex] for n some integer, in that case, [itex]X(x)= Csin(nx)[/itex] for unknown constant C and unknown integer n.
Now, go back to your equation for T: T''= \lambda T= -n
2T. What is the general solution to that?
( Your condition "[itex]u_t(\pi/4, 0)= 0" looks very peculiar- usually it is something like [itex]u(x, 0)= 0[/itex] for all x, but I think [itex]u_t(\pi/4, 0)= 0[/itex] will give the same thing. [itex]u_t(\pi/4, 0)= X(\pi/4)T'(0)= 0[/itex] and [itex]X(\pi/4)\ne 0[/itex] so you must have [itex]T'(\pi/4)= 0[/itex].<br />
<br />
<br />
<blockquote data-attributes="" data-quote="" data-source=""
class="bbCodeBlock bbCodeBlock--expandable bbCodeBlock--quote js-expandWatch">
<div class="bbCodeBlock-content">
<div class="bbCodeBlock-expandContent js-expandContent ">
the problem is that i don't know how to use the given equations to find the final result<br />
<br />
for example<br />
<br />
[URL]http://latex.codecogs.com/gif.latex?u%280,t%29%20=%200%20=%3E%20%28c_{1}%20+%20c_{2}%29%28C_{1}e^{t\sqrt{\lambda}}%20+%20C_{2}e^{-t\sqrt{\lambda}}%29%20=%200%20=%3E%20?[/URL]<br />
<br />
my professor says that we assume that there is a t for which <br />
<br />
[URL]http://latex.codecogs.com/gif.latex?%28C_{1}e^{t\sqrt{\lambda}}%20+%20C_{2}e^{-t\sqrt{\lambda}}%29%20\neq%200[/URL]<br />
<br />
and he finds that c<sub>1</sub> = -c<sub>2</sub><br />
<br />
for<br />
<br />
u(π,t) he says<br />
<br />
[URL]http://latex.codecogs.com/gif.latex?u%28%CF%80,t%29%20=%20%28c_{1}e^{\pi%20\sqrt{\lambda}}%20-%20c_{1}e^{-\pi\sqrt{\lambda}}%29%28C_{1}e^{t\sqrt{\lambda}}%20+%20C_{2}e^{-t\sqrt{\lambda}}%29%20=%200%20=%3E%20%28c_{1}e^{\pi%20\sqrt{\lambda}}%20-%20c_{1}e^{-\pi\sqrt{\lambda}}%29%20=%200%20=%3E%20\lambda%20=%20-n^{2},%20n%20=%200,1,2,3...[/URL]<br />
<br />
so we have<br />
<br />
[URL]http://latex.codecogs.com/gif.latex?%28e^{\pi%20\sqrt{\lambda}}%20-%20e^{-\pi\sqrt{\lambda}}%29%20=%200%20=%3E%20e^{in\pi}%20-%20e^{-in\pi]}%20=%200%20=%3E%202sin%28n\pi%29%20=%200[/URL]<br />
<br />
and from this he says the following which i can't understand<br />
<br />
from u(x,0) = sin(2x) we get<br />
<br />
λ = -4 and c<sub>2</sub> = 1/2<br />
<br />
how did we get this?<br />
<br />
Thanks in advance
</div>
</div>
</blockquote>[/itex]