Let
F(x)=k=0∑∞(4k+1)(4k+2)(4k+3)(4k+4)x4k+4.
This power series converges for ∣x∣≤1 and our goal is to
compute F(1).
Differentiating 4 times, we get
F(IV)(x)=k=0∑∞x4k=1−x41.
Since F(0)=F′(0)=F′′(0)=F′′′(0)=0 and F is continuous at
1−0 by Abel's continuity theorem, integrating 4 times we get
F′′′(y)=F′′′(0)+∫0yF(IV)(x)dx=∫0y1−x4dx=21arctany+41log(1+y)−41log(1−y),
F′′(z)=F′′(0)+∫0zF′′′(y)dx=∫0z(21arctany+41log(1+y)−41log(1−y))dy==21(zarctanz−∫0z1+y2ydy)+41((1+z)log(1+z)−∫0zdy)+41((1−z)log(1−z)+∫0zdy)==21zarctanz−41log(1+z2)+41(1+z)log(1+z)+41(1−z)log(1−z),
F′(t)=∫0t(21zarctanz−41log(1+z2)+41(1+z)log(1+z)+41(1−z)log(1−z))dt==41((1+t2)arctant−t)−41(tlog(1+t2)−2t+2arctant)++81((1+t)2log(1+t)−t−21t2)−81(−(1−t)2log(1−t)+t−21t2)==41(−1+t2)arctant−41tlog(1+t2)+81(1+t)2log(1+t)−81(1−t)2log(1−t),
F(1)=∫01(41(−1+t2)arctant−41tlog(1+t2)+81(1+t)2log(1+t)−81(1−t)2log(1−t))dt==[12−3t+t3arctant+241−3t2log(1+t2)+24(1+t)3log(1+t)+24(1−t)3log(1−t)]01=4ln2−24π.
Remark. The computation can be shorter if we change the order of
integrations.
F(1)=∫t=01∫z=0t∫y=0z∫x=0y1−x41dxdydzdt=∫x=011−x41(∫y=x1∫z=y1∫t=z1dtdzdy)dx==∫x=011−x41⋅6(1−x)3dx=[−61arctanx−121log(1+x2)+31log(1+x)]01=4ln2−24π.