I'm trying to compute the numerical solution for a PDE using the numerical solution of a first-order ODE. I looked at the INTEGRATE keyword of the pdsolve function, but I'm not sure how to use it. Can someone show an example of how to use this keyword?