boxM function

Box's M-test