import matplotlib.pyplot as plt
from preliz import MvNormal
_, axes = plt.subplots(2, 2, figsize=(9, 9), sharex=True, sharey=True)
mus = [[0., 0], [3, -2], [0., 0], [0., 0]]
sigmas = [np.eye(2), np.eye(2), np.array([[2, 2], [2, 4]]), np.array([[2, -2], [-2, 4]])]
for mu, sigma, ax in zip(mus, sigmas, axes.ravel()):
    MvNormal(mu, sigma).plot_pdf(marginals=False, ax=ax)