A = (0, 0, 0)
B = (5, 0, 0)
C = (0, 4, 0)
E = (5, 4, 0)
poly1, a, b, e, c = Polygon(A, B, E, C)
S = (2.5, 2, 8)
d, faceABS, faceBES, faceCES, faceACS, arêteAS, arêteBS, arêteES, arêteCS = Pyramid(poly1, S)
h = 3.01 (slider from 0 to 8, step 0.01)
H: (2.5, 2, h)
f = Plane(H, xOyPlane)
g = IntersectPath(f, faceCES)
D, F = Intersect(f, faceCES)
i = IntersectPath(f, faceACS)
G, I = Intersect(f, faceACS)
j = IntersectPath(f, faceABS)
J, K = Intersect(f, faceABS)
k = IntersectPath(f, faceBES)
L, M = Intersect(f, faceBES)
poly2, f_1, i_1, l, m, g_1 = Polygon(F, I, L, M, G)
o, T, U, V, W, Z, faceFIUT, faceILVU, faceLMWV, faceGMWZ, faceFGZT, faceTUVWZ, arêteFT, arêteIU, arêteLV, arêteMW, arêteGZ, arêteTU, arêteUV, arêteVW, arêteWZ, arêteTZ = Prism(poly2, (-h))
texte1: "Volume=" + (LaTeX[(poly2 * h)]) + " cm³"