gausscum:Stream Integer := concat([0,1], repeating [0]) gauss:=distributionByClassicalCumulants gausscum gaussmom k == odd? k => return 0 k2:= exquo(k,2) return (factorial(k)/(2^k2*factorial k2))::Integer gauss:=distributionByMoments [gaussmom k for k in 1..] gaussexp:= series(exp(z^2/2), z=0) coefficients gaussexp gaussmom := [c*factorial n for n in 1.. for c in coefficients gaussexp] R:= UP(z,Integer) W:=LODO1 R z:R:='z d:W:= D() d*z-z*d g6:=(d+z)^6 coefficient(g6, 0) coefficient(coefficient(g6, 0), 0) g:= distributionByMoments([coefficient(coefficient((d+z)^k, 0), 0) for k in 1..]) an:Stream Integer := repeating [0] bn:Stream Integer := [k for k in 1..] gauss:=distributionByJacobiParameters(an, bn) gauss:=gaussianDistribution 1 levymom:= rest rest freeCumulants gauss; levy:= distributionByMoments levymom hankelDeterminants(levy) jacobiParameters levy