Nan's profile - activity

oh,this page cannot be post the symbol like the star which means multiplication...

the four ode should be:

sage: de1=diff(E,t)==-k1*E*S+k2*ES+ki*ES
sage: de2=diff(ES,t)==k1*E*S-ki*ES-k2*ES
sage: de3=diff(S,t)==-k1*E*S+ki*ES
sage: de4=diff(P,t)==k2*ES

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.