[F(x):=,S(k,x):=k*SUM(p^(k-1),p,1,x-1)]

M(k,x,h):=INT(S(2*k,(y-x)/h)*DIF(F(y),y,2*k),y,x,x+h)/(2*k)!

;User=Simp(User)
S(2*1,x)=x*(x-1)

;Simp(User')
S(2*1,(y-x)/h)=(x-y)*(x-y+h)/h^2

;Sub(User)
M(1,x,h)=INT((x-y)*(x-y+h)/h^2*DIF(F(y),y,2*1),y,x,x+h)/(2*1)!

M(1,x,h)=INT((x-y)*(x-y+h)*F''(y),y,x,x+h)/(2*h^2)