Скачиваний:
329
Добавлен:
04.03.2014
Размер:
2 Кб
Скачать
uses crt;
var
  c,d,h : real;
  a,b,x : real;
  i,j   : integer;
  r,r1,r2 : real;
  n       : word;
  e,iz    : real;
{*********}
function f(x:real):real;
begin
f:=cos(3*x)*sin(5*x);
end;
{*********}
function t1(n:integer):real;
var q : real;
begin
q:=0;
for j:=1 to n do begin
                 h:=(d-c)/(2*n);
                 a:=c+(j-1)*h*2;
                 b:=c+(j)*h*2;
                 x:=(a+b)/2;
                 q:=q+2*h*f(x);
                 end;
t1:=q;
end;
{*********}
function t2(n:integer):real;
var q : real;
begin
q:=0;
for j:=1 to n do begin
                 h:=(d-c)/(2*n);
                 a:=c+(j-1)*h*2;
                 b:=c+j*h*2;
                 x:=(a+b)/2;
                 {q:=q+h*(f(x-h/sqrt(3))+f(x+h/sqrt(3)));}
                 q:=q+h/3*(f(a)+4*f(x)+f(b));
                 end;
t2:=q;
end;
{*********}
{for i:=1 to n do q:=q+h/9*(5*f(x[i]-h*sqrt(3/5))+8*f(x[i])+5*f(x[i]+h*sqrt(3/5)));
for i:=1 to n do q:=q+2*h*f(a[i]);
for i:=1 to n do q:=q+2*h*f(b[i]);
for i:=1 to n do q:=q+h*f(a[i]+b[i]);
for i:=1 to n do q:=q+h/3*(f(a[i])+4*f(x[i])+f(b[i]));}
{********}
begin
clrscr;
Writeln('True F(b)-F(a)=0.4');
c:=1;
d:=4;
n:=0;
e:=0.000001;
repeat
inc(n);
r1:=(t1(2*n)-t1(n))/(3);   {s=0 4,5 s=1 1,6 s=3 2,7 s=5 3}
until abs(r1)<=e;
writeln('N ==> ',n);
writeln('I1 ==> ',t1(n));
writeln('r ==> ',r1);
writeln('I1* ==> ',t1(2*n)+r1);
writeln('I1-I1* ==> ',t1(n)-t1(2*n)+r1);
writeln;
n:=0;
repeat
inc(n);
r2:=(t2(2*n)-t2(n))/(15);
until abs(r2)<=e;
writeln('N ==> ',n);
writeln('I1 ==> ',t2(2*n));
writeln('r ==> ',r2);
writeln('I1* ==> ',t2(2*n)+r2);
writeln('I1-I1* ==> ',t2(n)-t2(2*n)+r2);
readln;
end.
Соседние файлы в папке Alex