#这是一个粗略的程序，假设上述猜测是正确的（后来在Hadjicostas（2019）中得到了证明）。它仅对n>=3有效，因为Bower对DHK[k]的特殊定义适用于k=1和k=2的情况。
使用（数字理论）；
ff:=proc（k，x）（1/2）*x^k*（加法（φ（n）/（1-x^n）^（k/n），n以除数（k）表示）/k-（x+1）/（1-x^2）^楼层（（1/2）*k+1））；终末程序；
=proc（k，x）添加（Moebius（d）*ff（k/d，x^d），d在除法（k）中）；终末程序；
vv:=proc（n）简化（sub（x=0，diff（gg（n，x），x$（2*n）））/阶乘（2*n））；终末程序；
