Question: plot solution of PDE

restart;
with(Plot);
params := {alpha = 2.5, k = 3, w = 2, beta[3] = 3, beta[4] = 1.7};
xi := beta[3]*(2*alpha*k*t + x)*sqrt(1/(36*alpha*beta[4] + 36*gamma*beta[4]));
params := {alpha = 2.5, k = 3, w = 2, beta[3] = 3, beta[4] = 1.7}

          xi := beta[3] (2 alpha k t + x) 

                                                 (1/2)
            /                 1                 \     
            |-----------------------------------|     
            \36 alpha beta[4] + 36 gamma beta[4]/     


sol1n := u(x, t) = U(xi)*exp((-sqrt(1/(36*alpha*beta[4] + 36*gamma*beta[4]))*x + w*t + theta)*I);
                     /                          
                     |                          
 sol1n := u(x, t) = U|beta[3] (2 alpha k t + x) 
                     \                          

                                        (1/2)\    /  /
   /                 1                 \     |    |  |
   |-----------------------------------|     | exp|I |
   \36 alpha beta[4] + 36 gamma beta[4]/     /    \  \
                                       (1/2)                \\
  /                 1                 \                     ||
 -|-----------------------------------|      x + w t + theta||
  \36 alpha beta[4] + 36 gamma beta[4]/                     //



plot3d(rhs(sol1n), x = 0 .. 250, t = 0 .. 4);

how plot the the solution of PDE of this kind of function?

Download plot.mw