Скачиваний:
27
Добавлен:
09.12.2013
Размер:
1.14 Кб
Скачать
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.
Соседние файлы в папке МОЕ
  • #
    09.12.20133.88 Кб25INTM.BAK
  • #
    09.12.201320.34 Кб25INTM.EXE
  • #
    09.12.20133.88 Кб27INTM.PAS
  • #
    09.12.20131.16 Кб25INT___.BAK
  • #
    09.12.201318.8 Кб25INT___.EXE
  • #
    09.12.20131.14 Кб27INT___.PAS
  • #
    09.12.20133.9 Кб26L2.BAK
  • #
    09.12.20133.91 Кб26L2.PAS