diff options
author | Erik Löffler <100943759+erik-loeffler@users.noreply.github.com> | 2022-08-19 16:31:59 +0200 |
---|---|---|
committer | Erik Löffler <100943759+erik-loeffler@users.noreply.github.com> | 2022-08-19 16:31:59 +0200 |
commit | 73181c9248434eff236418fd95c7d8bb97b792b9 (patch) | |
tree | 3474048fafecda4dff38f0115096381e090ba535 /buch/papers/kugel/figures/tikz/curvature-1d.py | |
parent | Commiting again. (diff) | |
parent | Merge pull request #8 from haddoucher/sturmliouville/redabranch (diff) | |
download | SeminarSpezielleFunktionen-73181c9248434eff236418fd95c7d8bb97b792b9.tar.gz SeminarSpezielleFunktionen-73181c9248434eff236418fd95c7d8bb97b792b9.zip |
Merge remote-tracking branch 'origin/master' into sturmliouville/erik-branch
Diffstat (limited to 'buch/papers/kugel/figures/tikz/curvature-1d.py')
-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") |