Polyhedron() constructor can return Empty or Universe
[linpy.git] / examples / bac2014.py
index 775be66..cc02126 100755 (executable)
@@ -1,9 +1,15 @@
 #!/usr/bin/env python3
 
+# This example is inspired from a math question in the French baccalauréat 2014,
+# consisting in computing the intersection of a plane with a line.
+
 from linpy import *
 
 x, y, z = symbols('x y z')
-DF = Eq(x, y) & Eq(z, 6 - 2*x)
-P = Eq(x + y - 2*z, 0)
+plane = Eq(x, y) & Eq(z, 6 - 2*x)
+line = Eq(x + y - 2*z, 0)
 
-print('DF âˆ© P =', DF & P)
+if __name__ == '__main__':
+    print('plane:       ', plane)
+    print('line:        ', line)
+    print('intersection:', plane & line)