macro
input x.1-x.jmlcol n

mconstant i j n jcol totskor luas
mcolumn di p qq rata2 x.1-x.jmlcol y.1-y.n skor
mmatrix vek vektran vmean mat cov hasil

Let n=count(x.1)
copy x.1-x.jmlcol mat
trans mat mat 
cova x.1-x.jmlcol cov
invers cov cov
copy mat y.1-y.n
Rmean y.1-y.n rata2 
copy rata2 vmean
Note macro is running......
do i=1:n
 copy y.i vek
 subtract vek rata2 vek
 trans vek vektran 
 multi vektran cov hasil
 multi hasil vek hasil
 copy hasil y.i
endDo

copy y.1-y.n mat
trans mat mat
copy mat di

do i=1:n 
 let p(i)=(i-1/2)/n
enddo

InvCdf p qq;
 chisquare jmlcol.
let di=sort(di)

plot qq*di;
Title " Plot Uji Multinormal".

do i=1:n 
 if qq(i)<=di(i)
  let skor(i)=1
 else
  Let skor(i)=0 
 endif
enddo

let totskor=sum(skor)
let luas=totskor/n

Note Ho : data berdistribusi multinormal
Note H1 : data tidak berdistribusi multinormal

# name luas='daerah dibawah kurva chi-square (%)='
print luas
if luas <= 0.5
 Note Tolak Ho 
 Note Dan dapat disimpulkan bahwa data tidak mengikuti distribusi multinormal
else
 Note gagal tolak Ho
 Note dan dapat disimpulkan bahwa data mengikuti distribusi multinormal
endIf

endMacro.