Example

> f:=x->3+sin(2*x+.5)/2;ddddf:=(D@@4)(f);

> for i from 0 to 5 do h:=2.^(-i): int1:=int(f(x),x=1..1.+2*h):int2:=(h/3)*(f(1.+2*h)+4*f(1.+h)+f(1.)):print(h, int2,abs(int2-int1),8.*h^5/90); od:

>