Добавил:
bagiwow
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:выч методы не разобраны! / LAB2-И~1 / 2ЛАБА~1 / МОЕ / INT___
.PAS Program Integral;
uses crt;
const eps=0.001;
x0=0;
xn=1;
t=0.7947;
var x1,x2,X3,Int,h,a,b:real;
function f(x:real):real;
begin
f:=(x*x)*sin(4*ln(x))/ln(x)
end;
function epsi(a,b:real):real;
var h, f4,x9:real;
begin
h:=(a-b)/4;
x9:=(b+a)/2;
f4:=((2*x9*sin(4*ln(x9))/ln(x9))+(4*x9*cos(4*ln(x9))/ln(x9))+(x9*sin(4*ln(x9))/ln(x9)*ln(x9)));
epsi:=(1/2721600)*exp(5*ln(b-a))*f4;
end;
function Theb(a1,b1:extended):extended;
begin
x1:=(b1+a1)/2-(b1-a1)/2*t;
x2:=(b1+a1)/2+(b1-a1)/2*t;
x3:=(b1+a1)/2;
Int:=Int+(0.5*(f(x1)+f(x2)+f(x3)));
end;
begin
clrscr;
a:=x0; b:=xn;
Int:=0;
if abs(epsi(a,b))<=eps*(b-a)/(xn-x0) then
begin
write('[a:=',a:0:3,',b:=',b:0:3,'] ');
Theb(a,b);
write('Integral:=',Int:0:6,' ');
write('Epsi:=',epsi(a,b):0:6,' ');
write('Eps:=',eps*(b-a)/(xn-x0):0:6);
writeln;
a:=b; b:=xn;
end
else b:=(a+b)/2;
writeln('Integral:=',Int:0:6);
readln;
end.
uses crt;
const eps=0.001;
x0=0;
xn=1;
t=0.7947;
var x1,x2,X3,Int,h,a,b:real;
function f(x:real):real;
begin
f:=(x*x)*sin(4*ln(x))/ln(x)
end;
function epsi(a,b:real):real;
var h, f4,x9:real;
begin
h:=(a-b)/4;
x9:=(b+a)/2;
f4:=((2*x9*sin(4*ln(x9))/ln(x9))+(4*x9*cos(4*ln(x9))/ln(x9))+(x9*sin(4*ln(x9))/ln(x9)*ln(x9)));
epsi:=(1/2721600)*exp(5*ln(b-a))*f4;
end;
function Theb(a1,b1:extended):extended;
begin
x1:=(b1+a1)/2-(b1-a1)/2*t;
x2:=(b1+a1)/2+(b1-a1)/2*t;
x3:=(b1+a1)/2;
Int:=Int+(0.5*(f(x1)+f(x2)+f(x3)));
end;
begin
clrscr;
a:=x0; b:=xn;
Int:=0;
if abs(epsi(a,b))<=eps*(b-a)/(xn-x0) then
begin
write('[a:=',a:0:3,',b:=',b:0:3,'] ');
Theb(a,b);
write('Integral:=',Int:0:6,' ');
write('Epsi:=',epsi(a,b):0:6,' ');
write('Eps:=',eps*(b-a)/(xn-x0):0:6);
writeln;
a:=b; b:=xn;
end
else b:=(a+b)/2;
writeln('Integral:=',Int:0:6);
readln;
end.