ut =50uxx
u(0,t) = u(π,t) = 0, t > 0
u(x,0) = { x , 0 < x < π / 2
{ 4 , π / 2 < x < π
Solution
Let u(x,t) = X(x)T(t)
Substitution into equation gives X(x)T’(t) = 50 X’’(x)T(t) => T’(t)/[50 T(t)] = X’’(x)/ X(x)
Variables t and x are independent. So T’(t)/[50 T(t)] = X’’(x)/ X(x) = -k2
NOTE: variant k2 at the right side is impossible because equation X’’(x) - k2 X(x) = 0 with the given boundary conditions have only trivial solution.
Therefore we’ll get two ordinary differential equations
T’(t) + k2 50 T(t)] = 0
X’’(x) + k2 X(x) = 0
From this equations T(t) = Aexp(-k250t), X(x) = Bsin(kx) + Ccos(kx) (A, B, C – arbitrary constants)
From boundary conditions u(0,t) = u(π,t) = 0, t > 0 => X(0) = X(π) = 0 and C = 0 X(x) = Bsin(k π) = 0 => k = ±n (n = 1,2, . . .).
Solution is
From initial condition
Multiplying by sin(mx) and integrating with bounds 0, π
Comments