From fde57297b3efbef28d09a532e1b3895d2b2ad917 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Patrik=20M=C3=BCller?= <patrik.mueller@ost.ch>
Date: Thu, 14 Jul 2022 15:03:28 +0200
Subject: Correct Makefile, add text to gamma.tex, separate python-scripts for
 each image

---
 buch/papers/laguerre/scripts/rel_error_mirror.py | 28 ++++++++++++++++++++++++
 1 file changed, 28 insertions(+)
 create mode 100644 buch/papers/laguerre/scripts/rel_error_mirror.py

(limited to 'buch/papers/laguerre/scripts/rel_error_mirror.py')

diff --git a/buch/papers/laguerre/scripts/rel_error_mirror.py b/buch/papers/laguerre/scripts/rel_error_mirror.py
new file mode 100644
index 0000000..05e68e4
--- /dev/null
+++ b/buch/papers/laguerre/scripts/rel_error_mirror.py
@@ -0,0 +1,28 @@
+if __name__ == "__main__":
+    import matplotlib.pyplot as plt
+    import numpy as np
+    import scipy.special
+
+    import gamma_approx as ga
+
+    xmin = -15
+    xmax = 15
+    ns = np.arange(2, 12, 2)
+    ylim = np.array([-11, 1])
+    x = np.linspace(xmin + ga.EPSILON, xmax - ga.EPSILON, 400)
+    gamma = scipy.special.gamma(x)
+    fig, ax = plt.subplots(num=1, clear=True, constrained_layout=True, figsize=(5, 2.5))
+    for n in ns:
+        gamma_lag = ga.eval_laguerre_gamma(x, n=n, func="mirror")
+        rel_err = ga.calc_rel_error(gamma, gamma_lag)
+        ax.semilogy(x, np.abs(rel_err), label=f"$n={n}$")
+    ax.set_xlim(x[0], x[-1])
+    ax.set_ylim(*(10.0 ** ylim))
+    ax.set_xticks(np.arange(xmin, xmax + ga.EPSILON, 5))
+    ax.set_xticks(np.arange(xmin, xmax), minor=True)
+    ax.set_yticks(10.0 ** np.arange(*ylim, 2))
+    ax.set_xlabel(r"$z$")
+    # ax.set_ylabel("Relativer Fehler")
+    ax.legend(ncol=1, loc="upper left", fontsize=ga.fontsize)
+    ax.grid(1, "both")
+    fig.savefig(f"{ga.img_path}/rel_error_mirror.pgf")
-- 
cgit v1.2.1