1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
| program Zadanie9;
var x,s1,s2,s,s3,s4,i,n,n1,n2,n3,z,FC1,FC2:real;
const
b0 = 14; k = -0.47;
x1 = -4; y1 = -9; r = 28;
x2 = 7.2; y2 = 8; a = 18; b = 27;
function FCline(xline: real): real;
begin
FCLine := k * xline + b0;
end;
function FCcircle(xcircle: real): real;
begin
FCcircle := sqrt(abs(sqr(xcircle - x1) - sqr(r)));
end;
function FCellips(xellips: real): real;
begin
FCellips := sqrt(abs(sqr(b) - sqr(b) /sqr(a) * sqr(xellips - x2) ));
end;
function FCellips1(xellips1: real): real;
begin
FCellips1 := sqrt(abs(sqr(b) - sqr(b) /sqr(a) * sqr(xellips1-x2 )));
end;
procedure calculate (n:real);
begin
x:=-9.4;
begin
repeat
FC1:=FCLine(x);
s1:=s1+FC1*35/n;
x:=x+35/n;
until x>=0;
x:=-10.79;
repeat
FC1:=-FCellips1(x)+y2;
FC2:=Fcellips(x)+y2;
x:=x+35/n;
s1:=s1+(FC2-FC1)*35/n;
until x>=-10;
x:=-10;
repeat
FC1:= FCellips(x)+y2;
s1:=s1+FC1*35/n;
x:=x+35/n;
until x>=-9.4;
x:=0;
repeat
FC1:=-FCcircle(x)-y1;
FC2:=FCellips(x)-y2;
s2:=s2+(FC2-FC1)*35/n;
x:=x+35/n;
until x>=20.66;
x:=20.65;
repeat
FC1:=FCline(x);
FC2:=FCellips(x)+y2;
S2:=S2+(FC2-FC1)*35/n;
x:=x+35/n;
until x>=24.8;
x:=24.81;
repeat
FC1:=FCellips(x)+y2;
FC2:=-Fcellips(x)+y2;
x:=x+35/n;
s2:=s2+(FC1-FC2)*35/n;
until x>=25.2; s1:=s1-57; s2:=s2-124; s3:=s1+n1*0.017; s4:=s2+n1*0.017;
writeln('| ',n,' | ',s1:0:3,' | ',s3:0:3,' | | ',s2:0:3,' | ',s4:0:3,' |');
s1:=0;
s2:=0;
s3:=0;
x:=0;
end;
end;
begin
repeat
s1:=0;
s2:=0;
writeln('Введите чило интеракций');
read(n1,n2,n3);
writeln('| N | Площадь плоскости А |Площадь плоскости А| | Площадь плоскости В |Площадь плоскости В|');
writeln('| |методом прямоугольника| методом трапеции | |методом прямоугольника| методом трапеции |');
Calculate(n1);
Calculate(n2);
Calculate(n3);
writeln('end?');
read(z);
until z=1;
end. |