replace internal MVCHNV FORTRAN FUNCTION with R's sqrt(qchisq(p, n, FALSE, FALSE). This fixes a problem where NaN was returned as reported by David Charles Airey
