aboutsummaryrefslogtreecommitdiffstats
path: root/buch/papers/kugel
diff options
context:
space:
mode:
authorAndreas Müller <andreas.mueller@ost.ch>2022-05-21 20:42:47 +0200
committerAndreas Müller <andreas.mueller@ost.ch>2022-05-21 20:42:47 +0200
commitd8d6a61a2ab45d9171a93e4a72d254a3ed5ef87f (patch)
treede63cee7199572347f97e0614dd4b05e5ea65469 /buch/papers/kugel
parentadd image code (diff)
downloadSeminarSpezielleFunktionen-d8d6a61a2ab45d9171a93e4a72d254a3ed5ef87f.tar.gz
SeminarSpezielleFunktionen-d8d6a61a2ab45d9171a93e4a72d254a3ed5ef87f.zip
fix some bugs
Diffstat (limited to 'buch/papers/kugel')
-rw-r--r--buch/papers/kugel/images/curvature.pov5
-rw-r--r--buch/papers/kugel/images/curvgraph.m8
2 files changed, 7 insertions, 6 deletions
diff --git a/buch/papers/kugel/images/curvature.pov b/buch/papers/kugel/images/curvature.pov
index 9dbaa86..3b15d77 100644
--- a/buch/papers/kugel/images/curvature.pov
+++ b/buch/papers/kugel/images/curvature.pov
@@ -72,7 +72,8 @@ arrow(<0,0,-2.1>, <0,0,2.1>, 0.01, White)
#include "curvature.inc"
#declare sigma = 1;
-#declare N0 = 0.5;
+#declare s = 1.4;
+#declare N0 = 0.4;
#declare funktion = function(r) {
(exp(-r*r/(sigma*sigma)) / sigma
-
@@ -80,7 +81,7 @@ arrow(<0,0,-2.1>, <0,0,2.1>, 0.01, White)
};
#declare hypot = function(xx, yy) { sqrt(xx*xx+yy*yy) };
-#declare Funktion = function(x,y) { funktion(hypot(x+1,y)) - funktion(hypot(x-1,y)) };
+#declare Funktion = function(x,y) { funktion(hypot(x+s,y)) - funktion(hypot(x-s,y)) };
#macro punkt(xx,yy)
<xx, Funktion(xx, yy), yy>
#end
diff --git a/buch/papers/kugel/images/curvgraph.m b/buch/papers/kugel/images/curvgraph.m
index b83c877..75effd6 100644
--- a/buch/papers/kugel/images/curvgraph.m
+++ b/buch/papers/kugel/images/curvgraph.m
@@ -11,7 +11,7 @@ global sigma2;
sigma2 = 1;
global s;
-s = 1;
+s = 1.4;
global cmax;
cmax = 0.9;
@@ -40,7 +40,7 @@ end
global N0;
N0 = f0(0)
-N0 = 0.5;
+N0 = 0.4;
function retval = f1(x,y)
global N0;
@@ -115,13 +115,13 @@ endfunction
fn = fopen("curvature.inc", "w");
punkte = zeros(4,3);
-for ix = (0:xsteps)
+for ix = (0:xsteps-1)
x = xmin + ix * hx;
punkte(1,1) = x;
punkte(2,1) = x;
punkte(3,1) = x + hx;
punkte(4,1) = x + hx;
- for iy = (0:ysteps)
+ for iy = (0:ysteps-1)
y = ymin + iy * hy;
punkte(1,2) = y;
punkte(2,2) = y + hy;