Computer Simulation and Modeling

Author: ForhanShahriarFahim

CPM_input2.csv

@@ -0,0 +1,8 @@
+"id","name","time","pre"
+0,A,8,
+1,B,2,
+2,C,5,0
+3,D,4,1
+4,E,6,2 3
+5,F,7,4
+6,G,7,5

Chemical-Reactor.py

@@ -14,9 +14,9 @@ def chemical_reactor(A, B, C, k1, k2, end_time, dt):
         substance_C.append(C)
         time_points.append(time)
         time += dt
-
-        A += (k2 * C - k1 * A * B) * dt
-        B += (k2 * C - k1 * A * B) * dt
+        val = k2 * C - k1 * A * B
+        A += val * dt
+        B += val * dt
         C += (2 * k1 * A * B - 2 * k2 * C) * dt
 
         plt.clf()