Add script to plot gamma graphs>

1 files changed, 19 insertions, 0 deletions

diff --git a/hacks/gamma_graphs.py b/hacks/gamma_graphs.py
new file mode 100644
index 000000000..166a03828
--- /dev/null
+++ b/hacks/gamma_graphs.py
@@ -0,0 +1,19 @@
+#!/usr/bin/python
+
+import matplotlib.pyplot as plt
+import numpy as np
+
+x = np.logspace(-3, 0, 100)
+
+plt.loglog(x, pow(x, 2.2))
+plt.loglog(x, pow(x, 2.4))
+
+linearised = []
+for xx in x:
+    if xx > 0.04045:
+        linearised.append(pow((xx + 0.055) / 1.055, 2.4))
+    else:
+        linearised.append(xx / 12.92)
+
+plt.loglog(x, linearised)
+plt.show()