contour plot join prior and posterior distribution of MCMC result with HPD

Can anyone contour plot like this -in attached image, using following prior distribution and the MCMC samples of two variables:

prior distribution: X~N(4,2),Y~N(4,2)

posterior distribution MCMC samples is attached

[attachment=219136,1102]

solid contour is posterior distribution density

dash contour is prior distribution density

number on the contour line is highest probability density(50%,90%)