Hi, I'm very new to maple and I'm trying to solve a system of ODEs but as of now it has taken over 5 hours to solve and is still evaluating. i feel as though I may have made some basic mistakes in the code which make the calculation much longer than it should be. I will try tyo explain the problem as well as I can
I'm trying to model a chemical reaction, and solve for the concentrations of two species involved in the reaction after a given time.
The rate of the reaction is given as: r=2900exp(-53300/RT)*Cno^0.62*Cnh3^-0.05; where R is a constant, T is temperature and Cno and Cnh3 are the concentrations of NO and NH3 respectively. I am interested in solving for Cno and Cnh3 and getting an expression for each of them. I have tried to set up the system of ODEs as follows:
with(LinearAlgebra):
with(DEtools):
r := 2900*(exp(1)^(-53300/R*T)*CNO(t)^0.62*CNH3(t)^-0.05:
ode := diff(CNO(t), t) = -1* r: (negative because they are decaying with time)
ode2 := diff(CNH3(t), t) = -1* r:
ics := CNO(1020) = 1.6, CNH3(1020) = 1.6; (sets up known initial conditions)
ode := subs(R = 8.314, T = 473, ode):
ode2 := subs(R = 8.314, T = 473, ode2):
sys_ode := (ode, ode2) :
dsolve([sys_ode, ics]);
I wonder if the problem has to do with the boundary conditions that I've set or not. Please help as I know I may have set it up very inefficiently which might be causing the problems. Thank you for your help