/*
  生体システム工学の基礎
  図4.69
  gcc -o mm mm.c -l
 */
#include <stdio.h>
#include <math.h>

int main(int argc, char *argv[])
{
  int i;
  double kf, kr, kcat;
  double s, e, es, p;
  double dt;

  s = 1.0;
  es = 0.0;
  e = 1e-1;
  p = 0.0;
  kf = 1e0;
  kr = 1e-2;
  kcat = 0.5;
  dt = 1e-3;
  
  for(i=0; i<60000; i++) {
    /* 時間，基質，酵素，基質酵素結合物，生成物 */
    printf("%lf %lf %lf %lf %lf\n", i*dt, s, e, es, p);
    s = s + dt * (-kf * e * s + kr * es);
    e = e + dt * (-kf * e * s + kr * es + kcat * es);
    es = es + dt * (kf * e * s - kr * es - kcat * es);
    p = p + dt * (kcat * es);
 }
}