From 83ebe72479fdf541ed05e82460eecac3af31d25d Mon Sep 17 00:00:00 2001
From: Thibaut Horel <thibaut.horel@gmail.com>
Date: Tue, 28 Oct 2014 00:40:01 -0400
Subject: [ps5] Add code to compute integrality gap

---
 ps5/gap.py | 11 +++++++++++
 1 file changed, 11 insertions(+)
 create mode 100644 ps5/gap.py

(limited to 'ps5')

diff --git a/ps5/gap.py b/ps5/gap.py
new file mode 100644
index 0000000..45925ba
--- /dev/null
+++ b/ps5/gap.py
@@ -0,0 +1,11 @@
+from math import cos, pi
+
+
+def gap(n):
+    qp = n - 1
+    sdp = n * (1 - cos(2 * pi / n * (n / 2))) / 2.
+    return qp / sdp
+
+
+for i in xrange(3, 1000, 2):
+    print "{0:3}: {1}".format(i, gap(i))
-- 
cgit v1.2.3-70-g09d2