Скачиваний:
329
Добавлен:
04.03.2014
Размер:
3 Кб
Скачать
program lab1;
uses crt;
label as,as1;
var
E,H,i1,i1tek,i2,i2tek,i3,i3tek,r:real;
J,n,m,k,c:integer;ch:char;a,b:real; v:char;
function  F(X:real):real;
begin
F:=x*x*x;end;
procedure PR(H:real;var i:real);
begin
I:=0;
for J:=1 to n do
I:=I+f(A+(J-0.5)*h);
I:=I*h;end;
procedure TR(h:real;var i:real);
begin
i:=(f(A)+f(B))/2;
for J:=1 to n-1 do
I:=i+F(a+j*h);
i:=i*h;end;
procedure SM(h:real;var i:real);
var
m:integer;
begin
I:=f(a)+f(b);
M:=2;
for J:=1 to (n-1) do
begin
 if M=4 then m:=2 else M:=4;
 i:=i+m*f(a+h*j);end;
 i:=h/3*i;end;
BEGIN
{goto as;
as:}
clrscr;
        WRITELN('1)ЊҐв®¤ Їаאַ㣮«м­ЁЄ®ў');
        writeln('2)ЊҐв®¤ ва ЇҐжЁ©');
        writeln('3)ЊҐв®¤ ‘Ё¬Їб®­ ');
        write('ўўҐ¤ЁвҐ ­®¬Ґа ¬Ґв®¤ ');
k:=-5;
repeat
case c of
1: begin


       writeln(' ‚ўҐ¤ЁвҐ ўҐае­Ё© ЇаҐ¤Ґ« Ё­вҐЈаЁа®ў ­Ёп ');
         write('a=');
       readln(A);

       writeln(' ‚ўҐ¤ЁвҐ ­Ё¦­Ё© ЇаҐ¤Ґ« Ё­вҐЈаЁа®ў ­Ёп ');
       write('b=');
       readln(B); writeln('ўўҐ¤ЁвҐ зЁб«® а §ЎЁҐ­Ё© ');
       write('m=');
       readln(m);

         write('‡Ђ„Ђ‰ ’Ћ—ЌЋ‘’њ E=');REadln(e);

 writeln('¬Ґв®¤ Їаאַ㣮«м­ЁЄ®ў ');
 writeln('       Ip               I              Ip-I  ');
 H:=sqrt(e);
 n:=trunc((b-a)/H)+1;
 h:=(b-a)/n;
 repeat
 for k:=1 to m do
 begin
 h:=h/2;n:=n*2;
 pr(h,i1tek);
 r:=abs(i1tek-i1);
 writeln(i1:16,i1tek:16,r:16);
i1:=i1tek;end;
 v:=readkey until v<>#13; end;
 2:begin


       writeln(' ‚ўҐ¤ЁвҐ ўҐае­Ё© ЇаҐ¤Ґ« Ё­вҐЈаЁа®ў ­Ёп ');
         write('a=');
       readln(A);

       writeln(' ‚ўҐ¤ЁвҐ ­Ё¦­Ё© ЇаҐ¤Ґ« Ё­вҐЈаЁа®ў ­Ёп ');
       write('b=');
       readln(B); writeln('ўўҐ¤ЁвҐ зЁб«® а §ЎЁҐ­Ё© ');
       readln(m);
         write('‡Ђ„Ђ‰ ’Ћ—ЌЋ‘’њ E=');REadln(e);
 h:=sqrt(e);
 n:=trunc((b-a)/h)+1;
 h:=(b-a)/n;
 writeln(' ¬Ґв®¤ ва ЇҐжЁ©');
 writeln('      Ipr          I         I-Ipr');
 repeat for k:=1 to m do
 begin
 h:=h/2;n:=n*2;tr(h,i2tek);R:=abs(i2-i2tek);
{ writeln('    Ipr       I         abs[I-Ipr]'); }
 writeln(i2:16,i2tek:16,r:16);

 i2:=i2tek;end;v:=readkey;

 until v<>#13;end;
 3:begin
 writeln(' ¬Ґв®¤ бЁ¬Їб®­   ');

  writeln(' ‚ўҐ¤ЁвҐ ўҐае­Ё© ЇаҐ¤Ґ« Ё­вҐЈаЁа®ў ­Ёп ');
         write('a=');
       readln(A);
       writeln(' ‚ўҐ¤ЁвҐ ­Ё¦­Ё© ЇаҐ¤Ґ« Ё­вҐЈаЁа®ў ­Ёп ');
       write('b=');
       readln(B); writeln('ўўҐ¤ЁвҐ зЁб«® а §ЎЁҐ­Ё© ');
       readln(m);
         write('‡Ђ„Ђ‰ ’Ћ—ЌЋ‘’њ E=');REadln(e);
 writeln('     Ip                 I              I-Ip    ');
 h:=sqrt(sqrt(e));
 n:=trunc((b-a)/h)+1;
 if 2*int(n/2)<>n then n:=n+1;
 H:=(b-a)/n;
 repeat for k:=1 to m do
   begin
   h:=h/2;n:=n*2;
   sm(h,i3tek);
   r:=abs(i3-i3tek);
  { writeln('     Ipr        I         abs[I-Ipr]'); }
   writeln(i3:16,'  ',i3tek:16,'  ',r:16);
 {  writeln(n); }
   i3:=i3tek;end;
   v:=readkey;
   until v<>#13;end;end;
   {WRITELN('‚ўҐ¤ЁвҐ ­®¬Ґа ¬Ґв®¤ ');}
   read(c);CLRSCR;
   until c=0;readln;
   END.
Соседние файлы в папке Alex