I was hoping that by solving the pde the normal way which is part of the question that it might give a clue as to how the change of variables method should should turn out.
Solving the original PDE I get
. Taking the positive solution of y
The IC's are u(x,0)=0, therefore
. Let therefore
I am not sure of this is right either...