（PARI）seq（n）={my（A=O（x*x^n））；Vec（serlaplace（log（总和（k=0，n，3^二项式（k，2）*x^k/k！）+A）+log（cos（x+A）+sin（x+A））））/2）}\\安德鲁·霍罗伊德2018年9月10日
（岩浆）
m： =30；
f： =func<x|（&+[3^二项式（n，2）*x^n/阶乘（n）：[0..m+3]]中的n）>；
R<x>：=PowerSeriesRing（基本原理（），m）；
系数（R！（拉普拉斯（Log（f（x）*（Cos（x）+Sin（x）））/2））//G.C.格鲁贝尔2023年4月29日
（SageMath）
米=30
定义f（x）：返回和（3^二项式（n，2）*x^n/范围（m+4）内n的阶乘（n））
P.<x>=PowerSeriesRing（QQ，prec）
返回P（log（f（x）*（cos（x）+sin（x）））/2）.egf_to_ogf（）.list（）
