Some new GP features hypergeom hypergeom now supports power series: ? hypergeom([1,1,1],[5/3,5/3],z+O(z^5)) %57 = 1+9/25*z+81/400*z^2+6561/48400*z^3+59049/59290 ? hypergeom([1,1,1],[5/3,5/3],1/2) %58 = 1.2579176749481097328724460713125550477 ? hypergeom([1,1,1],[5/3,5/3],1/2+z+O(z^5)) %59 = 1.257917674948109733+0.7611467927103072052*z % +0.7887370510268457097*z^2+1.007980535530049765*z^ % +1.437897250391074849*z^4+O(z^5)