Priblizhennye_vychislenia_kratnykh_integralov / LABA52 ячейки

program cells;
uses printer, crt;
var a,b,c,d,hx,hy,xc,yc,v:real;
x,y:array[0..100] of real;
i,j,m,n:integer;
function f(x,y:real):real;
begin
f:=x*x+sqr(y-1);
end;
begin
clrscr;
write('a=');
readln(a);
write('b=');
readln(b);
write('c=');
readln(c);
write('d=');
readln(d);
write('n=');
readln(n);
write('m=');
readln(m);
hx:=(b-a)/n;
hy:=(d-c)/m;
x[0]:=a;
y[0]:=c;
v:=0;
for i:=1 to n do
begin
x[i]:=a+i*hx;
xc:=(x[i]+x[i-1])/2;
for j:=1 to m do
begin
y[i]:=c+j*hy;
yc:=(y[i]+y[i-1])/2;
v:=v+f(xc,yc);
end;
end;
v:=v*hx*hy;
writeln('v=', v:2:5);
end.