function [itg,error]=iehardy(a,b,f)
c=(a+b)/2;
i=hardy(a,b,f);
i1=hardy(a,c,f);
i2=hardy(c,b,f);
itg=i1+i2;
error=(256/255)*abs(itg-i);