diff options
author | LordMcFungus <mceagle117@gmail.com> | 2022-07-22 21:28:45 +0200 |
---|---|---|
committer | GitHub <noreply@github.com> | 2022-07-22 21:28:45 +0200 |
commit | 23f17598c1742c70f442b94044a20aa821022c5a (patch) | |
tree | a945540ee6a4e86b37df2f01e3a91584b4797c4f /buch/papers/laguerre/scripts/rel_error_mirror.py | |
parent | Merge pull request #2 from AndreasFMueller/master (diff) | |
parent | Merge pull request #25 from JODBaer/master (diff) | |
download | SeminarSpezielleFunktionen-23f17598c1742c70f442b94044a20aa821022c5a.tar.gz SeminarSpezielleFunktionen-23f17598c1742c70f442b94044a20aa821022c5a.zip |
Merge pull request #3 from AndreasFMueller/master
update
Diffstat (limited to '')
-rw-r--r-- | buch/papers/laguerre/scripts/rel_error_mirror.py | 38 |
1 files changed, 38 insertions, 0 deletions
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..7348d5e --- /dev/null +++ b/buch/papers/laguerre/scripts/rel_error_mirror.py @@ -0,0 +1,38 @@ +if __name__ == "__main__": + import matplotlib as mpl + import matplotlib.pyplot as plt + import numpy as np + import scipy.special + + import gamma_approx as ga + + mpl.rcParams.update( + { + "mathtext.fontset": "stix", + "font.family": "serif", + "font.serif": "TeX Gyre Termes", + } + ) + + 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") + fig.savefig(f"{ga.img_path}/rel_error_mirror.pdf") |