diff options
author | Andreas Müller <andreas.mueller@ost.ch> | 2022-08-17 20:38:23 +0200 |
---|---|---|
committer | GitHub <noreply@github.com> | 2022-08-17 20:38:23 +0200 |
commit | 5ac182855b95d6cadcb748406663c4cb9875995e (patch) | |
tree | 71597f1ea3a91d6ecbaccde367a5382cfcb69081 /buch/papers/kugel/figures/tikz/curvature-1d.py | |
parent | remove superfluous file (diff) | |
parent | kugel: Minor changes (diff) | |
download | SeminarSpezielleFunktionen-5ac182855b95d6cadcb748406663c4cb9875995e.tar.gz SeminarSpezielleFunktionen-5ac182855b95d6cadcb748406663c4cb9875995e.zip |
Merge pull request #58 from NaoPross/master
kugel: Eigenvalue problem and Legendre functions
Diffstat (limited to '')
-rw-r--r-- | buch/papers/kugel/figures/tikz/curvature-1d.py | 32 |
1 files changed, 32 insertions, 0 deletions
diff --git a/buch/papers/kugel/figures/tikz/curvature-1d.py b/buch/papers/kugel/figures/tikz/curvature-1d.py new file mode 100644 index 0000000..4710fc8 --- /dev/null +++ b/buch/papers/kugel/figures/tikz/curvature-1d.py @@ -0,0 +1,32 @@ +import numpy as np +import matplotlib.pyplot as plt + + +@np.vectorize +def fn(x): + return (x ** 2) * 2 / 100 + (1 + x / 4) + np.sin(x) + +@np.vectorize +def ddfn(x): + return 2 * 5 / 100 - np.sin(x) + +x = np.linspace(0, 10, 500) +y = fn(x) +ddy = ddfn(x) + +cmap = ddy - np.min(ddy) +cmap = cmap * 1000 / np.max(cmap) + +plt.plot(x, y) +plt.plot(x, ddy) +# plt.plot(x, cmap) + +plt.show() + +fname = "curvature-1d.dat" +np.savetxt(fname, np.array([x, y, cmap]).T, delimiter=" ") + +# with open(fname, "w") as f: +# # f.write("x y cmap\n") +# for xv, yv, cv in zip(x, y, cmap): +# f.write(f"{xv} {yv} {cv}\n") |