I am trying to find the covariance of a multivariate normal on some interval [a,b] via integration. The function is "int_a^b int_a^b f(x, y) (x - mu_X) (y - mu_Y) dx dy&quo
