Скачиваний:
288
Добавлен:
04.03.2014
Размер:
2.98 Кб
Скачать
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