Скачиваний:
312
Добавлен:
04.03.2014
Размер:
1.71 Кб
Скачать
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