with（LREtools）：with（FormalPowerSeries）：#需要Maple 2022
M： =（1-z-sqrt（1-2*z-3*z^2））/
ogf:=（M*T^2）/（2*T-1）：需求：=查找RE（ogf，z，u（n））：
初始化：=[1，1，3，9，24，66]：iseq:=序列（u（i-1）=初始化[i]，i=1..nops（初始化））：
rmin:=sub（n=n-4，最小递归（req，u（n），{iseq}）[1]）；#马塔尔复发
a： =gfun:-rectproc（{rmin，iseq}，u（n），记住）：
seq（a（n），n=0..27）#乔治·菲舍尔2022年11月4日
ogf:=（1-x-sqrt（-3*x^2-2*x+1））/（2*x^2*（3*x^2+2*sqrt（-3*x^2-2*x+1）+2*x-1））：
a:=FindSeq（ogf）：序列（a（n），n=0..28）#彼得·卢什尼2022年11月4日
