|
R[n_]：=b加[1，总计[b^#R[#，q]&/@范围[n]]]
Vp[n_]：=总计[2 v[#+2]q^（#+2）和/@范围[n]]
H[n_]：=展开[1/2*r^2+Vp[n]]
RRules[n_]：=与[{H=系列[ReplaceAll[H[n]，{q->R[n]q，R->R[n]}]，{b，0，n+2}]}，函数[{rules}，
嵌套[Rule[#[[1]]，ReplaceAll[#[2]]，rules]]&/@#&，rules，n]][
展平[R[#，q]->展开[-ReplaceAll[系数[H，b^（#+2）]，{R[#、q]->0}]&/@范围[n]]]
xDot[n_]：=展开[正常@系列[替换全部[Q^2 D[D[Q[t]，t]/Q[t]t，t]，{D[Q[t]，t]->R[n]P，Q[t]->R[n]Q，R->R[n]，D[Q[t]，{t，2}]
->替换全部[D[-（q^2/2+Vp[n]），q]，q->R[n]q]}]，{b，0，n}]/。R规则[n]/。{P^2->1-Q^2}]
dt[n_]：=展开[正常@系列[1/x点[n]，{b，0，n}]]
基[n_]：=次数[Times@@（v/@#），Q^Total[#]，2]和/@（整数分区[n]/.x_Integer:>x+2）
三角形行[n_，fun_]：=系数[fun，b^n#]和/@basis[n]
带有[{dt10=dt[10]}，TriangleRow[#，dt10]/。v[_]->0&&@范围[10]]
| 