I have the following PDE system steaming from Flash Photolysis:

when I use pdsolve(pdesys,[J,C]) I get:

The solution appears to be either C(x,t) = 0 or J(x,t) = 0. This are the obvious solutions (0 = 0). I have the analytical solution to this PDE system where neither C(x,t) nor J(x,t) are 0.

How to solve this system in maple? Thanks.

How do I create a function from dsolve() result?

For example, dsolve() outputs:

u(t) = u₀e^{a t}

then I would like to create function:

u := t → u₀e^{a t}

but I want to do it using dsolve() output, not typing it manually. Thanks.