An alternate solution procedure is to define a new unknown:
You must derive the problem for v:
The boundary condition is simply:
To get the PDE for , use
Similarly, for the initial conditions:
After finding , I still need to find
from the definition of
: