Добавил:
bagiwow
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:выч методы не разобраны! / LAB2-И~1 / 2ЛАБА~1 / МОЕ / L2
.PAS {$N+,E+}
{Џа®Ја ¬¬ ўлзЁб«пҐв ЁвҐЈа «« ¬Ґв®¤®¬ —ҐЎл襢 ў ¤ ЇвЁў®¬ «Ј®аЁв¬Ґ}
Uses Crt;
Const b=1;a=0;e=0.001;t=0.7071;
function f(x:extended):extended;
begin
f:=(x*x)/((1-x*x*x*x*x)*(1-x*x*x*x*x)*(1-x*x*x*x*x)+2)
end;
{Ћб®ў п Їа®Ја ¬¬ }
Var
xl,xm,xr,x,Theb1,Theb2 :extended;{‡ 票Ґ x Єа пе ®в१Є Ё Ї®бҐаҐ¤ЁҐ}
S,P,Q,h,ee:extended;{‡ 票Ґ ЁвҐЈа « }
ei,hi,Qi:array[1..1000]of extended;
xl1,xm1,xr1,x1 :extended;{‡ 票Ґ x Єа пе ®в१Є Ё Ї®бҐаҐ¤ЁҐ}
S1,P1,Q1,h1,ee1:extended;{‡ 票Ґ ЁвҐЈа « }
ei1,hi1,Qi1:array[1..1000]of extended;
i,j:byte;
function Theb(a1,b1:extended):extended;
var x1,x2,x3:extended;
begin
x1:=(b1+a1)/2-(b1-a1)*t/2;
x2:=(b1+a1)/2+(b1-a1)*t/2;
x3:=(b1+a1)/2;
Theb2:=(h*0.5*((b-a)/2)*(f(x1)+f(x2)+f(x3)));
Theb1:=Theb1+Theb2;
Theb:=Theb2;
end;
BEGIN
Theb1:=0;Theb2:=0;
ClrScr;
Writeln('‚лзЁб«ҐЁҐ ЁвҐЈа « ЎҐ§ Ў ЄЁа®ў Ёп');
WriteLn('j':5, '|', 'h(i)':17,'Q(i)':21,'e(i)':17);
WriteLn('-----------------------------------------------------------------------------');
xl:=a;xr:=b;i:=0;
repeat{–ЁЄ« Ї®ЁбЄ ®в१Є ,㤮ў«Ґвў®апо饣® Ґа ўҐбвўг}
repeat h:=xr-xl;{ђ §¬Ґа ®в१Є }
P:=Theb(xl, xr); {1*h}
xm:=(xl+xr)/2; {‘ҐаҐ¤Ё ®в१Є }
Q:=Theb(xl, xm)+Theb(xm, xr); {2*h}
x:=xr; {е-а ўҐ a}
xr:=xm; {xr-а ўҐ бҐаҐ¤ЁҐ ®в१Є }
ee:=Abs(P-Q)/63 {ўлз-Ґ¬ Ї®ЈаҐи®бвм н⮬ гз бвЄҐ}
until ee<={h*}e/(b-a); {…б«Ё Ї®ЈаҐи®бвм ¬ҐмиҐ}
i:=i+1;
Qi[i]:=Q; {‡ 票Ґ ЁвҐЈа « ®в१ЄҐ}
hi[i]:=h; {а §¬Ґа нв®Ј® ®в१Є }
ei[i]:=ee; {Ї®ЈаҐи®бвм ®в१ЄҐ}
xl:=xm; {ЇаЁбў Ґў Ґ¬ з «г ®в१Є ®ў®Ґ § 票Ґ}
xr:=b {ЇаЁбў ЁўЁҐ¬ Є®жг ®в१Є § 票Ґ b}
until xl>=b-h; {ЇҐаҐе®¤ Є б«Ґ¤го饬㠮в१Єг}
h:=0;S:=0;ee:=0;
for j:=1to i do begin if (j div 5)=0 then WriteLn(j:5, '|', hi[j], '|', Qi[j], '|', ei[j]);
S:=S+Qi[j];
h:=h+hi[j];
ee:=ee+ei[j]
end;
WriteLn('-----------------------------------------------------------------------------');
WriteLn(' ':5, ' ', h:5:5, ' | ', S:7:7, ' | ' , ee:5:5);
WriteLn;
WriteLn('‡ ¤ Ї®ЈаҐи®бвм : ',e:5:5);
WriteLn('‡ 票Ґ ЁвҐЈа « : ',S:7:7);readln;
xl1:=a;xr1:=b;i:=0;
repeat{–ЁЄ« Ї®ЁбЄ ®в१Є ,㤮ў«Ґвў®апо饣® Ґа ўҐбвўг}
repeat h1:=xr1-xl1;{ђ §¬Ґа ®в१Є }
P1:=Theb(xl1, xr1);
xm1:=(xl1+xr1)/2;
Q1:=Theb(xl1, xm1)+Theb(xm1, xr1);
x1:=xr1;
xr1:=xm1;
ee1:=Abs(P1-Q1)/63
until ee1<={h1*}e/(b-a);
i:=i+1;
Qi1[i]:=Q1;
hi1[i]:=h1;
ei1[i]:=ee1;
xl1:=x1;
xr1:=b
until xl1>=b-h1;
h1:=0;S1:=0;ee1:=0;
Writeln('‚лзЁб«ҐЁҐ ЁвҐЈа « c Ў ЄЁа®ў ЁҐ¬');
WriteLn('j':5, '|', 'h(i)':17,'Q(i)':21,'e(i)':17);
WriteLn('-------------------------------------------------------------------------------');
for j:=1to i do begin
S1:=S1+Qi1[j];
h1:=h1+hi1[j];
ei1[j+1]:=ei1[j+1]+ei1[j];
ee1:=ee1+ei1[j];
WriteLn(j:5, '|', hi1[j], '|', Qi1[j], '|', ei1[j]*10);
end;
WriteLn('-------------------------------------------------------------------------------');
WriteLn(' ':5, ' ', h1:5:5, ' | ', S1:5:5, ' | ' , ee1*10:5:5);
WriteLn('‡ ¤ Ї®ЈаҐи®бвм : ',e:5:5);
WriteLn('‡ 票Ґ ЁвҐЈа « : ',S1:5:5);readln
END.
{Џа®Ја ¬¬ ўлзЁб«пҐв ЁвҐЈа «« ¬Ґв®¤®¬ —ҐЎл襢 ў ¤ ЇвЁў®¬ «Ј®аЁв¬Ґ}
Uses Crt;
Const b=1;a=0;e=0.001;t=0.7071;
function f(x:extended):extended;
begin
f:=(x*x)/((1-x*x*x*x*x)*(1-x*x*x*x*x)*(1-x*x*x*x*x)+2)
end;
{Ћб®ў п Їа®Ја ¬¬ }
Var
xl,xm,xr,x,Theb1,Theb2 :extended;{‡ 票Ґ x Єа пе ®в१Є Ё Ї®бҐаҐ¤ЁҐ}
S,P,Q,h,ee:extended;{‡ 票Ґ ЁвҐЈа « }
ei,hi,Qi:array[1..1000]of extended;
xl1,xm1,xr1,x1 :extended;{‡ 票Ґ x Єа пе ®в१Є Ё Ї®бҐаҐ¤ЁҐ}
S1,P1,Q1,h1,ee1:extended;{‡ 票Ґ ЁвҐЈа « }
ei1,hi1,Qi1:array[1..1000]of extended;
i,j:byte;
function Theb(a1,b1:extended):extended;
var x1,x2,x3:extended;
begin
x1:=(b1+a1)/2-(b1-a1)*t/2;
x2:=(b1+a1)/2+(b1-a1)*t/2;
x3:=(b1+a1)/2;
Theb2:=(h*0.5*((b-a)/2)*(f(x1)+f(x2)+f(x3)));
Theb1:=Theb1+Theb2;
Theb:=Theb2;
end;
BEGIN
Theb1:=0;Theb2:=0;
ClrScr;
Writeln('‚лзЁб«ҐЁҐ ЁвҐЈа « ЎҐ§ Ў ЄЁа®ў Ёп');
WriteLn('j':5, '|', 'h(i)':17,'Q(i)':21,'e(i)':17);
WriteLn('-----------------------------------------------------------------------------');
xl:=a;xr:=b;i:=0;
repeat{–ЁЄ« Ї®ЁбЄ ®в१Є ,㤮ў«Ґвў®апо饣® Ґа ўҐбвўг}
repeat h:=xr-xl;{ђ §¬Ґа ®в१Є }
P:=Theb(xl, xr); {1*h}
xm:=(xl+xr)/2; {‘ҐаҐ¤Ё ®в१Є }
Q:=Theb(xl, xm)+Theb(xm, xr); {2*h}
x:=xr; {е-а ўҐ a}
xr:=xm; {xr-а ўҐ бҐаҐ¤ЁҐ ®в१Є }
ee:=Abs(P-Q)/63 {ўлз-Ґ¬ Ї®ЈаҐи®бвм н⮬ гз бвЄҐ}
until ee<={h*}e/(b-a); {…б«Ё Ї®ЈаҐи®бвм ¬ҐмиҐ}
i:=i+1;
Qi[i]:=Q; {‡ 票Ґ ЁвҐЈа « ®в१ЄҐ}
hi[i]:=h; {а §¬Ґа нв®Ј® ®в१Є }
ei[i]:=ee; {Ї®ЈаҐи®бвм ®в१ЄҐ}
xl:=xm; {ЇаЁбў Ґў Ґ¬ з «г ®в१Є ®ў®Ґ § 票Ґ}
xr:=b {ЇаЁбў ЁўЁҐ¬ Є®жг ®в१Є § 票Ґ b}
until xl>=b-h; {ЇҐаҐе®¤ Є б«Ґ¤го饬㠮в१Єг}
h:=0;S:=0;ee:=0;
for j:=1to i do begin if (j div 5)=0 then WriteLn(j:5, '|', hi[j], '|', Qi[j], '|', ei[j]);
S:=S+Qi[j];
h:=h+hi[j];
ee:=ee+ei[j]
end;
WriteLn('-----------------------------------------------------------------------------');
WriteLn(' ':5, ' ', h:5:5, ' | ', S:7:7, ' | ' , ee:5:5);
WriteLn;
WriteLn('‡ ¤ Ї®ЈаҐи®бвм : ',e:5:5);
WriteLn('‡ 票Ґ ЁвҐЈа « : ',S:7:7);readln;
xl1:=a;xr1:=b;i:=0;
repeat{–ЁЄ« Ї®ЁбЄ ®в१Є ,㤮ў«Ґвў®апо饣® Ґа ўҐбвўг}
repeat h1:=xr1-xl1;{ђ §¬Ґа ®в१Є }
P1:=Theb(xl1, xr1);
xm1:=(xl1+xr1)/2;
Q1:=Theb(xl1, xm1)+Theb(xm1, xr1);
x1:=xr1;
xr1:=xm1;
ee1:=Abs(P1-Q1)/63
until ee1<={h1*}e/(b-a);
i:=i+1;
Qi1[i]:=Q1;
hi1[i]:=h1;
ei1[i]:=ee1;
xl1:=x1;
xr1:=b
until xl1>=b-h1;
h1:=0;S1:=0;ee1:=0;
Writeln('‚лзЁб«ҐЁҐ ЁвҐЈа « c Ў ЄЁа®ў ЁҐ¬');
WriteLn('j':5, '|', 'h(i)':17,'Q(i)':21,'e(i)':17);
WriteLn('-------------------------------------------------------------------------------');
for j:=1to i do begin
S1:=S1+Qi1[j];
h1:=h1+hi1[j];
ei1[j+1]:=ei1[j+1]+ei1[j];
ee1:=ee1+ei1[j];
WriteLn(j:5, '|', hi1[j], '|', Qi1[j], '|', ei1[j]*10);
end;
WriteLn('-------------------------------------------------------------------------------');
WriteLn(' ':5, ' ', h1:5:5, ' | ', S1:5:5, ' | ' , ee1*10:5:5);
WriteLn('‡ ¤ Ї®ЈаҐи®бвм : ',e:5:5);
WriteLn('‡ 票Ґ ЁвҐЈа « : ',S1:5:5);readln
END.