bq:=proc（q:：integer）局部m；如果q mod 2=0，则返回（0）；否则m：=（q-1）/2；返回值（总和（二项式（q，s）*二项式式（s+m-1，m），s=0..q）/（q*2^（m+1）））；fi；结束：H:=proc（maxord:：integer）局部结果，r，B；结果：=0；B:=0；对于r从2到maxord乘2做B:=B+bq（r/2）*y^r；od：对于从1到maxord/2的r，执行结果：=结果+B^r；od：返回（结果）；end:maxord:=60:Hser：=展开（H（maxord））：对于n，从2到maxord乘以2，进行简化（coefftayl（Hser，y=0，n））；od#R.J.马塔尔2006年4月24日
