Export meshdata and solution from pdetool

I need to run the pdetool analysis several times for different values of pde coefficients. Repeating this procedure for different meshes is not convenient using the pdetool GUI.

Is there a way to automate the whole process via a .m file? The .m file should read the pde coefficients via the Matlab command prompt and should return the mesh (p, e, t) and solution (u) after the pde solution is complete.

If I understood correctly, the geometry does not change and only the pde coefficients change. If this is the case, I would recommend a one-time exporting of the mesh and the BCs from pdetool and saving it to disk. Then, rather than go through pdetool for defining coefficients, using the command-line function: assempde/parabolic/hyperolic, programmatically enter coefficients and solve. Of course you will also need to load the mesh and BCs from the matfile that you've saved to prior.
This way you need to use pdetool exactly once.