matplotlib:plt.subplots

import numpy as np
import matplotlib
matplotlib.use(“TkAgg”)
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.gridspec as gridspec

f,ax=plt.subplots(2,2,sharex=True,sharey=True)
ax[0,0].scatter([1,2],[1,2])
ax[0,1].scatter([1,2],[1,2])
ax[1,0].scatter([1,2],[1,2])
ax[1,1].scatter([1,2],[1,2])

plt.show()