How to approximate a 2d distribution with mu and covariance in python

Question

I want to create a heatmap, depending on just the mu and the covariance. The background is, that I want to approximate a specific complex distribution with just these two parame