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