Revised vNormalIntegral function. NN/2-1 should be multiple of 8.
Let (NN/2-1)/8 = LV, when LV=20, NN=322, which yields sufficient
precision for the integral. However, LV can be a tunable parameter
to accelerate the program. (Bigger LV may cause faster execution).