Given $X_i \sim N(\mu_i, \sigma_i^2)$, using the FOSM approach, SOSM approach, and one coding approach (MatLAB, Excel, R, etc), find the first and second moments of
where
$X = X_1^2 + X_2^2 - \frac{2}{X_1 X_2}$
$X_1 \sim N(16, 3.4^2)$ and $X_2 \sim N(21, 1.8^2)$.