function I = trapez(a,b,fa,fb)
  I = (1/2)* ( fa + fb) * (b-a);
end