This python script calculates generalised harmonic mean , H_{n,s} , i.e. 1+(1/2)^s+(1/3)^s+...+(1/n)^s  , with s = 1 for non-generalised