function dydt = sub_simulasi (t,y);
k1 = 0.003815; k2 = 0.000123;

cA = y; cH2O = 0;
cME = 0; cMOH = y;
dydt  = -k1*cA*cMOH - (k2*cME*cH2O);
dydt = dydt';