To derive let's write what #Var(XY)#:

#Var[XY] = E[(XY)^2] – {E[XY]}^2#

# Var(XY) = color(red)(E[X^2Y^2]) – color(blue)((E[X]E[Y])^2)#

#= color(red)(E[X^2]E[Y^2]) + Cov(X^2,Y^2)#

#= {color(blue)(E[X]E[Y]) + Cov(X,Y)}^2#

#=color(blue)({E[X]E[Y]}^2 + 2{E[X]E[Y]})Cov(X,Y) + Cov^2(X,Y)#

#Var(XY)=color(red)(E[X^2]E[Y^2]) +Cov(X^2,Y^2)- {color(blue)(E^2[X]E^2[Y]) + 2E[X]E[Y]Cov(X,Y) + Cov^2(X,Y)}#

Now if #X and Y# were independent the covariance will vanish which implies that correlation is also zero. However, in this case your random variables are correlated, thus the covariance stays on the above equation.

Now if you want to take it further you can sub the below for the expectation and variance operators. But i will stop here

#E[X]=mu_x; E[Y]=mu_y#

#E[XY]=E[X]E[Y]=mu_xmu_y#

#Var(X) = sigma_x; Var(Y)=sigma_y#

Assume normal distribution and you can show that:

#Var(XY)=mu_X^2*sigma_Y+ mu_Y^2*sigma_X+2*mu_Xmu_Y*Cov(X,Y)
+sigma_Xsigma_Y+Cov^2(X,Y)#

Please Double check, and good luck