#REPLACE NAMES OF THE VARIABLES USED HERE TO ONES YOU USED IN THE SIMULATION
modelC2.resetAll()
resultC2=modelC2.simulate(start=0, end=1000, points=10000)
plt.figure(1, (8, 8))

ax = plt.subplot(2, 2, 1)
plt.plot(resultC2['time'], resultC2['[AcP]'],label='AcP') # Acetyl Phosphate
plt.plot(resultC2['time'], resultC2['[HOAc]'], label='HOAc') # Acetate

plt.legend()
plt.xlabel('Time, min')
plt.ylabel('Concentration, mM')

ax = plt.subplot(2, 2, 2)
plt.plot(resultC2['time'], resultC2['[ACoA]'],label='ACoA',color="red") # acetyl CoA
plt.legend()
plt.xlabel('Time, min')
plt.ylabel('Concentration, mM')

ax = plt.subplot(2, 2, 3)
plt.plot(resultC2['time'], resultC2['[Pta]'],label='Pta',color="purple") # Pta

plt.legend()
plt.xlabel('Time, min')
plt.ylabel('Concentration, mM')

ax = plt.subplot(2, 2, 4)
plt.plot(resultC2['time'], resultC2['[Acs]'],label='Acs',color="violet") # Acs
plt.plot(resultC2['time'], resultC2['[LacI]'], label='LacI',color="darkred") # LacI

plt.legend()
plt.xlabel('Time, min')
plt.ylabel('Concentration, mM')

plt.tight_layout()
plt.show()