macro
qq x.1-x.p
mconstant i n p t chis
mcolumn d x.1-x.p dd pi q ss tt
mmatrix s sinv ma mb mc md me
let n=count(x.1)
cova x.1-x.p s
invert s sinv
do i=1:p
  let x.i=x.i-mean(x.i)
enddo 
do i=1:n
  copy x.1-x.p ma;
    use i.
  transpose ma mb
  multiply ma sinv mc
  multiply mc mb md
  copy md tt
  let t=tt(1)
  let d(i)=t
enddo
set pi
  1:n
end
let pi=(pi-0.5)/n
sort d dd
invcdf pi q;
chis p.
plot q*dd
invcdf 0.5 chis;
chis p.
let ss=dd<chis
let t=sum(ss)/n
print t 
if t>0.5
  note distribusi data multinormal
endif
if t<=0.5
  note distribusi data bukan multinormal
endif
endmacro