add fortran code for computing 2nd moment of mvn