выч методы не разобраны! / Лаба№3 / Мишины / RELOAD

program lab2;
uses crt;
const n=4;a=1.2;b=0.4;
t:array[0..n] of real=(0.9062,0.5385,0,-0.5385,-0.9062);
c:array[0..n] of real=(0.2369,0.4786,0.5689,0.4786,0.2369);
var i,j:integer;
    s,bank,Q,P,H,H1,hs,hch,Eps1,hh1,Q1,Eps,Eps2,H2,Q2,Eps3,x,y,z:real;

  {¤ ­­ п дг­ЄжЁп}
  function f(x:real):real;
 begin
   f:=sqrt(x)*cos(sqr(x));
 end;

  {д®а¬г«  ѓ гбб }
  function gauss(dawn,up:real):real;
    begin
     s:=0;
     for i:=0 to n do
     s:=s+c[i]*f((dawn+up)/2+(up-dawn)/2*t[i]);
     s:=(up-dawn)/2*s;
     gauss:=s;
    end;

  {д®а¬г«  ва ЇҐжЁ©}
  function trap(dawn,up:real):real;
  var h1:real;
    begin
     h1:=(up-dawn);
     trap:=h1/6*(f(dawn)+f(up)+4*f(dawn+h1));
    end;

  {Џ®ЁбЄ Ї®¤е®¤п饣® ®в१Є }
procedure step(lim1,lim2:real;v:integer; var z:real);
  begin
   if v=1 then
       repeat
        h:=lim2-lim1;
        q:=trap(lim1,lim2);
        p:=trap(lim1,lim1+h/2)+trap(lim1+h/2,lim2);
        z:=lim2;eps1:=abs(p-q)/31;
        lim2:=lim1+(lim2-lim1)/2;
        until eps1<=(h/(b-a)*Eps)
         else repeat
        h:=lim2-lim1;
        q:=gauss(lim1,lim2);
        p:=gauss(lim1,lim1+h/2)+gauss(lim1+h/2,lim2);
        z:=lim2;eps1:=abs(p-q)/31;
        lim2:=lim1+(lim2-lim1)/2;
        until eps1<=(h/(b-a)*Eps);
  end;

  Begin
   clrscr;
   x:=a;y:=b;i:=0;
   H1:=0;hs:=0;Q1:=0;
   Eps2:=0;bank:=0;
   write('Џ®ЈаҐи­®бвм=');
   read(Eps);
   clrscr;

   textcolor(yellow);
   writeln('”®а¬г«  ва ЇҐжЁЁ:');
   writeln;
  repeat
   x:=x-bank;
   step(x,y,1,z);
   hs:=Eps1*(b-a)/Eps;
   bank:=h-hs;
   q:=trap(x,z-bank);
   i:=i+1;
   x:=z;y:=b;
   H1:=H1+hs;
   Q1:=Q1+q;
   Eps2:=Eps2+Eps1;textcolor(lightgreen);
   writeln('i= ' ,i,' h= ',hs:3,' Q= ',q:3,' e=',Eps1:3);
  until x=y  ;
   step(b-bank,b,1,z);
   H1:=H1+h;
   Q1:=Q1+q;
   Eps2:=Eps2+Eps1;textcolor(lightred);
   writeln('Ў ­Є: h=',h:2:2,' Q=',q:2:2,' e=',eps1);
   writeln('б㬬Ёа®ў ­ЁҐ: h=',h1:2:2,' Q=',q1,' e=',eps2);

   readln;readln;
   x:=a;y:=b;
   j:=0;
   H2:=0;Q2:=0;
   hch:=0;
   Eps3:=0;bank:=0;textcolor(yellow);
   writeln('”®а¬г«  ѓ гбб :');
   writeln;
  repeat
   x:=x-bank;
   step(x,y,2,z);
   j:=j+1;
   hch:=Eps1*(b-a)/Eps;
   bank:=h-hch;
   q:=gauss(x,z-bank);
   x:=z;
   y:=b;
   H2:=H2+hch;Q2:=Q2+q;
   textcolor(lightgreen);
   Eps3:=Eps3+Eps1;
   writeln('i=' ,j:2,' h= ',hch:2:2,' Q=',q:2:2,' e=',Eps1);
  until x=y  ;
   step(b-bank,b,2,z);
   j:=j+1;
   H2:=H2+h;
   Q2:=Q2+q;
   Eps3:=Eps3+Eps1;
   textcolor(lightred);
   writeln('Ў ­Є: h=',h:10:5,' Q=',q:10:5,' e=',eps1);
   writeln('б㬬Ёа®ў ­ЁҐ: h=',h2:10:5,' Q=',q2,' e=',eps3);

   readln;
 end.