Hola
Dadas dos v.a \( X_1,X_2 \) con distribución conjunta continua, se cumple que son independientes si y sólo si se cumple que:
\( f(x_1,x_2)=f_1(x_1)f_2(x_2) \)
donde \( f,f_1,f_2 \) son respectivamente las funciones de densidad de \( (X_1,X_2),X_1 y X_2 \). Esto es inmediato de la definición de variables independientes (e incluso a veces para variables continuas se da como definición).
Pues bien simplemente usando la expresión de la densidad de una distribución normal bivariada comprueba que si la matriz de covarianzas es diagonal, tal densidad se expresa como producto de las dos densidades marginales de cada variable. Es bastante inmediato esencialmente usando que \( e^{a+b}=e^ae^b \).
Saludos.