The density of a bivariate normal distribution with mean equal to (0,0) and correlation matrix equal to
Ļ = [1.0 0.5; 0.5 1.0]
is
f(x,y) = 1/(Ļā3) exp[-2(x^2 + y^2 - xy)/3].
Write R code that will create objects x and y that are each equal to the sequence {-3.00, -2.95, ..., 2.95, 3.00}. Then create a matrix z and populate it with the values of f(x,y) for each combination of x and y values. A nested for loop can be used to do this. Once this has been completed create a contour plot of the data. The lines on the plot should be red.