Runge-Kutta i excel (koblede differentiallignigner)

Hej Eksperter...

Jeg har et problem, jeg skal løse et differentialligningssystem i Excel via Runge-Kutta metoden med skridtet h=0.001 og finde E,S,ES,Z til tiden t=0.5

Den kemiske reaktion:
E+S -(c1f)-> ES -(c2)-> Z+E
hvor ES kan "løbe" tilbage ES -(c1b)-> E+S
(parentesen angiver hastigheden)

Følgende system betragtes:
dE/dt = c1b * ES + c2 * ES - c1f * E * S (defineres som funktion RKfunc1)
dS/dt = c1b * ES - c1 * E * S) (defineres som funktion RKfunc2)
dES/dt = c1 * E * S - c1b * ES - c2 * ES (defineres som funktion RKfunc3)
dZ/dt = c2 * ES (defineres som funktion RKfunc4)

Men hvordan defineres konstanterne k1E, K2E, K3E, K4E … K1Z, K2Z, K3Z, K4Z?

Jeg har defineret differentialligningerne i VBA med de overstående navne som funktioner af (E,S,ES,Z)

Håber der er nogle Matematik genier der kan hjælpe
selvom det ser meget uoverskueligt ud i notepad form :)