[impersonalis]

inter:=0
repeat   
         rez:=rez+inter; //<=========
         xi:=dx*i;  
         xi1:=dx*(i+1);  
         h:=(xi1-xi)/4;  
         y0:=y(xi);  
         y1:=y(xi+h);  
         y2:=y(xi+2*h);  
         y3:=y(xi+3*h);  
         y4:=y(xi1);  
         inter:=(2*h/45)*(7*y0+32*y1+12*y2+32*y3+7*y4);  
         //тут нельзя складывать - вдруг прирост уже и так меньше? тогда мы получим приближение лучше, чем заказывали
    until(Abs(inter) <= E); 


тогда