dear, Jakubi,
thanks for observing that,but in fact it produced no goods.
the corrected program is :
> with(PDEtools);

rho[f0] := 650;
epsilon0 := .3;
epsilonb := .9;
R := 8.314;
k1 := A[1]*exp(-E[1]/(R*T(z, t)));
k2 := A[2]*exp(-E[2]/(R*T(z, t)));
k3 := A[3]*exp(-E[3]/(R*T(z, t)));
A[1] := 1200;
A[2] := 1200;
A[3] := 800;
E[1] := 25000;
E[2] := 22100;
E[3] := 5000;
Cp := 3*T(z, t)+12*10^(-3)*T(z, t)^2; k := 2*10^(-1)*T(z, t);
h := .345*T(z, t); d := 0.3e-2*T(z, t); e := .9; s := 1.38*10^(-23);
ODE := {diff(epsilon(t), t) = k3*(epsilonb-epsilon(t)), diff(rho[TG](t), t) = k1*rho[f](t), diff(rho[c](t), t) = k2*rho[f](t), diff(rho[f](t), t) = -(k1+k2)*rho[f](t)};
ICE := {epsilon(0) = epsilon0, rho[TG](0) = 0, rho[c](0) = 0, rho[f](0) = rho[f0]};
PDE := {diff(Cp*T(z, t), t) = diff(k*(diff(T(z, t), z)), z)+255000*k1*Y(z, t)*rho[f](t), diff(Y(z, t), t) = diff(d*(diff(Y(z, t), z)), z)-k1*Y(z, t)};
IBC := {T(z, 0) = 298, Y(z, 0) = 0, (D[1](T))(0, t) = 0, (D[1](T))(.1, t) = 0.345e-2*T(.1, t)*(1323-T(.1, t))/(20*10^(-1)*T(.1, t))+e*s*(1323^4-T(.1, t)^4), (D[1](Y))(0, t) = 0, (D[1](Y))(.1, t) = 450/(.3*T(.1, t))};
sys := {ODE, PDE};
CON := {IBC, ICE};
pdsolve(sys, CON, numeric);

Error, (in pdsolve/numeric/process_PDEs) number of dependent variables and number of PDE must be the same
 

probaly there is not a typing error but something in the solving  procedure or the pdsolve options
i cant really understand what this message that i get is trying to say.several examples in maple help show similar systems but icant get a result...
help me please!!!!!
 

dear, Jakubi,
thanks for observing that,but in fact it produced no goods.
the corrected program is :
> with(PDEtools);

rho[f0] := 650;
epsilon0 := .3;
epsilonb := .9;
R := 8.314;
k1 := A[1]*exp(-E[1]/(R*T(z, t)));
k2 := A[2]*exp(-E[2]/(R*T(z, t)));
k3 := A[3]*exp(-E[3]/(R*T(z, t)));
A[1] := 1200;
A[2] := 1200;
A[3] := 800;
E[1] := 25000;
E[2] := 22100;
E[3] := 5000;
Cp := 3*T(z, t)+12*10^(-3)*T(z, t)^2; k := 2*10^(-1)*T(z, t);
h := .345*T(z, t); d := 0.3e-2*T(z, t); e := .9; s := 1.38*10^(-23);
ODE := {diff(epsilon(t), t) = k3*(epsilonb-epsilon(t)), diff(rho[TG](t), t) = k1*rho[f](t), diff(rho[c](t), t) = k2*rho[f](t), diff(rho[f](t), t) = -(k1+k2)*rho[f](t)};
ICE := {epsilon(0) = epsilon0, rho[TG](0) = 0, rho[c](0) = 0, rho[f](0) = rho[f0]};
PDE := {diff(Cp*T(z, t), t) = diff(k*(diff(T(z, t), z)), z)+255000*k1*Y(z, t)*rho[f](t), diff(Y(z, t), t) = diff(d*(diff(Y(z, t), z)), z)-k1*Y(z, t)};
IBC := {T(z, 0) = 298, Y(z, 0) = 0, (D[1](T))(0, t) = 0, (D[1](T))(.1, t) = 0.345e-2*T(.1, t)*(1323-T(.1, t))/(20*10^(-1)*T(.1, t))+e*s*(1323^4-T(.1, t)^4), (D[1](Y))(0, t) = 0, (D[1](Y))(.1, t) = 450/(.3*T(.1, t))};
sys := {ODE, PDE};
CON := {IBC, ICE};
pdsolve(sys, CON, numeric);

Error, (in pdsolve/numeric/process_PDEs) number of dependent variables and number of PDE must be the same
 

probaly there is not a typing error but something in the solving  procedure or the pdsolve options
i cant really understand what this message that i get is trying to say.several examples in maple help show similar systems but icant get a result...
help me please!!!!!
 

dear jakubi,

thanks for replying to my question.

the system tha that i have to solve is :
 

> with(PDEtools);

rho[f0] := 650:
epsilon0 := .3:
epsilonb := .9:
R := 8.314:
k1 := A[1]*exp(-E[1]/(R*T(t, z))):
k2 := A[2]*exp(-E[2]/(R*T(t, z))):
k3 := A[3]*exp(-E[3]/(R*T(t, z))):
A[1] := 1200:
A[2] := 1200:
A[3] := 800:
E[1] := 25000:
E[2] := 22100:
E[3] := 5000:
Cp := 3*T(z, t)+12*10^(-3)*T(z, t)^2:
k := 2*10^(-1)*T(z, t):
h := .345*T(z, t):
d := 0.3e-2*T(z, t):
e := .9:
s := 1.38*10^(-23):

ODE := {diff(epsilon(t), t) = k3*(epsilonb-epsilon(t)), diff(rho[TG](t), t) = k1*rho[f](t), diff(rho[c](t), t) = k2*rho[f](t), diff(rho[f](t), t) = -(k1+k2)*rho[f](t)};

ICE := {epsilon(0) = epsilon0, rho[TG](0) = 0, rho[c](0) = 0, rho[f](0) = rho[f0]};

PDE := {diff(Cp*T(z, t), t) = diff(k*(diff(T(z, t), z)), z)+255000*k1*Y(z, t)*rho[f](t), diff(Y(z, t), t) = diff(d*(diff(Y(z, t), z)), z)-k1*Y(z, t)};
 

IBC := {T(z, 0) = 298, Y(z, 0) = 0, (D[1](T))(0, t) = 0, (D[1](T))(.1, t) = 0.345e-2*T(.1, t)*(1323-T(.1, t))/(20*10^(-1)*T(.1, t))+e*s*(1323^4-T(.1, t)^4), (D[1](Y))(0, t) = 0, (D[1](Y))(.1, t) = 450/(.3*T(.1, t))};
sys := {ODE, PDE};
CON := {IBC, ICE};
pdsolve(sys, CON, numeric);

Error, (in pdsolve/numeric/process_PDEs) number of dependent variables and number of PDE must be the same

this is what i get while trying to run the program.
i need desperately your help!!!!!
thank you
 

dear jakubi,

thanks for replying to my question.

the system tha that i have to solve is :
 

> with(PDEtools);

rho[f0] := 650:
epsilon0 := .3:
epsilonb := .9:
R := 8.314:
k1 := A[1]*exp(-E[1]/(R*T(t, z))):
k2 := A[2]*exp(-E[2]/(R*T(t, z))):
k3 := A[3]*exp(-E[3]/(R*T(t, z))):
A[1] := 1200:
A[2] := 1200:
A[3] := 800:
E[1] := 25000:
E[2] := 22100:
E[3] := 5000:
Cp := 3*T(z, t)+12*10^(-3)*T(z, t)^2:
k := 2*10^(-1)*T(z, t):
h := .345*T(z, t):
d := 0.3e-2*T(z, t):
e := .9:
s := 1.38*10^(-23):

ODE := {diff(epsilon(t), t) = k3*(epsilonb-epsilon(t)), diff(rho[TG](t), t) = k1*rho[f](t), diff(rho[c](t), t) = k2*rho[f](t), diff(rho[f](t), t) = -(k1+k2)*rho[f](t)};

ICE := {epsilon(0) = epsilon0, rho[TG](0) = 0, rho[c](0) = 0, rho[f](0) = rho[f0]};

PDE := {diff(Cp*T(z, t), t) = diff(k*(diff(T(z, t), z)), z)+255000*k1*Y(z, t)*rho[f](t), diff(Y(z, t), t) = diff(d*(diff(Y(z, t), z)), z)-k1*Y(z, t)};
 

IBC := {T(z, 0) = 298, Y(z, 0) = 0, (D[1](T))(0, t) = 0, (D[1](T))(.1, t) = 0.345e-2*T(.1, t)*(1323-T(.1, t))/(20*10^(-1)*T(.1, t))+e*s*(1323^4-T(.1, t)^4), (D[1](Y))(0, t) = 0, (D[1](Y))(.1, t) = 450/(.3*T(.1, t))};
sys := {ODE, PDE};
CON := {IBC, ICE};
pdsolve(sys, CON, numeric);

Error, (in pdsolve/numeric/process_PDEs) number of dependent variables and number of PDE must be the same

this is what i get while trying to run the program.
i need desperately your help!!!!!
thank you