Program TabRGR;
{Tim Rubchenko, 04-113, 14.12.2010}
{Programma tabulirovanja funkcii Y=f(A,B,X), gde argument X izmenjaetsja ot   
nachalnogo Xn s shagom Dx prinimaet N razlich.
   znachenij, parametrA izmenjaetsja
   ot An do Ak s shagom Da,a parametr B
    chislenno raven integralu f-cii
     ln(x)*cos(x),vychisl.prizazdannyhznachenijah c i d i pogreshnosti Eps}
const Nmax=20;
        type
        Tmy=Array[1..Nmax,1..Nmax] of real;
        Tmx=Array[1..Nmax] of real;
        Ter=Array[1..Nmax,1..Nmax] of integer;
Var
Ma,Mx:Tmx;
My:Tmy;
Xn,Dx,B,An,Ak,Da,C,D,Eps,Int:real;

K,N,Er,Km:integer;
Mer:Ter;

{Text of function}
Function F(X:real): real;
begin
   F:=ln(x)*cos(x);
   end;

   {Text procedure of integrate}
Procedure Integrate(C,D,Eps:real; Km:Integer; Var Int:Real; Var Er:Integer);
Var
     X,Int1,Y1,Yn,S4,S2,H,Del:real;
     K,I,J:Integer;
begin
       In1:=0;
       Y1:=F(c);
       Yn:=F(d);
       K:=2;
       Er:=1;
       J:=1;
       While J<=Km do
       begin
            In:=0;
            H:=(d-c)/k;
            S4:=0;
            S2:=0;
            I:=1;
            X:=c+H;
              
       repeat
             If odd(I) then S4:=S4+F(x)
             else S2:=S2+F(x);
             Int:=H*(Y1+Yn+4*S4+2*S2)/3;
             Inc(I);
             X:=X+H;      
       Until I>=K;
             del:=Abs(Int1-Int);
             If del<Eps then begin
             Er:=Er-1;
             J:=Km;
             Exit;
             end
      else                
   Begin           
            Int1:=Int;
            k:=k+2;
            J:=J+1;
            end;
            end;
       end;

{Text procedure of tabulation}
Procedure Tab(B,Xn,Dx,An,Ak,Da:real; N:integer;  var Ma,Mx:Tmx; var My:Tmy; var Mer:Ter);
var
     I,J:Integer;
     A,X,Y:Real;
     begin
       I:=1;
          A:=An;
          While A<=Ak do begin
           
          X:=Xn;
          For J:=1 to N do begin
          Mer[I,J]:=0;
           if b<>0 then   
             
             
           Y:=A*sqr(x)+B*x+cos(3.14*(A/B))
           else Mer[I,J]:=1;
             
             
           X:=Xn+Dx;
           Mx[J]:=X;
           end;
             

     My[I,J]:=Y;
Ma[I]:=A;
A:=A+Da
end;
end;

{Text procedure ResOut}
Procedure ResOut(Var Mx:Tmx; var Ma:Tmx; Var My:Tmy; Var Mer:Ter; K:integer);
Var
     I,J:integer;
begin
       for I:=1 to N do {Vyvod result vypolnenija programmy}
        begin
             writeln('A[',I,']=',Ma[I]:8:4);
             writeln('   X','          Y');
              for J:=1 to K do
              If Mer[I,J]=1 then
               writeln(Mx[I]:6:3, ' ERROR B DAHHblX')
               else
               writeln(Mx[I]:6:3, '    ',My[I,J]:6:3);
        end;
end;

begin {Start of general program}
   writeln('BBEDuTE 3HA4EHu9I nEPEMEHHbIX:Xn,Dx,N,Ak,An,Da');
   readln(Xn,Dx,N,Ak,An,Da);
Writeln('Xn=',Xn:4:2,' Dx=',Dx:4:2,' N=',N, ' Ak=',Ak,' An=',An,' Da=',Da);
writeln('BBEDuTE nPEDEJIbI uHTerpupoBaHu9 C,D,norpewHocTb Eps',
   ', npeDeJIbHoe 4ucJIo u,ukJIoB Km');
readln(C,D,Eps,Km);
writeln('C=',C:6:3,' D=',D:6:3,' Eps=',Eps:7:5,' Km=',Km:3);
   Integrate(C,D,Eps,Km,Int,Er); {Vyzov programmy nahozhdenija integrala}
if Er=1 then
     begin
          writeln('UHTErPAJl HE HAUDEH 3a ', Km:2,' iteracij');
          exit
     end;
writeln('Integral PABEH B=',B:4:2);
Tab(B,Xn,Dx,An,Ak,Da,N,Ma,Mx,My,Mer); {Vyzov podprogrammy tabulirovania}
K:=trunc((Ak-An)/Da+1);
ResOut(Mx,Ma,My,Mer,K);  {Vyzov podprogrammyvyvoda dannyh}
readln;
end.