read `QuickSortAnalysis.txt`:

t0:=time():

M:=10000:
K:=50000:
n1:=200:

m:=4:

print(`Simulation Average, Variance and Higher Moments vs. their exact values`):

gu:=MT(M,K,n1,m):

print(`This is a simulation of Quicksort on random lists of length`, n1):
print(`using`, K, `randomly chosen sequence of that length`):

print(``):
print(`The simulation average, followed by the theoretical average are`, gu[1][1],gu[2][1]):
print(``):
print(`The simulation variance, followed by the theoretical variance are`, gu[1][2],gu[2][2]):

for m1 from 3 to m do

print(`The simulation scaled`, m1, `-th moment  followed by the theoretical value are`, gu[1][m1],gu[2][m1]):

od:

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

print(`This took`, time()-t0, `seconds. `):

quit: