Отчет по лабораторной 4
.docxПермский государственный технический университет
Кафедра Информационных технологий и автоматизированных систем
Отчет по
Лабораторной работе №4
«Вычисление площади»
По дисциплине «Информатика»
Выполнил студент группы АУЦ-10
Рычин В.
Проверил ассистент кафедры ИТАС
Баранов А.А.
Пермь 2010
Цель: Разработать программу вычисления заданной площади, используя два различных метода вычислений: прямоугольников, трапеций.
1. Метод прямоугольников:
Заменим каждую полоску криволинейной трапеции прямоугольником с высотой, равной значению функции на левом конце промежутка. Площадь криволинейной фигуры заменится площадью сумм прямоугольников.
Формулы прямоугольников*: 1). h[Y(0)+Y(1)+ ... +Y(N-1)]
2). h[Y(1)+Y(2)+ ... +Y(N)]
2. Метод трапеций:
Если каждую полоску криволинейной трапеции заменить обычной трапецией, то площадь фигуры будет равна сумме площадей трапеций.
Формула трапеций*: h[(Y(0)+Y(N))/2+Y(1)+Y(2)+ ... +Y(N-1)]
Программа:
uses crt;
var x,y,y1,y2,y3,y4,y5,dx,dy,dy1,dy2,s,s1,s2,s3,
s4,s5,s6,s7,
s8,s9,s10,s11:real; i,a,b,c:integer;
begin;
clrscr;
writeln('BBeDu 4ucJLo Pa36uenuu');
readln(a,b,c);
begin
x:=0; dx:=28/a; s1:=0; s:=0;
for i:=0 to a-1 do
begin
if x<=15 then
y1:=abs(sqrt((1-sqr(x+5)/400)*1089)-2)
else y1:=abs(8.6-0.57*x);
y2:=abs(sqrt((900-sqr(x+2)))-5);
dy:=y2-y1;
s1:=dx*dy;
s:=s+s1;
x:=x+dx;
end;
x:=0; dx:=28/a; s2:=0; s3:=0;
for i:=0 to a-1 do
begin
if x<=15 then
y1:=abs(sqrt((1-sqr(x+5)/400)*1089)-2)
else y1:=abs(8.6-0.57*x);
y2:=abs(sqrt(900-sqr(x+2))-5);
y3:=abs(sqrt((1-sqr(dx+5)/400)*1089)-2);
y4:=abs(sqrt(900-sqr(dx+2))-5);
y5:=abs(8.6-0.57*dx);
dy1:=y2-y1;
if x<=15 then
dy2:=y4-y3
else dy2:=y4-y5;
dy:=(dy1+dy2)/2;
s2:=dx*dy;
s3:=s3+s2;
x:=x+dx;
end;
end;
begin
x:=0; dx:=28/b; s4:=0; s5:=0;
for i:=0 to b-1 do
begin
if x<=15 then
y1:=abs(sqrt((1-sqr(x+5)/400)*1089)-2)
else y1:=abs(8.6-0.57*x);
y2:=abs(sqrt(900-sqr(x+2))-5);
dy:=y2-y1;
s4:=dx*dy;
s5:=s5+s4;
x:=x+dx;
end;
x:=0; dx:=28/b; s6:=0; s7:=0;
for i:=0 to b-1 do
begin
if x<=15 then
y1:=abs(sqrt((1-sqr(x+5)/400)*1089)-2)
else y1:=abs(8.6-0.57*x);
y2:=abs(sqrt(900-sqr(x+2))-5);
y3:=abs(sqrt((1-sqr(dx+5)/400)*1089)-2);
y4:=abs(sqrt(900-sqr(dx+2))-5);
y5:=abs(8.6-0.57*dx);
dy1:=y2-y1;
if x<=15 then
dy2:=y4-y3
else dy2:=y4-y5;
dy:=(dy1+dy2)/2;
s6:=dx*dy;
s7:=s7+s6;
x:=x+dx;
end;
end;
begin
x:=0; dx:=28/c; s8:=0; s9:=0;
for i:=0 to c-1 do
begin
if x<=15 then
y1:=abs(sqrt((1-sqr(x+5)/400)*1089)-2)
else y1:=abs(8.6-0.57*x);
y2:=abs(sqrt(900-sqr(x+2))-5);
dy:=y2-y1;
s8:=dx*dy;
s9:=s9+s8;
x:=x+dx;
end;
x:=0; dx:=28/c; s10:=0; s11:=0;
for i:=0 to c-1 do
begin
if x<=15 then
y1:=abs(sqrt((1-sqr(x+5)/400)*1089)-2)
else y1:=abs(8.6-0.57*x);
y2:=abs(sqrt(900-sqr(x+2))-5);
y3:=abs(sqrt((1-sqr(dx+5)/400)*1089)-2);
y4:=abs(sqrt(900-sqr(dx+2))-5);
y5:=abs(8.6-0.57*dx);
dy1:=y2-y1;
if x<=15 then
dy2:=y4-y3
else dy2:=y4-y5;
dy:=(dy1+dy2)/2;
s10:=dx*dy;
s11:=s11+s10;
x:=x+dx;
end;
end;
writeln('N',' ','Metod pryamoug',' ','Metod trapecy');
writeln(a,' ',s:6:4,' ',s3:6:4);
writeln(b,' ',s5:6:4,' ',s7:6:4);
writeln(c,' ',s9:6:4,' ',s11:6:4);
end.