read `MVNM.txt`:

K:=35;

print(`The list of the first`, K, `even diagonal moments of the trivariate normal distrubution with covariance matrix`):
print(``):
print(matrix([[1,c12,c13],[c12,1,c23],[c13,c23,1]]));
print(`are `):
print(``):

lprint([seq(MOM3(c12,c13,c23,[2*i,2*i,2*i]),i=1..K)]):

print(`-------------------------`):


print(`This took `, time(), `moments `):

quit: