It is very likely that they used wolfram alpha, Mathematica, or MATLAB to solve for the two ODE's: F(z) and G(t). You can let u(0,t) = h(t). Then once u(z,t) is solved for, its not hard to find something like that. But again, I would need to see what their PDE was, and the boundary conditions. I...