k:=2;
m:=4;
for j from 0 to m/2 do
  rule:=FiniteDifferenceFormula(m,k,j);
  err:=taylor(rule-(D@@k)(f)(0),h=0,m+3);
end do;