I want to solve the chemistry reaction rate. Here I got the 4 differential equatiobs.
dE/dt=-k1ES+k2ES+kiES dES/dt=k1ES-kiES-k2ES dS/dt=-k1ES+kiES dP/dt=k2ES
Here is my code?
k1,k2,ki,t=var('k1 k2 ki t') E=function('E',t) ES=function('ES',t) S=function('S',t) P=function('P',t)
des=(diff(E,t)==-k1ES+k2ES+kiES,diff(ES,t)==k1ES-kiES-k2ES diff(S,t)==-k1ES+kiES diff(P,t)==k2ES)
desolve_system(des,[E,ES,S,P],ivar=t,ics=[0,1,1,1,1])
The above is my code. But I just cannot execute it. I am a rookie with this kind of coding style.
Thanks for answering.
