From 56cc6c1fbae271c16c78935384b52e047cdd6f27 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Patrik=20M=C3=BCller?= Date: Thu, 19 May 2022 16:11:27 +0200 Subject: Error correction & add gamma integrand plot --- buch/papers/laguerre/scripts/integrand.py | 34 +++++++++++++++++++++++++++++++ 1 file changed, 34 insertions(+) create mode 100644 buch/papers/laguerre/scripts/integrand.py (limited to 'buch/papers/laguerre/scripts/integrand.py') diff --git a/buch/papers/laguerre/scripts/integrand.py b/buch/papers/laguerre/scripts/integrand.py new file mode 100644 index 0000000..89b9256 --- /dev/null +++ b/buch/papers/laguerre/scripts/integrand.py @@ -0,0 +1,34 @@ +#!/usr/bin/env python3 +# -*- coding:utf-8 -*- +"""Plot for integrand of gamma function with shifting terms.""" + +import os +from pathlib import Path + +import matplotlib.pyplot as plt +import numpy as np + +EPSILON = 1e-12 +xlims = np.array([-3, 3]) + +root = str(Path(__file__).parent) +img_path = f"{root}/../images" +os.makedirs(img_path, exist_ok=True) + +t = np.logspace(*xlims, 1001)[:, None] +z = np.arange(-5, 5)[None] + 0.5 + + +r = t ** z + +fig, ax = plt.subplots(num=1, clear=True, constrained_layout=True, figsize=(6, 4)) +ax.semilogx(t, r) +ax.set_xlim(*(10.**xlims)) +ax.set_ylim(1e-3, 40) +ax.set_xlabel(r"$t$") +ax.set_ylabel(r"$t^z$") +ax.grid(1, "both") +labels = [f"$z={zi:.1f}$" for zi in np.squeeze(z)] +ax.legend(labels, ncol=2, loc="upper left") +fig.savefig(f"{img_path}/integrands.pdf") +# plt.show() -- cgit v1.2.1