Priblizhennye_vychislenia_kratnykh_integralov / лаба 5

«Уфимский государственный нефтяной технический университет»

Кафедра: Прикладная математика и механика

Лабораторная работа №5

по информатике

«Приближенные вычисления кратных интегралов»

Вариант№13

Выполнил: ст. гр. БДС 11-01 Закиров Р.Х.

Проверил: Виноградов Д.А.

Уфа 2012

Здание:

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.

Результат:

Границы области интегрирования:

а=0

b=1

c=0

d=1

Число разбиений по осям:

n=100

m=00

Значение интеграла методом ячеек

v=0.42039