 |
Функция :: Интеграл (теория)
Функция :: Интеграл
uses
crt,mat;
type
float=extended;
const
xn=0.001;
xk=0.2;
e=0.001;
var
x,h,f_left,f_sred,f_prav:float;
I1,I2:float;
function f (x:float):float;
begin
f:=((sin(16*x))/(3*x))+4*sqr(arctan(1))+kor(log(1,2),2);
end;
Begin
clrscr;
I1:=0;
I2:=0;
h:=(xk-xn)/10;
x:=xn;
while x
begin
I1:=I1+f(x)*h;
x:=x+h
end;
Repeat
I2:=I1;
I1:=0;
h:=h/2;
x:=xn;
while x
begin
I1:=I1+f(x)*h;
x:=x+h
end;
Until abs(I2-I1)
f_left:=I1;
I1:=0;
I2:=0;
h:=(xk-xn)/5;
x:=xn;
while x
begin
I1:=I1+f(x+h/2)*h;
x:=x+h
end;
Repeat
I2:=I1;
I1:=0;
h:=h/2;
x:=xn;
while x
begin
I1:=I1+f(x+h/2)*h;
x:=x+h
end;
Until abs(I2-I1)
f_sred:=I1;
I1:=0;
I2:=0;
h:=(xk-xn)/5;
x:=xn;
while x
begin
I1:=I1+f(x+h)*h;
x:=x+h
end;
repeat
I2:=I1;
I1:=0;
h:=h/2;
x:=xn;
while x
begin
I1:=I1+f(x+h)*h;
x:=x+h;
end;
until abs(I2-I1)
f_prav:=I1;
writeln('integral f(x) - levie treugolniki ot ',xn:0:2,' do ',xk:0:2,' = ',f_left:0:2);
writeln('integral f(x) - srednie treugolniki ot ',xn:0:2,' do ',xk:0:2,' = ',f_sred:0:2);
writeln('integral f(x) - pravie treugolniki ot ',xn:0:2,' do ',xk:0:2,' = ',f_prav:0:2);
readln;
end.
|
