diff options
author | haddoucher <reda.haddouche@ost.ch> | 2022-08-22 14:43:20 +0200 |
---|---|---|
committer | haddoucher <reda.haddouche@ost.ch> | 2022-08-22 14:43:20 +0200 |
commit | d80e30b37d3b51fc4d47229fb3e88610fbc7a476 (patch) | |
tree | 8c8808681616d4ed3cb9ac5e088c28df4139a761 /buch/papers/kugel/figures/tikz/curvature-1d.py | |
parent | Einleitung (diff) | |
download | SeminarSpezielleFunktionen-d80e30b37d3b51fc4d47229fb3e88610fbc7a476.tar.gz SeminarSpezielleFunktionen-d80e30b37d3b51fc4d47229fb3e88610fbc7a476.zip |
neuste Version
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") |