read(`ExpSP.txt`):

r:=2:
a:=2*r+2:
b:=2*r+14:
f:=[seq(GFsimu(Hnr(n,r),t,1000),n=a..b)];     # generating functions
av:=[seq(subs(t=1,diff(f[n],t)),n=1..b-a)];   # averages