function area=trapes1(f,a,b)
ya=feval(f,a);
yb=feval(f,b);
area=(b-a)*(ya+yb)/2;