Edit to use LU decomposition and not exact 1st 2 moments