How can I solve this PDE kinetics equation with Symbolic Toolbox?

Hello , I want to solve a PDE from a kinetics equation as depicted next:

which

I made the following code ,but dont know what i have to change to make it run.

syms y(t,z)
U=10;
E_l=1e-8;
cyanide=159;
oxygen=8;
r=0.706e-3*(cyanide^1.38)*(oxygen^0.161)*(y-0.17)^2.19

r(t, z) = 

ode = -(diff(y,t))-U*(diff(y,z,1))+E_l*(diff(y,z,2))-r == 0

ode(t, z) = 

ySol = dsolve(ode) 

Error using dsolve>mupadDsolve
Symbolic ODEs must have exactly one independent variable.

Error in dsolve (line 203)
sol = mupadDsolve(args, options);

Thanks in advance for the help.

Answer 1

[Torsten]

PDEs cannot be solved with the symbolic toolbox.

Use "pdepe" as numerical solver instead.