function  TRAPEZ=trap1(I,a,b)

if a>b    
error("b must be greater than a") 
end 
h=(b-a); 
TRAPEZ=h*(feval(I,a)+feval(I,b))/2; 


end