<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">Le 29/04/2016 15:35, Jens Simon Strom a
écrit :<br>
</div>
<blockquote cite="mid:57236325.609@hslmg.de" type="cite">
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
<pre style="font-family:Courier New;font-style:normal;font-size:12.0;"><span style="color:rgb(50,185,185);">Hallo Rafael,
your script is perfect, mine isn't. Where is the bug?</span></pre>
</blockquote>
in contourd2di() called by contour2d().<br>
<br>
As a workaround, you may use:<br>
<pre style="font-family:Monospaced;font-style:normal;font-size:14.0;"><span style="color:rgb(0,0,0);">R</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">1.</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(100,174,100);font-style:italic;">//Radius Erde</span>
<span style="color:rgb(0,0,0);">a</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">50</span><span style="color:rgb(0,0,0);">;</span><span style="color:rgb(100,174,100);font-style:italic;">//Abstand Mond - Erde</span>
<span style="color:rgb(0,0,0);">g0</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(100,174,100);font-style:italic;">//Schwerebeschleunigung des Mondes im Zentrum der Erde</span>
<span style="color:rgb(0,0,0);">VLev</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(74,85,219);">[</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(188,143,143);">0.02</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(188,143,143);">0.002</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(188,143,143);">0.02</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(100,174,100);font-style:italic;">//Werte der Äquipotenziallinien-Kandidaten </span>
<span style="color:rgb(0,0,0);">n</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">400</span><span style="color:rgb(0,0,0);">;</span><span style="color:rgb(100,174,100);font-style:italic;">//Anzahl der x- und y-Werte</span>
<span style="color:rgb(174,92,176);">clf</span>
<span style="color:rgb(0,0,0);">f</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(174,92,176);">gcf</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">f</span><span style="color:rgb(92,92,92);">.</span><span style="color:rgb(170,170,170);">color_map</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(50,185,185);">zeros</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">100</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">3</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(100,174,100);font-style:italic;">//Colormap</span>
<span style="color:rgb(0,0,0);">
an</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(174,92,176);">linspace</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">0</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(218,112,214);">%pi</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">400</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">x</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(0,0,0);">a</span><span style="color:rgb(92,92,92);">+</span><span style="color:rgb(0,0,0);">R</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(50,185,185);">sin</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">an</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">y</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(0,0,0);">R</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(50,185,185);">cos</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">an</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(174,92,176);">plot</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">x</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">y</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">-w</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(0,0,0);">ca</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(174,92,176);">gca</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">ca</span><span style="color:rgb(92,92,92);">.</span><span style="color:rgb(170,170,170);">isoview</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">on</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">p</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(0,0,0);">ca</span><span style="color:rgb(92,92,92);">.</span><span style="color:rgb(170,170,170);">children</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">.</span><span style="color:rgb(170,170,170);">children</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">p</span><span style="color:rgb(92,92,92);">.</span><span style="color:rgb(170,170,170);">thickness</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(188,143,143);">4</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">p</span><span style="color:rgb(92,92,92);">.</span><span style="color:rgb(170,170,170);">data</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">3</span><span style="color:rgb(74,85,219);">)</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">0.1</span><span style="color:rgb(0,0,0);">;
</span>
<span style="color:rgb(0,0,0);">x</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(174,92,176);">linspace</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">a</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(0,0,0);">R</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">a</span><span style="color:rgb(92,92,92);">+</span><span style="color:rgb(0,0,0);">R</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">n</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">y</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(174,92,176);">linspace</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(0,0,0);">R</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">R</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">n</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(100,174,100);font-style:italic;">//Erdgebiet</span>
<span style="color:rgb(74,85,219);">[</span><span style="color:rgb(0,0,0);">X</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">Y</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(174,92,176);">ndgrid</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">x</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">y</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">V</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(0,0,0);">g0</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">X</span><span style="color:rgb(92,92,92);">+</span><span style="color:rgb(0,0,0);">a</span><span style="color:rgb(92,92,92);">^</span><span style="color:rgb(188,143,143);">2.</span><span style="color:rgb(92,92,92);">./</span><span style="color:rgb(50,185,185);">sqrt</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">X</span><span style="color:rgb(92,92,92);">.^</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(92,92,92);">+</span><span style="color:rgb(0,0,0);">Y</span><span style="color:rgb(92,92,92);">.^</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">-</span><spa
n style="color:rgb(188,143,143);">2</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">a</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">pos</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">X</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(0,0,0);">a</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">.^</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(92,92,92);">+</span><span style="color:rgb(0,0,0);">Y</span><span style="color:rgb(92,92,92);">.^</span><span style="color:rgb(188,143,143);">2</span> <span style="color:rgb(92,92,92);">></span> <span style="color:rgb(0,0,0);">R</span><span style="color:rgb(92,92,92);">^</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">V</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">pos</span><span style="color:rgb(74,85,219);">)</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(92,92,92);">-</span><span style="color:rgb(218,112,214);">%inf</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(50,185,185);">xset</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">fpf</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);"> </span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(74,85,219);">)</span> <span style="color:rgb(100,174,100);font-style:italic;">//Keine Levelanzeige, Reaktion erst nach dem 2. Lauf</span>
<span style="color:rgb(0,0,0);">philev</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(174,92,176);">contour2d</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">x</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">y</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">V</span><span style="color:rgb(0,0,0);">,</span> <span style="color:rgb(0,0,0);">VLev</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(100,174,100);font-style:italic;">//Plot und dargestellte Äquipotenziallinien</span>
<span style="color:rgb(0,0,0);"></span>... and report the bug.
Samuel
</pre>
</body>
</html>