function y=update(f,a,h,n,val)
%   sumx=0;
%    for i=1:n
%      x=a+(i-0.5)*h;
%     sumx=sumx + h*feval(f,x);
%   end;
    v=a+0.5*h:h:a+(n-0.5)*h;
    sumx=h*sum(feval(f,v));
    y=0.5*(val+sumx);