The binormal probability density function

PDF[MultinormalDistribution[{20,10},{{25,9},{9,9}}],{x,y}]

ContourPlot[ PDF[MultinormalDistribution[{20,10},{{25,9},{9,9}}],{x,y}],{x,10,30},{y,5,15}]

Manipulate[
 Module[{v},v={{σx^2,ρ σx σy},{ρ σx σy,σy^2}};ContourPlot[ PDF[MultinormalDistribution[{μx,μy},v],{x,y}],{x,μx-3 5,μx+3 5},{y,μy-3 3,μy+3 3}]],
{{μx,20},0,50},{{μy,10},0,50},
{{σx,5},0.5,10},{{σy,3},0.5,10},
Style["Notice that note only ρ, but also σx and σy affect the inclination of main axes - keep an eye on the coordinate scales",12],
{{ρ,0.6},-0.999,0.999,Appearance->"Labeled"}]