1. #REPLACE NAMES OF THE VARIABLES USED HERE TO ONES YOU USED IN THE SIMULATION
2. modelC2.resetAll()
3. resultC2=modelC2.simulate(start=0, end=1000, points=10000)
4. plt.figure(1, (8, 8))
5.  
6. ax = plt.subplot(2, 2, 1)
7. plt.plot(resultC2['time'], resultC2['[AcP]'],label='AcP') # Acetyl Phosphate
8. plt.plot(resultC2['time'], resultC2['[HOAc]'], label='HOAc') # Acetate
9.  
10. plt.legend()
11. plt.xlabel('Time, min')
12. plt.ylabel('Concentration, mM')
13.  
14. ax = plt.subplot(2, 2, 2)
15. plt.plot(resultC2['time'], resultC2['[ACoA]'],label='ACoA',color="red") # acetyl CoA
16. plt.legend()
17. plt.xlabel('Time, min')
18. plt.ylabel('Concentration, mM')
19.  
20. ax = plt.subplot(2, 2, 3)
21. plt.plot(resultC2['time'], resultC2['[Pta]'],label='Pta',color="purple") # Pta
22.  
23. plt.legend()
24. plt.xlabel('Time, min')
25. plt.ylabel('Concentration, mM')
26.  
27. ax = plt.subplot(2, 2, 4)
28. plt.plot(resultC2['time'], resultC2['[Acs]'],label='Acs',color="violet") # Acs
29. plt.plot(resultC2['time'], resultC2['[LacI]'], label='LacI',color="darkred") # LacI
30.  
31. plt.legend()
32. plt.xlabel('Time, min')
33. plt.ylabel('Concentration, mM')
34.  
35. plt.tight_layout()
36. plt.show()