Consideremos el vector aleatorio normalmente distribuido $$X \sim \mathcal{N}(\mu, \Sigma)$$ ¿Cuál es la distribución de $Y = f(X)$ ?
En general $f$ este es un problema desafiante pero para el caso lineal afín $$f(x)_i = c_i + L_{ij}x_j$$ con $c$ un vector y $L$ una matriz. Sabemos que esto tiene una bonita forma cerrada. De hecho, $Y$ se distribuye de nuevo como una normal multivariante. $$Y \sim \mathcal{N}(c + L\mu, \;L\Sigma L^T)$$
Consideremos ahora el siguiente caso más sencillo. Consideremos que $f$ no es lineal, sino cuadrática. Es decir $$f(x)_i = c_i + L_{ij}x_j + H_{ijk} x_jx_k$$ con $c$ un vector, $L$ una matriz y $H$ un tensor de rango tres.
¿Existe una expresión de forma cerrada para la densidad de $Y$ ?