help simpson
  returns an approximation for the integral of f(x) over [a,b]
  using Simpson's rule



simpson('cos',pi/4,pi/2)

ans =

    0.5859

diary off;