(PARI)我的(x='x+O('x^40));Vec((1-x)/(1-5*x-2*x^2+8*x^3))\\G.C.格鲁贝尔2019年10月27日
(Magma)R<x>:=PowerSeriesRing(整数(),40);系数(R!((1-x)/(1-5*x-2*x^2+8*x^3))//G.C.格鲁贝尔2019年10月27日
(鼠尾草)
P.<x>=PowerSeriesRing(ZZ,prec)
返回P((1-x)/(1-5*x-2*x^2+8*x^3)).list()
(间隙)a:=[1,4,22];;对于[4..40]中的n,做a[n]:=5*a[n-1]+2*a[n-2]-8*a[n-3];od;a#G.C.格鲁贝尔2019年10月27日