d= sqrt(sqr(x2-x1)+sqr(y2-y1)); { расстояние между центрами }

Program area

Implicit none

Real: r1,r2,x1,x2,y1,y2

Real:d,seg1,seg2,s

Real:alfa1,alfa2

Real: s1,s2,sec1,sec2

Real: x,h1,h2,b1,b2

Print*,’ okruzhnost 1’

Read*, r1,x1,y1

Print*,’okruzhnost 2’

Read*,r2,x2,y2

d= sqrt(sqr(x2-x1)+sqr(y2-y1)); { расстояние между центрами }

if R1+R2<=d then S=0 { если круги не пересекаются...}

else if d<=R1-R2 then { если один круг полностью в другом }

if R1>R2 then S= pi*sqr(R2)

else S= pi*sqr(R1)

else

x:= R1+R2 - d

if R1<R2 then

R1 = R1+R2

R2 = R1-R2

R1 = R1-R2

end

h1 = R1 - (x/2)*(R2/R1)

h2 = R2 - (x/2)*(R1/R2)

b1 = 2*sqrt(sqr(R1)-sqr(h1))

b2 = 2*sqrt(sqr(R2)-sqr(h2))

S1 = b1*h1/2

S2 = b2*h2/2

alpha1 = 2*arctan(b1/(2*h1))

alpha1 = alpha1*180/pi

alpha2 = 2*arctan(b2/(2*h2))

alpha2 = alpha2*180/pi

Sec1 = pi*sqr(R1)*alpha1/360

Sec2 = pi*sqr(R2)*alpha2/360

if alpha1>180 then S1 = -S1

if alpha2>180 then S2 = -S2

Seg1 = Sec1 - S1;

Seg2 = Sec2 - S2;

S= Seg1 + Seg2;

Print*,S

end

Дата добавления: 2015-11-04