Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
6
Добавлен:
06.03.2016
Размер:
746 б
Скачать
program posint;
uses printer, crt;
var a,b,c,d,x,y,hx,hy,v:real;
s:array[0..100] of real; i,j,n,m:integer;
function f(x,y:real):real;
begin
f:=1-sqrt(x*x+y*y);
end;
function fi(x:real):real;
begin
fi:=0
end;
function psi(x:real):real;
begin
psi:=ln(x);
end;
begin
clrscr;
write('a='); readln(a);
write('b='); readln(b);
write('hx='); readln(hx); write('hy='); readln(hy);
n:=trunc((b-a)/hx);
for i:=0 to n do
begin
x:=a+i*hx;
c:=fi(x); d:=psi(x);
m:=trunc((d-c)/hy);
s[i]:=(f(x,c)+f(x,d))/2;
for j:=1 to m-1 do
begin
y:=c+j*hy;
s[i]:=s[i]+f(x,y);
end;
s[i]:=s[i]*hy;
end;
v:=(s[0]+s[n])/2;
for i:=1 to n-1 do v:=v+s[i];
v:=v*hx;
writeln('v=', v:2:5);
end.
Соседние файлы в папке Priblizhennye_vychislenia_kratnykh_integralov