Question: First order differential system of equation Combustion

I would like to solve the a sytem of 5 differential equations for a combution process and i dont know how to do.

 

my enties are as follows

 

with(plots);
PDEtools[declare](prime = t);
print(`output redirected...`); # input placeholder
derivatives with respect to t of functions of one variable will

   now be displayed with '
MW = 29;
print(`output redirected...`); # input placeholder
                            MW = 29
cp = 1200;
print(`output redirected...`); # input placeholder
                           cp = 1200
hf = 4*10^7;
print(`output redirected...`); # input placeholder
                         hf = 40000000
AF = 16;
print(`output redirected...`); # input placeholder
                            AF = 16
eq1 := (D(F))(t) = 6.19*10^9*exp(-15098/T)*F(t)^.1*(.2*O(t))^1.65;
print(`output redirected...`); # input placeholder
                            8    /  15098\     0.1     1.65
    D(F)(t) = 4.349007352 10  exp|- -----| F(t)    O(t)    
                                 \    T  /                 
eq2 := (D(O))(t) = 16*(D(F))(t);
print(`output redirected...`); # input placeholder
                      D(O)(t) = 16 D(F)(t)
eq3 := (D(Pr))(t) = -17*(D(F))(t);
print(`output redirected...`); # input placeholder
                     D(Pr)(t) = -17 D(F)(t)
eq4 := (D(T))(t) = (D(F))(t)*hf*(8314*T(t))/((cp-8314)*P(t));
print(`output redirected...`); # input placeholder
                           8314 D(F)(t) hf T(t)
                 D(T)(t) = --------------------
                             (cp - 8314) P(t)  
eq5 := (D(P))(t) = P(t)*(D(T))(t)/T(t)-P(0)*(D(T))(t)/T(0);
print(`output redirected...`); # input placeholder
                       P(t) D(T)(t)   P(0) D(T)(t)
             D(P)(t) = ------------ - ------------
                           T(t)           T(0)    
sol := dsolve(subs([eq1, eq2, eq3, eq4, eq5, O(0) = 152.32*10^(-3), P(0) = 101325, T(0) = 753, F(0) = 9.52*10^(-3), Pr(0) = 0]), numeric, [F(t), Pr(t), O(t), T(t), P(t)]);
%;
Error, (in dsolve/numeric/process_input) dependent variables cannot be assigned a procedure. The following do not meet this requirement: [F(t), Pr(t), O(t), T(t), P(t)]

 

Help

Please Wait...