(PARI)N=20;x='x+O('x^N);Vec(serlaplace(对数(总和(k=0,N,3^二项式(k,2)*x^k/k!)))\\Seiichi Manyama先生2019年5月18日
(岩浆)
m: =30;
f: =func<x|(&+[3^二项式(n,2)*x^n/阶乘(n):[0..m+3]]中的n)>;
R<x>:=PowerSeriesRing(基本原理(),m);
系数(R!(拉普拉斯(Log(f(x))))//G.C.格鲁贝尔2023年4月28日
(SageMath)
米=30
定义f(x):返回和(3^二项式(n,2)*x^n/范围(m+4)内n的阶乘(n))
P.<x>=PowerSeriesRing(QQ,prec)
return P(log(f(x))).egf_to_ogf().list()