G.393 rectangle contained rational — Code

X.19

The rectangle contained by rational segments commensurable in length is rational.
A B C D
"""Programmatic construction for Book X, Proposition 19 (X.19)."""

        from __future__ import annotations

        from geometor.model import Model


        def construct_g1() -> Model:
            """Construct the figure 0051.g1 using geometor.model."""
            model = Model("X.19")

            # Vertices of square AD (A, B, D, D_l) and rectangle AC (A, B, C, C_l)
            A = model.set_point(25.0, -141.0, classes=["given"], ID="A")
            B = model.set_point(151.0, -141.0, classes=["given"], ID="B")
            C = model.set_point(151.0, -58.5, classes=["given"], ID="C")
            D = model.set_point(151.0, -16.0, ID="D")
            C_l = model.set_point(25.0, -58.5, ID="C_l")
            D_l = model.set_point(25.0, -16.0, ID="D_l")

            # Base lines
            model.construct_line(A, B)
            model.construct_line(B, D)
            model.construct_line(D_l, D)
            model.construct_line(A, D_l)
            model.construct_line(C_l, C)

            # Segments
            model.set_segment(A, B)
            model.set_segment(B, C)
            model.set_segment(C, D)
            model.set_segment(B, D)
            model.set_segment(C_l, C)
            model.set_segment(D_l, D)
            model.set_segment(A, C_l)
            model.set_segment(C_l, D_l)
            model.set_segment(A, D_l)

            # Section on BD
            model.set_section([B, C, D])

            # Polygons: rectangle AC and square AD
            model.set_polygon([A, B, C, C_l])
            model.set_polygon([A, B, D, D_l])
            model.set_polygon([C_l, C, D, D_l])

            return model


        if __name__ == "__main__":
            construct_g1()