Metamath Proof Explorer


Theorem bpolyval

Description: The value of the Bernoulli polynomials. (Contributed by Scott Fenton, 16-May-2014)

Ref Expression
Assertion bpolyval ( ( 𝑁 ∈ ℕ0𝑋 ∈ ℂ ) → ( 𝑁 BernPoly 𝑋 ) = ( ( 𝑋𝑁 ) − Σ 𝑘 ∈ ( 0 ... ( 𝑁 − 1 ) ) ( ( 𝑁 C 𝑘 ) · ( ( 𝑘 BernPoly 𝑋 ) / ( ( 𝑁𝑘 ) + 1 ) ) ) ) )

Proof

Step Hyp Ref Expression
1 fvex ( ♯ ‘ dom 𝑐 ) ∈ V
2 oveq2 ( 𝑛 = ( ♯ ‘ dom 𝑐 ) → ( 𝑋𝑛 ) = ( 𝑋 ↑ ( ♯ ‘ dom 𝑐 ) ) )
3 oveq1 ( 𝑛 = ( ♯ ‘ dom 𝑐 ) → ( 𝑛 C 𝑚 ) = ( ( ♯ ‘ dom 𝑐 ) C 𝑚 ) )
4 oveq1 ( 𝑛 = ( ♯ ‘ dom 𝑐 ) → ( 𝑛𝑚 ) = ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) )
5 4 oveq1d ( 𝑛 = ( ♯ ‘ dom 𝑐 ) → ( ( 𝑛𝑚 ) + 1 ) = ( ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) + 1 ) )
6 5 oveq2d ( 𝑛 = ( ♯ ‘ dom 𝑐 ) → ( ( 𝑐𝑚 ) / ( ( 𝑛𝑚 ) + 1 ) ) = ( ( 𝑐𝑚 ) / ( ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) + 1 ) ) )
7 3 6 oveq12d ( 𝑛 = ( ♯ ‘ dom 𝑐 ) → ( ( 𝑛 C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( 𝑛𝑚 ) + 1 ) ) ) = ( ( ( ♯ ‘ dom 𝑐 ) C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) + 1 ) ) ) )
8 7 sumeq2sdv ( 𝑛 = ( ♯ ‘ dom 𝑐 ) → Σ 𝑚 ∈ dom 𝑐 ( ( 𝑛 C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( 𝑛𝑚 ) + 1 ) ) ) = Σ 𝑚 ∈ dom 𝑐 ( ( ( ♯ ‘ dom 𝑐 ) C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) + 1 ) ) ) )
9 2 8 oveq12d ( 𝑛 = ( ♯ ‘ dom 𝑐 ) → ( ( 𝑋𝑛 ) − Σ 𝑚 ∈ dom 𝑐 ( ( 𝑛 C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( 𝑛𝑚 ) + 1 ) ) ) ) = ( ( 𝑋 ↑ ( ♯ ‘ dom 𝑐 ) ) − Σ 𝑚 ∈ dom 𝑐 ( ( ( ♯ ‘ dom 𝑐 ) C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) + 1 ) ) ) ) )
10 1 9 csbie ( ♯ ‘ dom 𝑐 ) / 𝑛 ( ( 𝑋𝑛 ) − Σ 𝑚 ∈ dom 𝑐 ( ( 𝑛 C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( 𝑛𝑚 ) + 1 ) ) ) ) = ( ( 𝑋 ↑ ( ♯ ‘ dom 𝑐 ) ) − Σ 𝑚 ∈ dom 𝑐 ( ( ( ♯ ‘ dom 𝑐 ) C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) + 1 ) ) ) )
11 oveq2 ( 𝑚 = 𝑘 → ( 𝑛 C 𝑚 ) = ( 𝑛 C 𝑘 ) )
12 fveq2 ( 𝑚 = 𝑘 → ( 𝑐𝑚 ) = ( 𝑐𝑘 ) )
13 oveq2 ( 𝑚 = 𝑘 → ( 𝑛𝑚 ) = ( 𝑛𝑘 ) )
14 13 oveq1d ( 𝑚 = 𝑘 → ( ( 𝑛𝑚 ) + 1 ) = ( ( 𝑛𝑘 ) + 1 ) )
15 12 14 oveq12d ( 𝑚 = 𝑘 → ( ( 𝑐𝑚 ) / ( ( 𝑛𝑚 ) + 1 ) ) = ( ( 𝑐𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) )
16 11 15 oveq12d ( 𝑚 = 𝑘 → ( ( 𝑛 C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( 𝑛𝑚 ) + 1 ) ) ) = ( ( 𝑛 C 𝑘 ) · ( ( 𝑐𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) )
17 16 cbvsumv Σ 𝑚 ∈ dom 𝑐 ( ( 𝑛 C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( 𝑛𝑚 ) + 1 ) ) ) = Σ 𝑘 ∈ dom 𝑐 ( ( 𝑛 C 𝑘 ) · ( ( 𝑐𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) )
18 dmeq ( 𝑐 = 𝑔 → dom 𝑐 = dom 𝑔 )
19 fveq1 ( 𝑐 = 𝑔 → ( 𝑐𝑘 ) = ( 𝑔𝑘 ) )
20 19 oveq1d ( 𝑐 = 𝑔 → ( ( 𝑐𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) = ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) )
21 20 oveq2d ( 𝑐 = 𝑔 → ( ( 𝑛 C 𝑘 ) · ( ( 𝑐𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) = ( ( 𝑛 C 𝑘 ) · ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) )
22 21 adantr ( ( 𝑐 = 𝑔𝑘 ∈ dom 𝑐 ) → ( ( 𝑛 C 𝑘 ) · ( ( 𝑐𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) = ( ( 𝑛 C 𝑘 ) · ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) )
23 18 22 sumeq12dv ( 𝑐 = 𝑔 → Σ 𝑘 ∈ dom 𝑐 ( ( 𝑛 C 𝑘 ) · ( ( 𝑐𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) = Σ 𝑘 ∈ dom 𝑔 ( ( 𝑛 C 𝑘 ) · ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) )
24 17 23 eqtrid ( 𝑐 = 𝑔 → Σ 𝑚 ∈ dom 𝑐 ( ( 𝑛 C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( 𝑛𝑚 ) + 1 ) ) ) = Σ 𝑘 ∈ dom 𝑔 ( ( 𝑛 C 𝑘 ) · ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) )
25 24 oveq2d ( 𝑐 = 𝑔 → ( ( 𝑋𝑛 ) − Σ 𝑚 ∈ dom 𝑐 ( ( 𝑛 C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( 𝑛𝑚 ) + 1 ) ) ) ) = ( ( 𝑋𝑛 ) − Σ 𝑘 ∈ dom 𝑔 ( ( 𝑛 C 𝑘 ) · ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) ) )
26 25 csbeq2dv ( 𝑐 = 𝑔 ( ♯ ‘ dom 𝑐 ) / 𝑛 ( ( 𝑋𝑛 ) − Σ 𝑚 ∈ dom 𝑐 ( ( 𝑛 C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( 𝑛𝑚 ) + 1 ) ) ) ) = ( ♯ ‘ dom 𝑐 ) / 𝑛 ( ( 𝑋𝑛 ) − Σ 𝑘 ∈ dom 𝑔 ( ( 𝑛 C 𝑘 ) · ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) ) )
27 10 26 eqtr3id ( 𝑐 = 𝑔 → ( ( 𝑋 ↑ ( ♯ ‘ dom 𝑐 ) ) − Σ 𝑚 ∈ dom 𝑐 ( ( ( ♯ ‘ dom 𝑐 ) C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) + 1 ) ) ) ) = ( ♯ ‘ dom 𝑐 ) / 𝑛 ( ( 𝑋𝑛 ) − Σ 𝑘 ∈ dom 𝑔 ( ( 𝑛 C 𝑘 ) · ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) ) )
28 18 fveq2d ( 𝑐 = 𝑔 → ( ♯ ‘ dom 𝑐 ) = ( ♯ ‘ dom 𝑔 ) )
29 28 csbeq1d ( 𝑐 = 𝑔 ( ♯ ‘ dom 𝑐 ) / 𝑛 ( ( 𝑋𝑛 ) − Σ 𝑘 ∈ dom 𝑔 ( ( 𝑛 C 𝑘 ) · ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) ) = ( ♯ ‘ dom 𝑔 ) / 𝑛 ( ( 𝑋𝑛 ) − Σ 𝑘 ∈ dom 𝑔 ( ( 𝑛 C 𝑘 ) · ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) ) )
30 27 29 eqtrd ( 𝑐 = 𝑔 → ( ( 𝑋 ↑ ( ♯ ‘ dom 𝑐 ) ) − Σ 𝑚 ∈ dom 𝑐 ( ( ( ♯ ‘ dom 𝑐 ) C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) + 1 ) ) ) ) = ( ♯ ‘ dom 𝑔 ) / 𝑛 ( ( 𝑋𝑛 ) − Σ 𝑘 ∈ dom 𝑔 ( ( 𝑛 C 𝑘 ) · ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) ) )
31 30 cbvmptv ( 𝑐 ∈ V ↦ ( ( 𝑋 ↑ ( ♯ ‘ dom 𝑐 ) ) − Σ 𝑚 ∈ dom 𝑐 ( ( ( ♯ ‘ dom 𝑐 ) C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) + 1 ) ) ) ) ) = ( 𝑔 ∈ V ↦ ( ♯ ‘ dom 𝑔 ) / 𝑛 ( ( 𝑋𝑛 ) − Σ 𝑘 ∈ dom 𝑔 ( ( 𝑛 C 𝑘 ) · ( ( 𝑔𝑘 ) / ( ( 𝑛𝑘 ) + 1 ) ) ) ) )
32 eqid wrecs ( < , ℕ0 , ( 𝑐 ∈ V ↦ ( ( 𝑋 ↑ ( ♯ ‘ dom 𝑐 ) ) − Σ 𝑚 ∈ dom 𝑐 ( ( ( ♯ ‘ dom 𝑐 ) C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) + 1 ) ) ) ) ) ) = wrecs ( < , ℕ0 , ( 𝑐 ∈ V ↦ ( ( 𝑋 ↑ ( ♯ ‘ dom 𝑐 ) ) − Σ 𝑚 ∈ dom 𝑐 ( ( ( ♯ ‘ dom 𝑐 ) C 𝑚 ) · ( ( 𝑐𝑚 ) / ( ( ( ♯ ‘ dom 𝑐 ) − 𝑚 ) + 1 ) ) ) ) ) )
33 31 32 bpolylem ( ( 𝑁 ∈ ℕ0𝑋 ∈ ℂ ) → ( 𝑁 BernPoly 𝑋 ) = ( ( 𝑋𝑁 ) − Σ 𝑘 ∈ ( 0 ... ( 𝑁 − 1 ) ) ( ( 𝑁 C 𝑘 ) · ( ( 𝑘 BernPoly 𝑋 ) / ( ( 𝑁𝑘 ) + 1 ) ) ) ) )