A = (1.5, 1.02)
a_1 = 4
B = Point(Circle[A, 4])
a = Segment(A, B)
poly1, b, c, d, e, C, D = Polygon(A, B, 4)
f = OrthogonalLine(A, poly1)
g = Circle(A, 6, a)
S = Intersect(f, g, 2)
h, faceABE, faceBCE, faceCDE, faceADE, arêteAE, arêteBE, arêteCE, arêteDE = Pyramid(poly1, S)
M = Point(arêteAE)
i = Plane(M, poly1)
N = Intersect(arêteBE, i)
O = Intersect(arêteCE, i)
P = Intersect(arêteDE, i)
poly2, m, n, o, p = Polygon(M, N, O, P)
Aire: (Distance[S, M], poly2)
r = Segment((Distance[S, M], 0), Aire)
s = Segment((0, poly2), Aire)
F: (Distance[S, M], 0)
G: (0, poly2)
Dist = Distance(S, M)
texte1: "SM=" + (LaTeX[Dist]) + ""
texte2: "Aire = " + (LaTeX[poly2]) + ""
rapport1: Distance[S, M] / 6
rapport2: poly2 / poly1
texte4: "Aire MNOP / Aire ABCD = " + (LaTeX[rapport2]) + ""
texte3: "SM/SA= " + (LaTeX[rapport1]) + ""
"SM/6= " + (LaTeX[rapport1]) + "": "SM/SA= " + (LaTeX[rapport1]) + ""