diff options
author | Andreas Müller <andreas.mueller@ost.ch> | 2022-05-21 20:42:47 +0200 |
---|---|---|
committer | Andreas Müller <andreas.mueller@ost.ch> | 2022-05-21 20:42:47 +0200 |
commit | d8d6a61a2ab45d9171a93e4a72d254a3ed5ef87f (patch) | |
tree | de63cee7199572347f97e0614dd4b05e5ea65469 /buch/papers/kugel/images | |
parent | add image code (diff) | |
download | SeminarSpezielleFunktionen-d8d6a61a2ab45d9171a93e4a72d254a3ed5ef87f.tar.gz SeminarSpezielleFunktionen-d8d6a61a2ab45d9171a93e4a72d254a3ed5ef87f.zip |
fix some bugs
Diffstat (limited to '')
-rw-r--r-- | buch/papers/kugel/images/curvature.pov | 5 | ||||
-rw-r--r-- | buch/papers/kugel/images/curvgraph.m | 8 |
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; |