рисование двумерных гауссовских распределений в matplotlib

Как мы можем построить (в python matplotlib) двумерные распределения Гаусса, учитывая их центры и ковариационные матрицы в виде пустых массивов?

Допустим, что наши параметры следующие:

center1=np.array([3,3])
center2=np.array([5,5])
cov1=np.array([ [1.,.5], [.5,.1]])
cov2=np.array([ [.2,.5], [.5,.2]])