Since we need to calculate the rates for two different values of $K_m$, we define a function to do it in the code below.
import pylab
# Michaelis-Menten constants (M), maximum rate(M.s-1)
Km1, Km2, Vmax = 0.04, 0.4, 0.01
# Grid of substrate concentration values (M)
S = pylab.linspace(0, 2, 1000)
def rate(S, Km, Vmax):
return Vmax * S / (Km + S)
pylab.plot(S, rate(S, Km1, Vmax))
pylab.plot(S, rate(S, Km2, Vmax))
pylab.show()