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

[U(y):=F(y),V(y):=y]

INT(V(y)*DIF(U(y),y),y,x,x+h)=V(x+h)*U(x+h)-V(x)*U(x)-INT(U(y)*DIF(V(y),y),y,~
x,x+h)

;Simp(#3)
INT(y*F'(y),y,x,x+h)=(x+h)*F(x+h)-x*F(x)-INT(F(y),y,x,x+h)

INT((x-y)*(x-y+h)*F''(y),y,x,x+h)=(2*x+h)*F(x+h)-(2*x+h)*F(x)-2*INT(y*F'(y),y~
,x,x+h)

;Sub(#5')
INT((x-y)*(x-y+h)*F''(y),y,x,x+h)=(2*x+h)*F(x+h)-(2*x+h)*F(x)-2*((x+h)*F(x+h)~
-x*F(x)-INT(F(y),y,x,x+h))

;Simp(#6')
INT((x-y)*(x-y+h)*F''(y),y,x,x+h)=-h*F(x+h)-h*F(x)+2*INT(F(y),y,x,x+h)

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

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

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

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