diff --git a/dependencies/asm3_class.js b/dependencies/asm3_class.js
index 2b85fe5..1b214db 100644
--- a/dependencies/asm3_class.js
+++ b/dependencies/asm3_class.js
@@ -5,7 +5,7 @@ class ASM3 {
 
     // Hydrolysis
     k_H:      3.,   // hydrolysis rate constant       [g X_S g-1 X_H d-1]
-    K_x:      1.,   // hydrolysis saturation constant [g X_S g-1 X_H] 
+    K_X:      1.,   // hydrolysis saturation constant [g X_S g-1 X_H] 
     // Heterotrophs
     k_STO:    5.,   // storage rate constant          [g S_S g-1 X_H d-1]
     nu_NO:    0.6,  // anoxic reduction factor        [-]
@@ -16,16 +16,16 @@ class ASM3 {
     mu_H_max: 2.,   // maximum specific growth rate   [d-1]
     K_NH:     0.01, // saturation constant S_NH3      [g NH3-N m-3]
     K_HCO:    0.1,  // saturation constant S_HCO      [mole HCO3 m-3]
-    b_H_O2:   0.2,  // aerobic respiration rate       [d-1]
+    b_H_O:    0.2,  // aerobic respiration rate       [d-1]
     b_H_NO:   0.1,  // anoxic respiration rate        [d-1]
-    b_STO_O2: 0.2,  // aerobic respitation rate X_STO [d-1]
+    b_STO_O:  0.2,  // aerobic respitation rate X_STO [d-1]
     b_STO_NO: 0.1,  // anoxic respitation rate X_STO  [d-1]
     // Autotrophs
     mu_A_max: 1.0,  // maximum specific growth rate   [d-1]
     K_A_NH:   1.,   // saturation constant S_NH3      [g NH3-N m-3] 
     K_A_O:    0.5,  // saturation constant S_0        [g O2 m-3]
     K_A_HCO:  0.5,  // saturation constant S_HCO      [mole HCO3 m-3]
-    b_A_O2:   0.15, // aerobic respiration rate       [d-1]
+    b_A_O:    0.15, // aerobic respiration rate       [d-1]
     b_A_NO:   0.05  // anoxic respiration rate        [d-1]
   }
 
@@ -51,13 +51,43 @@ class ASM3 {
     i_TS:     0.90, // TSS content X_H / X_A          [g TS g-1 X_H / X_A]
     i_TS:     0.60  // TSS content X_STO (PHB based)  [g TS g-1 X_STO]
   }
+  // S_O, S_I, S_S, S_NH, S_N2, S_NO, S_HCO, X_I, X_S, X_H, X_STO, X_A, X_TS
+  state = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
 
-  constructor() {
+  constructor(state) {
+    this.state = state;
+  }
+
+  _monod(c, K){
+    return c / (K + c);
+  }
+
+  _inv_monod(c, K){
+    return K / (K + c);
   }
 
   compute_rates(state) {
     const rates = new Array(12);
-    rates[0] = this.parameters[]
+    const [S_O, S_I, S_S, S_NH, S_N2, S_NO, S_HCO, X_I, X_S, X_H, X_STO, X_A, X_TS] = state;
+    
+    // Hydrolysis
+    rates[0] = this.kin_params.k_H * this._monod(X_S / X_H, this.kin_params.K_X) * X_H;
+
+    // Heterotrophs
+    rates[1] = this.kin_params.k_STO * this._monod(S_O, this.kin_params.K_O) * this._monod(S_S, this.kin_params.K_S) * X_H;
+    rates[2] = this.kin_params.k_STO * this.kin_params.nu_NO * this._inv_monod(S_O, this.kin_params.K_O) * this._monod(S_NO, this.kin_params.K_NO) * this._monod(S_S, this.kin_params.K_S) * X_H;
+    rates[3] = this.kin_params.mu_H_max * this._monod(S_O, this.kin_params.K_O) * this._monod(S_NH, this.kin_params.K_NH) * this._monod(S_HCO, this.kin_params.K_HCO) * this._monod(X_STO/X_H, K_STO) * X_H;
+    rates[4] = this.kin_params.mu_H_max * this.kin_params.nu_NO * this._inv_monod(S_O, this.kin_params.K_O) * this._monod(S_NO, this.kin_params.NO) * this._monod(S_NH, this.kin_params.K_NH) * this._monod(S_HCO, this.kin_params.K_HCO) * this._monod(X_STO/X_H, this.kin_params.K_STO) * X_H;
+    rates[5] = this.kin_params.b_H_O * this._monod(S_O, this.kin_params.K_O) * X_H;
+    rates[6] = this.kin_params.b_H_NO * this._inv_monod(S_O, this.kin_params.K_O) * this._monod(S_NO, this.kin_params.K_NO) * X_H;
+    rates[7] = this.kin_params.b_STO_O * this._monod(S_O, this.kin_params.K_O) * X_H;
+    rates[8] = this.kin_params.b_STO_NO * this._inv_monod(S_O, this.kin_params.K_O) * this._monod(S_NO, this.kin_params.K_NO) * X_STO;
+    
+    // Autotrophs
+    rates[9] = this.kin_params.mu_A_max * this._monod(S_O, this.kin_params.K_A_O) * this._monod(S_NH, this.kin_params.K_A_NH) * this._monod(S_HCO, this.kin_params.K_A_HCO) * X_A;
+    rates[10] = this.kin_params.b_A_O * this._monod(S_O, this.kin_params.K_O) * X_A;
+    rates[11] = this.kin_params.b_A_NO * this._inv_monod(S_O, this.kin_params.K_A_O) * this._monod(S_NO, this.kin_params.K_NO) * X_A;
+
     return rates;
   }
 }
\ No newline at end of file