```run-python
import numpy as np
import matplotlib.pyplot as plt
def f(x,t):
x1, x2 = x
return [x1*(4-2*x1-x2), x2*(9-3*x1-3*x2)]
x1min, x1max, x1step = 0, 3, 0.25
x2min, x2max, x2step = 0, 3, 0.25
x1points = np.arange(x1min, x1max+x1step, x1step)
x2points = np.arange(x2min, x2max+x2step, x2step)
x1, x2 = np.meshgrid(x1points, x2points)
Df1, Df2 = f([x1,x2],0) #input for t doesn't matter since the equation is autonomous
fig, ax = plt.subplots(figsize=(9,9))
ax.quiver(x1,x2,Df1,Df2)
ax.axis([0,3,0,3])
ax.set_aspect('equal')
plt.show()
```