A = (-2, 2)
B = (-2, -2)
C = (2, -2)
D = (2, 2)
poly1, b, c, d, a = Polygon(B, C, D, A)
S: (x(A), y(A), 6)
e, faceBCS, faceCDS, faceADS, faceABS, arêteBS, arêteCS, arêteDS, arêteAS = Pyramid(poly1, S)
α = Angle(B, A, S)
β = Angle(D, A, S)
γ = Angle(D, A, B, space)
M = Point(arêteAS)
f = Plane(M, poly1)
poly2, poly3, E, F, G, H, J, K, L, N, g, h, i, j, k, l, m, n = IntersectPath(f, e)
I = (-3.2, 4.88)
M': (Distance[S, M], poly2)