function I=Trapeze(f,a,b)
  I=((b-a)/2)*(f(a)+f(b));