A = (3, -1, 0)
B = (3, 3, 0)
C = (-1, 3, 0)
D = (-1, -1, 0)
v = (0, 0, 6)
E: A + v
F: B + v
G: C + v
H: D + v
e, I, J, K, FlächeABIE, FlächeBCJI, FlächeCDKJ, FlächeADKE, FlächeEIJK, KanteAE, KanteBI, KanteCJ, KanteDK, KanteEI, KanteIJ, KanteJK, KanteEK = Prism(Polygon[A, B, C, D], E)
L = Point(e)
A_1: A + v / 2
B_1: B + v / 2
C_1: C + (5 / 6 * v)
D_1: D + (5 / 6 * v)
a, M, N, O, FlächeABMD, Fläche2, Fläche3, Fläche1, FlächeBMNO, KanteAB, KanteDM, Kante2, Kante1, KanteBM, KanteMN, KanteNO, KanteBO = Prism(Polygon[A, D, D_1, A_1], B)