Добавил:
bagiwow
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:выч методы не разобраны! / LAB2-И~1 / 2ЛАБА~1 / INTEGRAL / INTEGRAL
.PAS Program Integral;
uses crt;
var x1, x2, Int, h, a, b, f4, e:real;
function f(_x_:real):real;
begin
f:=sqr(_x_)/(sqr(1-sqr(sqr(_x_))*_x_)*(1-sqr(sqr(_x_))*_x_)+2);
end;
begin
clrscr;
a:=0; b:=1; h:=(b-a)/4;
x1:=0.5*(1-1/sqrt(3));
x2:=0.5*(1+1/sqrt(3));
Int:=0.5*(f(x1)+f(x2));
f4:=(f(b)-4*f(b-h)+6*f(b-2*h)-4*f(b-3*h)+f(a))/(sqr(sqr(h)));
e:=1/11520*sqr(sqr(b-a))*(b-a)*f4;
writeln('‡ 票Ґ ЁвҐЈа « :',Int:0:5);
writeln('Џ®ЈаҐи®бвм:',e:0:5);
writeln('‡ ¤ п Ї®ЈаҐи®бвм: 0.001');
readln;
end.
uses crt;
var x1, x2, Int, h, a, b, f4, e:real;
function f(_x_:real):real;
begin
f:=sqr(_x_)/(sqr(1-sqr(sqr(_x_))*_x_)*(1-sqr(sqr(_x_))*_x_)+2);
end;
begin
clrscr;
a:=0; b:=1; h:=(b-a)/4;
x1:=0.5*(1-1/sqrt(3));
x2:=0.5*(1+1/sqrt(3));
Int:=0.5*(f(x1)+f(x2));
f4:=(f(b)-4*f(b-h)+6*f(b-2*h)-4*f(b-3*h)+f(a))/(sqr(sqr(h)));
e:=1/11520*sqr(sqr(b-a))*(b-a)*f4;
writeln('‡ 票Ґ ЁвҐЈа « :',Int:0:5);
writeln('Џ®ЈаҐи®бвм:',e:0:5);
writeln('‡ ¤ п Ї®ЈаҐи®бвм: 0.001');
readln;
end.
Соседние файлы в папке INTEGRAL