jsxgraph 1.10.0 → 1.10.1

package/distrib/docs/symbols/src/src_math_statistics.js.html

@@ -995,421 +995,428 @@
 <span class='line'>988</span>      * @memberof JXG.Math.Statistics
 <span class='line'>989</span>      *
 <span class='line'>990</span>      * @example
-<span class='line'>991</span>      *  let board = JXG.JSXGraph.initBoard('JXGbox',
-<span class='line'>992</span>      *       { boundingbox: [-1.7, .5, 30, -.03], axis: true});
+<span class='line'>991</span>      * let board = JXG.JSXGraph.initBoard('jxgbox',
+<span class='line'>992</span>      *     { boundingbox: [-1.7, .5, 30, -.03], axis: true });
 <span class='line'>993</span>      *
-<span class='line'>994</span>      *   let runs = [
-<span class='line'>995</span>      *       [0.5, 20, 'blue'],
-<span class='line'>996</span>      *       [0.7, 20, 'green'],
-<span class='line'>997</span>      *       [0.5, 40, 'red'],
-<span class='line'>998</span>      *   ]
+<span class='line'>994</span>      * let runs = [
+<span class='line'>995</span>      *     [0.5, 20, 'blue'],
+<span class='line'>996</span>      *     [0.7, 20, 'green'],
+<span class='line'>997</span>      *     [0.5, 40, 'red'],
+<span class='line'>998</span>      * ];
 <span class='line'>999</span>      *
-<span class='line'>1000</span>      *   let labelY = .4
-<span class='line'>1001</span>      *   runs.forEach((run,i) => {
-<span class='line'>1002</span>      *       board.create('segment',[[7,labelY-(i/50)],[9,labelY-(i/50)]],{strokeColor:run[2]})
-<span class='line'>1003</span>      *       board.create('text',[10,labelY-(i/50),`p=${run[0]}, n=${run[1]}`])
+<span class='line'>1000</span>      * let labelY = .4;
+<span class='line'>1001</span>      * runs.forEach((run, i) => {
+<span class='line'>1002</span>      *     board.create('segment', [[7, labelY - (i / 50)], [9, labelY - (i / 50)]], { strokeColor: run[2] });
+<span class='line'>1003</span>      *     board.create('text', [10, labelY - (i / 50), `p=${run[0]}, n=${run[1]}`]);
 <span class='line'>1004</span>      *
-<span class='line'>1005</span>      *       let x = Array(50000).fill(0).map(() => JXG.Math.Statistics.randomBinomial(run[1],run[0]))
-<span class='line'>1006</span>      *       let res = JXG.Math.Statistics.histogram(x, { bins: 40, density: true, cumulative: false, range: [0, 40] });
-<span class='line'>1007</span>      *       board.create('curve', [res[1], res[0]], { strokeColor: run[2]});
-<span class='line'>1008</span>      *   })
-<span class='line'>1009</span>      * }
-<span class='line'>1010</span>      *
-<span class='line'>1011</span>      *
-<span class='line'>1012</span>      * &lt;/pre>&lt;div id="JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302-3" class="jxgbox" style="width: 300px; height: 300px;">&lt;/div>
-<span class='line'>1013</span>      * &lt;script type="text/javascript">
-<span class='line'>1014</span>      * {
-<span class='line'>1015</span>      *  let board = JXG.JSXGraph.initBoard('JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302-3',
-<span class='line'>1016</span>      *       { boundingbox: [-1.7, .5, 30, -.03], axis: true});
-<span class='line'>1017</span>      *
-<span class='line'>1018</span>      *   let runs = [
-<span class='line'>1019</span>      *       [0.5, 20, 'blue'],
-<span class='line'>1020</span>      *       [0.7, 20, 'green'],
-<span class='line'>1021</span>      *       [0.5, 40, 'red'],
-<span class='line'>1022</span>      *   ]
-<span class='line'>1023</span>      *
-<span class='line'>1024</span>      *   let labelY = .4
-<span class='line'>1025</span>      *   runs.forEach((run,i) => {
-<span class='line'>1026</span>      *       board.create('segment',[[7,labelY-(i/50)],[9,labelY-(i/50)]],{strokeColor:run[2]})
-<span class='line'>1027</span>      *       board.create('text',[10,labelY-(i/50),`p=${run[0]}, n=${run[1]}`])
-<span class='line'>1028</span>      *
-<span class='line'>1029</span>      *       let x = Array(50000).fill(0).map(() => JXG.Math.Statistics.randomBinomial(run[1],run[0]))
-<span class='line'>1030</span>      *       let res = JXG.Math.Statistics.histogram(x, { bins: 40, density: true, cumulative: false, range: [0, 40] });
-<span class='line'>1031</span>      *       board.create('curve', [res[1], res[0]], { strokeColor: run[2]});
-<span class='line'>1032</span>      *   })
-<span class='line'>1033</span>      * }
-<span class='line'>1034</span>      * &lt;/script>&lt;pre>
-<span class='line'>1035</span>      *
-<span class='line'>1036</span>      */</span><span class="WHIT">
-<span class='line'>1037</span> </span><span class="WHIT">    </span><span class="NAME">randomBinomial</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">n</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1038</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">y</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">c</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1039</span> </span><span class="WHIT">            </span><span class="NAME">a</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">b</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">N1</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1040</span> 
-<span class='line'>1041</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1042</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">NaN</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1043</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1005</span>      *     let x = Array(50000).fill(0).map(() => JXG.Math.Statistics.randomBinomial(run[1], run[0]));
+<span class='line'>1006</span>      *     let res = JXG.Math.Statistics.histogram(x, {
+<span class='line'>1007</span>      *         bins: 40,
+<span class='line'>1008</span>      *         density: true,
+<span class='line'>1009</span>      *         cumulative: false,
+<span class='line'>1010</span>      *         range: [0, 40]
+<span class='line'>1011</span>      *     });
+<span class='line'>1012</span>      *     board.create('curve', [res[1], res[0]], { strokeColor: run[2] });
+<span class='line'>1013</span>      * });
+<span class='line'>1014</span>      *
+<span class='line'>1015</span>      *
+<span class='line'>1016</span>      * &lt;/pre>&lt;div id="JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302-3" class="jxgbox" style="width: 300px; height: 300px;">&lt;/div>
+<span class='line'>1017</span>      * &lt;script type="text/javascript">
+<span class='line'>1018</span>      * {
+<span class='line'>1019</span>      *  let board = JXG.JSXGraph.initBoard('JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302-3',
+<span class='line'>1020</span>      *       { boundingbox: [-1.7, .5, 30, -.03], axis: true});
+<span class='line'>1021</span>      *
+<span class='line'>1022</span>      *   let runs = [
+<span class='line'>1023</span>      *       [0.5, 20, 'blue'],
+<span class='line'>1024</span>      *       [0.7, 20, 'green'],
+<span class='line'>1025</span>      *       [0.5, 40, 'red'],
+<span class='line'>1026</span>      *   ]
+<span class='line'>1027</span>      *
+<span class='line'>1028</span>      *   let labelY = .4
+<span class='line'>1029</span>      *   runs.forEach((run,i) => {
+<span class='line'>1030</span>      *       board.create('segment',[[7,labelY-(i/50)],[9,labelY-(i/50)]],{strokeColor:run[2]})
+<span class='line'>1031</span>      *       board.create('text',[10,labelY-(i/50),`p=${run[0]}, n=${run[1]}`])
+<span class='line'>1032</span>      *
+<span class='line'>1033</span>      *       let x = Array(50000).fill(0).map(() => JXG.Math.Statistics.randomBinomial(run[1],run[0]))
+<span class='line'>1034</span>      *       let res = JXG.Math.Statistics.histogram(x, { bins: 40, density: true, cumulative: false, range: [0, 40] });
+<span class='line'>1035</span>      *       board.create('curve', [res[1], res[0]], { strokeColor: run[2]});
+<span class='line'>1036</span>      *   })
+<span class='line'>1037</span>      * }
+<span class='line'>1038</span>      * &lt;/script>&lt;pre>
+<span class='line'>1039</span>      *
+<span class='line'>1040</span>      */</span><span class="WHIT">
+<span class='line'>1041</span> </span><span class="WHIT">    </span><span class="NAME">randomBinomial</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">n</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1042</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">y</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">c</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1043</span> </span><span class="WHIT">            </span><span class="NAME">a</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">b</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">N1</span><span class="PUNC">;</span><span class="WHIT">
 <span class='line'>1044</span> 
-<span class='line'>1045</span> </span><span class="WHIT">        </span><span class="COMM">// Edge cases</span><span class="WHIT">
-<span class='line'>1046</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1047</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1048</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1049</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1050</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">n</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1051</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1052</span> 
-<span class='line'>1053</span> </span><span class="WHIT">        </span><span class="COMM">// Now, we can assume 0 &lt; p &lt; 1.</span><span class="WHIT">
-<span class='line'>1054</span> 
-<span class='line'>1055</span> </span><span class="WHIT">        </span><span class="COMM">// Fast path for common cases</span><span class="WHIT">
-<span class='line'>1056</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1057</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1058</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1059</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1060</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="PUNC">(</span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">?</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">:</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1061</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1062</span> 
-<span class='line'>1063</span> </span><span class="WHIT">        </span><span class="COMM">// Exploit symmetry</span><span class="WHIT">
-<span class='line'>1064</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NUMB">0.5</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1065</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">this.randomBinomial</span><span class="PUNC">(</span><span class="NAME">n</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1066</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1067</span> 
-<span class='line'>1068</span> </span><span class="WHIT">        </span><span class="COMM">// General case: n > 1, p &lt;= 0.5</span><span class="WHIT">
-<span class='line'>1069</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">100</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1070</span> </span><span class="WHIT">            </span><span class="COMM">// n small:</span><span class="WHIT">
-<span class='line'>1071</span> </span><span class="WHIT">            </span><span class="COMM">// Algorithm BG (Devroye) from:</span><span class="WHIT">
-<span class='line'>1072</span> </span><span class="WHIT">            </span><span class="COMM">// https://dl.acm.org/doi/pdf/10.1145/42372.42381</span><span class="WHIT">
-<span class='line'>1073</span> </span><span class="WHIT">            </span><span class="COMM">// Time O(np) so suitable for np small only.</span><span class="WHIT">
-<span class='line'>1074</span> </span><span class="WHIT">            </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1075</span> </span><span class="WHIT">            </span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1076</span> 
-<span class='line'>1077</span> </span><span class="WHIT">            </span><span class="NAME">c</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.log</span><span class="PUNC">(</span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1078</span> </span><span class="WHIT">            </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">c</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1079</span> </span><span class="WHIT">                </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1080</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1081</span> 
-<span class='line'>1082</span> </span><span class="WHIT">            </span><span class="KEYW">do</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1083</span> </span><span class="WHIT">                </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1084</span> </span><span class="WHIT">                </span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.floor</span><span class="PUNC">(</span><span class="NAME">Math.log</span><span class="PUNC">(</span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">c</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1085</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">while</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">n</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1086</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1087</span> </span><span class="WHIT">            </span><span class="COMM">// n large:</span><span class="WHIT">
-<span class='line'>1088</span> </span><span class="WHIT">            </span><span class="COMM">// Knuth, TAOCP 2, p 131</span><span class="WHIT">
-<span class='line'>1089</span> </span><span class="WHIT">            </span><span class="NAME">a</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">Math.floor</span><span class="PUNC">(</span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NUMB">0.5</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1090</span> </span><span class="WHIT">            </span><span class="NAME">b</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">a</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1091</span> </span><span class="WHIT">            </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">this.randomBeta</span><span class="PUNC">(</span><span class="NAME">a</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">b</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1092</span> </span><span class="WHIT">            </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">>=</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1093</span> </span><span class="WHIT">                </span><span class="NAME">N1</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">this.randomBinomial</span><span class="PUNC">(</span><span class="NAME">a</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1094</span> </span><span class="WHIT">                </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">N1</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1095</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1096</span> </span><span class="WHIT">                </span><span class="NAME">N1</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">this.randomBinomial</span><span class="PUNC">(</span><span class="NAME">b</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1097</span> </span><span class="WHIT">                </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">a</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">N1</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1098</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1099</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1100</span> </span><span class="WHIT">        </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1101</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1102</span> 
-<span class='line'>1103</span> </span><span class="WHIT">    </span><span class="COMM">/**
-<span class='line'>1104</span>      * Generate values for a random variable in geometric distribution with probability &lt;i>p&lt;/i>.
-<span class='line'>1105</span>      * See &lt;a href="https://en.wikipedia.org/wiki/Geometric_distribution">https://en.wikipedia.org/wiki/Geometric_distribution&lt;/a>.
-<span class='line'>1106</span>      *
-<span class='line'>1107</span>      * @param {Number} p (0 &lt;= p &lt;= 1)
-<span class='line'>1108</span>      * @returns Number
-<span class='line'>1109</span>      * @memberof JXG.Math.Statistics
-<span class='line'>1110</span>      */</span><span class="WHIT">
-<span class='line'>1111</span> </span><span class="WHIT">    </span><span class="NAME">randomGeometric</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1112</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">u</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1113</span> 
-<span class='line'>1114</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1115</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">NaN</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1116</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1117</span> </span><span class="WHIT">        </span><span class="COMM">// Knuth, TAOCP 2, p 131</span><span class="WHIT">
-<span class='line'>1118</span> </span><span class="WHIT">        </span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1119</span> 
-<span class='line'>1120</span> </span><span class="WHIT">        </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">Math.ceil</span><span class="PUNC">(</span><span class="NAME">Math.log</span><span class="PUNC">(</span><span class="NAME">u</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">Math.log</span><span class="PUNC">(</span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1121</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1122</span> 
-<span class='line'>1123</span> </span><span class="WHIT">    </span><span class="COMM">/**
-<span class='line'>1124</span>      * Generate values for a random variable in Poisson distribution with mean &lt;i>mu&lt;/i>.
-<span class='line'>1125</span>      * See &lt;a href="https://en.wikipedia.org/wiki/Poisson_distribution">https://en.wikipedia.org/wiki/Poisson_distribution&lt;/a>.
-<span class='line'>1126</span>      *
-<span class='line'>1127</span>      * @param {Number} mu (0 &lt; mu)
-<span class='line'>1128</span>      * @returns Number
-<span class='line'>1129</span>      * @memberof JXG.Math.Statistics
-<span class='line'>1130</span>      */</span><span class="WHIT">
-<span class='line'>1131</span> </span><span class="WHIT">    </span><span class="NAME">randomPoisson</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">mu</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1132</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">e</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.exp</span><span class="PUNC">(</span><span class="PUNC">-</span><span class="NAME">mu</span><span class="PUNC">)</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1133</span> </span><span class="WHIT">            </span><span class="NAME">N</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1134</span> </span><span class="WHIT">            </span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1135</span> </span><span class="WHIT">            </span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1136</span> </span><span class="WHIT">            </span><span class="NAME">x</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1137</span> </span><span class="WHIT">            </span><span class="NAME">alpha</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">7</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NUMB">8</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1138</span> 
-<span class='line'>1139</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">mu</span><span class="WHIT"> </span><span class="PUNC">&lt;=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1140</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">NaN</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1141</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1045</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1046</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">NaN</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1047</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1048</span> 
+<span class='line'>1049</span> </span><span class="WHIT">        </span><span class="COMM">// Edge cases</span><span class="WHIT">
+<span class='line'>1050</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1051</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1052</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1053</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1054</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">n</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1055</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1056</span> 
+<span class='line'>1057</span> </span><span class="WHIT">        </span><span class="COMM">// Now, we can assume 0 &lt; p &lt; 1.</span><span class="WHIT">
+<span class='line'>1058</span> 
+<span class='line'>1059</span> </span><span class="WHIT">        </span><span class="COMM">// Fast path for common cases</span><span class="WHIT">
+<span class='line'>1060</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1061</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1062</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1063</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1064</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="PUNC">(</span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">?</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">:</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1065</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1066</span> 
+<span class='line'>1067</span> </span><span class="WHIT">        </span><span class="COMM">// Exploit symmetry</span><span class="WHIT">
+<span class='line'>1068</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NUMB">0.5</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1069</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">this.randomBinomial</span><span class="PUNC">(</span><span class="NAME">n</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1070</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1071</span> 
+<span class='line'>1072</span> </span><span class="WHIT">        </span><span class="COMM">// General case: n > 1, p &lt;= 0.5</span><span class="WHIT">
+<span class='line'>1073</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">100</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1074</span> </span><span class="WHIT">            </span><span class="COMM">// n small:</span><span class="WHIT">
+<span class='line'>1075</span> </span><span class="WHIT">            </span><span class="COMM">// Algorithm BG (Devroye) from:</span><span class="WHIT">
+<span class='line'>1076</span> </span><span class="WHIT">            </span><span class="COMM">// https://dl.acm.org/doi/pdf/10.1145/42372.42381</span><span class="WHIT">
+<span class='line'>1077</span> </span><span class="WHIT">            </span><span class="COMM">// Time O(np) so suitable for np small only.</span><span class="WHIT">
+<span class='line'>1078</span> </span><span class="WHIT">            </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1079</span> </span><span class="WHIT">            </span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1080</span> 
+<span class='line'>1081</span> </span><span class="WHIT">            </span><span class="NAME">c</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.log</span><span class="PUNC">(</span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1082</span> </span><span class="WHIT">            </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">c</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1083</span> </span><span class="WHIT">                </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1084</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1085</span> 
+<span class='line'>1086</span> </span><span class="WHIT">            </span><span class="KEYW">do</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1087</span> </span><span class="WHIT">                </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1088</span> </span><span class="WHIT">                </span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.floor</span><span class="PUNC">(</span><span class="NAME">Math.log</span><span class="PUNC">(</span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">c</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1089</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">while</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">n</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1090</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1091</span> </span><span class="WHIT">            </span><span class="COMM">// n large:</span><span class="WHIT">
+<span class='line'>1092</span> </span><span class="WHIT">            </span><span class="COMM">// Knuth, TAOCP 2, p 131</span><span class="WHIT">
+<span class='line'>1093</span> </span><span class="WHIT">            </span><span class="NAME">a</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">Math.floor</span><span class="PUNC">(</span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NUMB">0.5</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1094</span> </span><span class="WHIT">            </span><span class="NAME">b</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">n</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">a</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1095</span> </span><span class="WHIT">            </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">this.randomBeta</span><span class="PUNC">(</span><span class="NAME">a</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">b</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1096</span> </span><span class="WHIT">            </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">>=</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1097</span> </span><span class="WHIT">                </span><span class="NAME">N1</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">this.randomBinomial</span><span class="PUNC">(</span><span class="NAME">a</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1098</span> </span><span class="WHIT">                </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">N1</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1099</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1100</span> </span><span class="WHIT">                </span><span class="NAME">N1</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">this.randomBinomial</span><span class="PUNC">(</span><span class="NAME">b</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1101</span> </span><span class="WHIT">                </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">a</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">N1</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1102</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1103</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1104</span> </span><span class="WHIT">        </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1105</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1106</span> 
+<span class='line'>1107</span> </span><span class="WHIT">    </span><span class="COMM">/**
+<span class='line'>1108</span>      * Generate values for a random variable in geometric distribution with probability &lt;i>p&lt;/i>.
+<span class='line'>1109</span>      * See &lt;a href="https://en.wikipedia.org/wiki/Geometric_distribution">https://en.wikipedia.org/wiki/Geometric_distribution&lt;/a>.
+<span class='line'>1110</span>      *
+<span class='line'>1111</span>      * @param {Number} p (0 &lt;= p &lt;= 1)
+<span class='line'>1112</span>      * @returns Number
+<span class='line'>1113</span>      * @memberof JXG.Math.Statistics
+<span class='line'>1114</span>      */</span><span class="WHIT">
+<span class='line'>1115</span> </span><span class="WHIT">    </span><span class="NAME">randomGeometric</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1116</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">u</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1117</span> 
+<span class='line'>1118</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">p</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1119</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">NaN</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1120</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1121</span> </span><span class="WHIT">        </span><span class="COMM">// Knuth, TAOCP 2, p 131</span><span class="WHIT">
+<span class='line'>1122</span> </span><span class="WHIT">        </span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1123</span> 
+<span class='line'>1124</span> </span><span class="WHIT">        </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">Math.ceil</span><span class="PUNC">(</span><span class="NAME">Math.log</span><span class="PUNC">(</span><span class="NAME">u</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">Math.log</span><span class="PUNC">(</span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">p</span><span class="PUNC">)</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1125</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1126</span> 
+<span class='line'>1127</span> </span><span class="WHIT">    </span><span class="COMM">/**
+<span class='line'>1128</span>      * Generate values for a random variable in Poisson distribution with mean &lt;i>mu&lt;/i>.
+<span class='line'>1129</span>      * See &lt;a href="https://en.wikipedia.org/wiki/Poisson_distribution">https://en.wikipedia.org/wiki/Poisson_distribution&lt;/a>.
+<span class='line'>1130</span>      *
+<span class='line'>1131</span>      * @param {Number} mu (0 &lt; mu)
+<span class='line'>1132</span>      * @returns Number
+<span class='line'>1133</span>      * @memberof JXG.Math.Statistics
+<span class='line'>1134</span>      */</span><span class="WHIT">
+<span class='line'>1135</span> </span><span class="WHIT">    </span><span class="NAME">randomPoisson</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">mu</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1136</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">e</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.exp</span><span class="PUNC">(</span><span class="PUNC">-</span><span class="NAME">mu</span><span class="PUNC">)</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1137</span> </span><span class="WHIT">            </span><span class="NAME">N</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1138</span> </span><span class="WHIT">            </span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1139</span> </span><span class="WHIT">            </span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1140</span> </span><span class="WHIT">            </span><span class="NAME">x</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1141</span> </span><span class="WHIT">            </span><span class="NAME">alpha</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">7</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NUMB">8</span><span class="PUNC">;</span><span class="WHIT">
 <span class='line'>1142</span> 
-<span class='line'>1143</span> </span><span class="WHIT">        </span><span class="COMM">// Knuth, TAOCP 2, p 132</span><span class="WHIT">
-<span class='line'>1144</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">mu</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">10</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1145</span> </span><span class="WHIT">            </span><span class="KEYW">do</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1146</span> </span><span class="WHIT">                </span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1147</span> </span><span class="WHIT">                </span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1148</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">while</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NAME">e</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1149</span> </span><span class="WHIT">            </span><span class="NAME">N</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1150</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1151</span> </span><span class="WHIT">            </span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.floor</span><span class="PUNC">(</span><span class="NAME">alpha</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NAME">mu</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1152</span> </span><span class="WHIT">            </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">this.randomGamma</span><span class="PUNC">(</span><span class="NAME">m</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1153</span> </span><span class="WHIT">            </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">mu</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1154</span> </span><span class="WHIT">                </span><span class="NAME">N</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">this.randomPoisson</span><span class="PUNC">(</span><span class="NAME">mu</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1155</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1156</span> </span><span class="WHIT">                </span><span class="NAME">N</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">this.randomBinomial</span><span class="PUNC">(</span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">mu</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1157</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1158</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1159</span> </span><span class="WHIT">        </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">N</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1160</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1161</span> 
-<span class='line'>1162</span> </span><span class="WHIT">    </span><span class="COMM">/**
-<span class='line'>1163</span>      * Generate values for a random variable in Pareto distribution with
-<span class='line'>1164</span>      * shape &lt;i>gamma&lt;/i> and scale &lt;i>k&lt;/i>.
-<span class='line'>1165</span>      * See &lt;a href="https://en.wikipedia.org/wiki/Pareto_distribution">https://en.wikipedia.org/wiki/Pareto_distribution&lt;/a>.
-<span class='line'>1166</span>      * Method: use inverse transformation sampling.
-<span class='line'>1167</span>      *
-<span class='line'>1168</span>      * @param {Number} gamma shape (0 &lt; gamma)
-<span class='line'>1169</span>      * @param {Number} k scale (0 &lt; k &lt; x)
-<span class='line'>1170</span>      * @returns Number
-<span class='line'>1171</span>      * @memberof JXG.Math.Statistics
-<span class='line'>1172</span>      */</span><span class="WHIT">
-<span class='line'>1173</span> </span><span class="WHIT">    </span><span class="NAME">randomPareto</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">gamma</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">k</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1174</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1175</span> 
-<span class='line'>1176</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">gamma</span><span class="WHIT"> </span><span class="PUNC">&lt;=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">&lt;=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1177</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">NaN</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1178</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1179</span> </span><span class="WHIT">        </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NAME">Math.pow</span><span class="PUNC">(</span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">u</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">gamma</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1180</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1181</span> 
-<span class='line'>1182</span> </span><span class="WHIT">    </span><span class="COMM">/**
-<span class='line'>1183</span>      * Generate values for a random variable in hypergeometric distribution.
-<span class='line'>1184</span>      * Samples are drawn from a hypergeometric distribution with specified parameters, &lt;i>good&lt;/i> (ways to make a good selection),
-<span class='line'>1185</span>      * &lt;i>bad&lt;/i> (ways to make a bad selection), and &lt;i>samples&lt;/i> (number of items sampled, which is less than or equal to &lt;i>good + bad&lt;/i>).
-<span class='line'>1186</span>      * &lt;p>
-<span class='line'>1187</span>      * Naive implementation with runtime &lt;i>O(samples)&lt;/i>.
-<span class='line'>1188</span>      *
-<span class='line'>1189</span>      * @param {Number} good ways to make a good selection
-<span class='line'>1190</span>      * @param {Number} bad ways to make a bad selection
-<span class='line'>1191</span>      * @param {Number} samples number of items sampled
-<span class='line'>1192</span>      * @returns
-<span class='line'>1193</span>      * @memberof JXG.Math.Statistics
-<span class='line'>1194</span>      */</span><span class="WHIT">
-<span class='line'>1195</span> </span><span class="WHIT">    </span><span class="NAME">randomHypergeometric</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">good</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">k</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1196</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">u</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1197</span> </span><span class="WHIT">            </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1198</span> </span><span class="WHIT">            </span><span class="COMM">// kk,</span><span class="WHIT">
-<span class='line'>1199</span> </span><span class="WHIT">            </span><span class="COMM">// n = good + bad,</span><span class="WHIT">
-<span class='line'>1200</span> </span><span class="WHIT">            </span><span class="NAME">d1</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">good</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">k</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1201</span> </span><span class="WHIT">            </span><span class="NAME">d2</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.min</span><span class="PUNC">(</span><span class="NAME">good</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="PUNC">)</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1202</span> </span><span class="WHIT">            </span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">d2</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1203</span> 
-<span class='line'>1204</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">good</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NAME">good</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1205</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">NaN</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1206</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1143</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">mu</span><span class="WHIT"> </span><span class="PUNC">&lt;=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1144</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">NaN</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1145</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1146</span> 
+<span class='line'>1147</span> </span><span class="WHIT">        </span><span class="COMM">// Knuth, TAOCP 2, p 132</span><span class="WHIT">
+<span class='line'>1148</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">mu</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">10</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1149</span> </span><span class="WHIT">            </span><span class="KEYW">do</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1150</span> </span><span class="WHIT">                </span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1151</span> </span><span class="WHIT">                </span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1152</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">while</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NAME">e</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1153</span> </span><span class="WHIT">            </span><span class="NAME">N</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1154</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1155</span> </span><span class="WHIT">            </span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.floor</span><span class="PUNC">(</span><span class="NAME">alpha</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NAME">mu</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1156</span> </span><span class="WHIT">            </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">this.randomGamma</span><span class="PUNC">(</span><span class="NAME">m</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1157</span> </span><span class="WHIT">            </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">mu</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1158</span> </span><span class="WHIT">                </span><span class="NAME">N</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">this.randomPoisson</span><span class="PUNC">(</span><span class="NAME">mu</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1159</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1160</span> </span><span class="WHIT">                </span><span class="NAME">N</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">this.randomBinomial</span><span class="PUNC">(</span><span class="NAME">m</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">mu</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1161</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1162</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1163</span> </span><span class="WHIT">        </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">N</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1164</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1165</span> 
+<span class='line'>1166</span> </span><span class="WHIT">    </span><span class="COMM">/**
+<span class='line'>1167</span>      * Generate values for a random variable in Pareto distribution with
+<span class='line'>1168</span>      * shape &lt;i>gamma&lt;/i> and scale &lt;i>k&lt;/i>.
+<span class='line'>1169</span>      * See &lt;a href="https://en.wikipedia.org/wiki/Pareto_distribution">https://en.wikipedia.org/wiki/Pareto_distribution&lt;/a>.
+<span class='line'>1170</span>      * Method: use inverse transformation sampling.
+<span class='line'>1171</span>      *
+<span class='line'>1172</span>      * @param {Number} gamma shape (0 &lt; gamma)
+<span class='line'>1173</span>      * @param {Number} k scale (0 &lt; k &lt; x)
+<span class='line'>1174</span>      * @returns Number
+<span class='line'>1175</span>      * @memberof JXG.Math.Statistics
+<span class='line'>1176</span>      */</span><span class="WHIT">
+<span class='line'>1177</span> </span><span class="WHIT">    </span><span class="NAME">randomPareto</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">gamma</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">k</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1178</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1179</span> 
+<span class='line'>1180</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">gamma</span><span class="WHIT"> </span><span class="PUNC">&lt;=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">&lt;=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1181</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">NaN</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1182</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1183</span> </span><span class="WHIT">        </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NAME">Math.pow</span><span class="PUNC">(</span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">u</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">gamma</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1184</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1185</span> 
+<span class='line'>1186</span> </span><span class="WHIT">    </span><span class="COMM">/**
+<span class='line'>1187</span>      * Generate values for a random variable in hypergeometric distribution.
+<span class='line'>1188</span>      * Samples are drawn from a hypergeometric distribution with specified parameters, &lt;i>good&lt;/i> (ways to make a good selection),
+<span class='line'>1189</span>      * &lt;i>bad&lt;/i> (ways to make a bad selection), and &lt;i>samples&lt;/i> (number of items sampled, which is less than or equal to &lt;i>good + bad&lt;/i>).
+<span class='line'>1190</span>      * &lt;p>
+<span class='line'>1191</span>      * Naive implementation with runtime &lt;i>O(samples)&lt;/i>.
+<span class='line'>1192</span>      *
+<span class='line'>1193</span>      * @param {Number} good ways to make a good selection
+<span class='line'>1194</span>      * @param {Number} bad ways to make a bad selection
+<span class='line'>1195</span>      * @param {Number} samples number of items sampled
+<span class='line'>1196</span>      * @returns
+<span class='line'>1197</span>      * @memberof JXG.Math.Statistics
+<span class='line'>1198</span>      */</span><span class="WHIT">
+<span class='line'>1199</span> </span><span class="WHIT">    </span><span class="NAME">randomHypergeometric</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">good</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">k</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1200</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">u</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1201</span> </span><span class="WHIT">            </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1202</span> </span><span class="WHIT">            </span><span class="COMM">// kk,</span><span class="WHIT">
+<span class='line'>1203</span> </span><span class="WHIT">            </span><span class="COMM">// n = good + bad,</span><span class="WHIT">
+<span class='line'>1204</span> </span><span class="WHIT">            </span><span class="NAME">d1</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">good</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">k</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1205</span> </span><span class="WHIT">            </span><span class="NAME">d2</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.min</span><span class="PUNC">(</span><span class="NAME">good</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="PUNC">)</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1206</span> </span><span class="WHIT">            </span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">d2</span><span class="PUNC">;</span><span class="WHIT">
 <span class='line'>1207</span> 
-<span class='line'>1208</span> </span><span class="WHIT">        </span><span class="COMM">// Naive method</span><span class="WHIT">
-<span class='line'>1209</span> </span><span class="WHIT">        </span><span class="COMM">// kk = Math.min(k, n - k);</span><span class="WHIT">
-<span class='line'>1210</span> </span><span class="WHIT">        </span><span class="COMM">// for (i = 0; i &lt; k; i ++) {</span><span class="WHIT">
-<span class='line'>1211</span> </span><span class="WHIT">        </span><span class="COMM">//     u = Math.random();</span><span class="WHIT">
-<span class='line'>1212</span> </span><span class="WHIT">        </span><span class="COMM">//     if (n * u &lt;= good) {</span><span class="WHIT">
-<span class='line'>1213</span> </span><span class="WHIT">        </span><span class="COMM">//         x += 1;</span><span class="WHIT">
-<span class='line'>1214</span> </span><span class="WHIT">        </span><span class="COMM">//         if (x === good) {</span><span class="WHIT">
-<span class='line'>1215</span> </span><span class="WHIT">        </span><span class="COMM">//             return x;</span><span class="WHIT">
-<span class='line'>1216</span> </span><span class="WHIT">        </span><span class="COMM">//         }</span><span class="WHIT">
-<span class='line'>1217</span> </span><span class="WHIT">        </span><span class="COMM">//         good -= 1;</span><span class="WHIT">
-<span class='line'>1218</span> </span><span class="WHIT">        </span><span class="COMM">//     }</span><span class="WHIT">
-<span class='line'>1219</span> </span><span class="WHIT">        </span><span class="COMM">//     n -= 1;</span><span class="WHIT">
-<span class='line'>1220</span> </span><span class="WHIT">        </span><span class="COMM">// }</span><span class="WHIT">
-<span class='line'>1221</span> </span><span class="WHIT">        </span><span class="COMM">// return x;</span><span class="WHIT">
-<span class='line'>1222</span> 
-<span class='line'>1223</span> </span><span class="WHIT">        </span><span class="COMM">// Implementation from</span><span class="WHIT">
-<span class='line'>1224</span> </span><span class="WHIT">        </span><span class="COMM">// Monte Carlo by George S. Fishman</span><span class="WHIT">
-<span class='line'>1225</span> </span><span class="WHIT">        </span><span class="COMM">// https://link.springer.com/book/10.1007/978-1-4757-2553-7</span><span class="WHIT">
-<span class='line'>1226</span> </span><span class="WHIT">        </span><span class="COMM">// page 218</span><span class="WHIT">
-<span class='line'>1227</span> </span><span class="WHIT">        </span><span class="COMM">//</span><span class="WHIT">
-<span class='line'>1228</span> </span><span class="WHIT">        </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">k</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1229</span> </span><span class="WHIT">        </span><span class="KEYW">while</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1230</span> </span><span class="WHIT">            </span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1231</span> </span><span class="WHIT">            </span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.floor</span><span class="PUNC">(</span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">d1</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">)</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1232</span> </span><span class="WHIT">            </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1233</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1234</span> </span><span class="WHIT">        </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">d2</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">y</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1235</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">good</span><span class="WHIT"> </span><span class="PUNC">&lt;=</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1236</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1237</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1238</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1239</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1240</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1241</span> 
-<span class='line'>1242</span> </span><span class="WHIT">    </span><span class="COMM">/**
-<span class='line'>1243</span>      * Compute the histogram of a dataset.
-<span class='line'>1244</span>      * Optional parameters can be supplied through a JavaScript object
-<span class='line'>1245</span>      * with the following default values:
-<span class='line'>1246</span>      * &lt;pre>
-<span class='line'>1247</span>      * {
-<span class='line'>1248</span>      *   bins: 10,          // Number of bins
-<span class='line'>1249</span>      *   range: false,      // false or array. The lower and upper range of the bins.
-<span class='line'>1250</span>      *                      // If not provided, range is simply [min(x), max(x)].
-<span class='line'>1251</span>      *                      // Values outside the range are ignored.
-<span class='line'>1252</span>      *   density: false,    // If true, normalize the counts by dividing by sum(counts)
-<span class='line'>1253</span>      *   cumulative: false
-<span class='line'>1254</span>      * }
-<span class='line'>1255</span>      * &lt;/pre>
-<span class='line'>1256</span>      * The function returns an array containing two arrays. The first array is of length bins+1
-<span class='line'>1257</span>      * containing the start values of the bins. The last entry contains the end values of the last bin.
-<span class='line'>1258</span>      * &lt;p>
-<span class='line'>1259</span>      * The second array contains the counts of each bin.
-<span class='line'>1260</span>      * @param {Array} x
-<span class='line'>1261</span>      * @param {Object} opt Optional parameters
-<span class='line'>1262</span>      * @returns Array [bin, counts] Array bins contains start values of bins, array counts contains
-<span class='line'>1263</span>      * the number of entries of x which are contained in each bin.
-<span class='line'>1264</span>      * @memberof JXG.Math.Statistics
-<span class='line'>1265</span>      *
-<span class='line'>1266</span>      * @example
-<span class='line'>1267</span>      *  let board = JXG.JSXGraph.initBoard('jxgbox',
-<span class='line'>1268</span>      *       { boundingbox: [-1.7, .5, 20, -.03], axis: true});
-<span class='line'>1269</span>      *  let board2 = JXG.JSXGraph.initBoard('jxgbox2',
-<span class='line'>1270</span>      *       { boundingbox: [-1.6, 1.1, 20, -.06], axis: true});
-<span class='line'>1271</span>      *
-<span class='line'>1272</span>      *   let runs = [
-<span class='line'>1273</span>      *       [0.5, 1.0, 'brown'],
-<span class='line'>1274</span>      *       [1.0, 2.0, 'red'],
-<span class='line'>1275</span>      *       [2.0, 2.0, 'orange'],
-<span class='line'>1276</span>      *       [3.0, 2.0, 'yellow'],
-<span class='line'>1277</span>      *       [5.0, 1.0, 'green'],
-<span class='line'>1278</span>      *       [9.0, 0.5, 'black'],
-<span class='line'>1279</span>      *       [7.5, 1.0, 'purple'],
-<span class='line'>1280</span>      *   ]
-<span class='line'>1281</span>      *
-<span class='line'>1282</span>      *   let labelY = .4
-<span class='line'>1283</span>      *   runs.forEach((run,i) => {
-<span class='line'>1284</span>      *       board.create('segment',[[7,labelY-(i/50)],[9,labelY-(i/50)]],{strokeColor:run[2]})
-<span class='line'>1285</span>      *       board.create('text',[10,labelY-(i/50),`k=${run[0]}, &theta;=${run[1]}`])
-<span class='line'>1286</span>      *
-<span class='line'>1287</span>      *       // density
-<span class='line'>1288</span>      *       let x = Array(50000).fill(0).map(() => JXG.Math.Statistics.randomGamma(run[0],run[1]))
-<span class='line'>1289</span>      *       let res = JXG.Math.Statistics.histogram(x, { bins: 50, density: true, cumulative: false, range: [0, 20] });
-<span class='line'>1290</span>      *       board.create('curve', [res[1], res[0]], { strokeColor: run[2], strokeWidth:2});
-<span class='line'>1291</span>      *
-<span class='line'>1292</span>      *       // cumulative density
-<span class='line'>1293</span>      *       res = JXG.Math.Statistics.histogram(x, { bins: 50, density: true, cumulative: true, range: [0, 20] });
-<span class='line'>1294</span>      *       res[0].unshift(0)  // add zero to front so cumulative starts at zero
-<span class='line'>1295</span>      *       res[1].unshift(0)
-<span class='line'>1296</span>      *       board2.create('curve', [res[1], res[0]], { strokeColor: run[2], strokeWidth:2 });
-<span class='line'>1297</span>      *   })
-<span class='line'>1298</span>      *
-<span class='line'>1299</span>      *
-<span class='line'>1300</span>      * &lt;/pre>&lt;div id="JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302" class="jxgbox" style="width: 300px; height: 300px; float:left;">&lt;/div>
-<span class='line'>1301</span>      * &lt;div style='float:left;'>&nbsp;&nbsp;&lt;/div>
-<span class='line'>1302</span>      * &lt;div id="JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302-2" class="jxgbox" style="width: 300px; height: 300px;">&lt;/div>
-<span class='line'>1303</span>      * &lt;script type="text/javascript">
-<span class='line'>1304</span>      * {
-<span class='line'>1305</span>      *  let board = JXG.JSXGraph.initBoard('JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302',
-<span class='line'>1306</span>      *       { boundingbox: [-1.7, .5, 20, -.03], axis: true});
-<span class='line'>1307</span>      *  let board2 = JXG.JSXGraph.initBoard('JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302-2',
-<span class='line'>1308</span>      *       { boundingbox: [-1.6, 1.1, 20, -.06], axis: true});
-<span class='line'>1309</span>      *
-<span class='line'>1310</span>      *   let runs = [
-<span class='line'>1311</span>      *       [0.5, 1.0, 'brown'],
-<span class='line'>1312</span>      *       [1.0, 2.0, 'red'],
-<span class='line'>1313</span>      *       [2.0, 2.0, 'orange'],
-<span class='line'>1314</span>      *       [3.0, 2.0, 'yellow'],
-<span class='line'>1315</span>      *       [5.0, 1.0, 'green'],
-<span class='line'>1316</span>      *       [9.0, 0.5, 'black'],
-<span class='line'>1317</span>      *       [7.5, 1.0, 'purple'],
-<span class='line'>1318</span>      *   ]
-<span class='line'>1319</span>      *
-<span class='line'>1320</span>      *   let labelY = .4
-<span class='line'>1321</span>      *   runs.forEach((run,i) => {
-<span class='line'>1322</span>      *       board.create('segment',[[7,labelY-(i/50)],[9,labelY-(i/50)]],{strokeColor:run[2]})
-<span class='line'>1323</span>      *       board.create('text',[10,labelY-(i/50),`k=${run[0]}, &theta;=${run[1]}`])
-<span class='line'>1324</span>      *
-<span class='line'>1325</span>      *       let x = Array(50000).fill(0).map(() => JXG.Math.Statistics.randomGamma(run[0],run[1]))
-<span class='line'>1326</span>      *       let res = JXG.Math.Statistics.histogram(x, { bins: 50, density: true, cumulative: false, range: [0, 20] });
-<span class='line'>1327</span>      *       board.create('curve', [res[1], res[0]], { strokeColor: run[2], strokeWidth:2});
+<span class='line'>1208</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">good</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NAME">good</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1209</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">NaN</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1210</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1211</span> 
+<span class='line'>1212</span> </span><span class="WHIT">        </span><span class="COMM">// Naive method</span><span class="WHIT">
+<span class='line'>1213</span> </span><span class="WHIT">        </span><span class="COMM">// kk = Math.min(k, n - k);</span><span class="WHIT">
+<span class='line'>1214</span> </span><span class="WHIT">        </span><span class="COMM">// for (i = 0; i &lt; k; i ++) {</span><span class="WHIT">
+<span class='line'>1215</span> </span><span class="WHIT">        </span><span class="COMM">//     u = Math.random();</span><span class="WHIT">
+<span class='line'>1216</span> </span><span class="WHIT">        </span><span class="COMM">//     if (n * u &lt;= good) {</span><span class="WHIT">
+<span class='line'>1217</span> </span><span class="WHIT">        </span><span class="COMM">//         x += 1;</span><span class="WHIT">
+<span class='line'>1218</span> </span><span class="WHIT">        </span><span class="COMM">//         if (x === good) {</span><span class="WHIT">
+<span class='line'>1219</span> </span><span class="WHIT">        </span><span class="COMM">//             return x;</span><span class="WHIT">
+<span class='line'>1220</span> </span><span class="WHIT">        </span><span class="COMM">//         }</span><span class="WHIT">
+<span class='line'>1221</span> </span><span class="WHIT">        </span><span class="COMM">//         good -= 1;</span><span class="WHIT">
+<span class='line'>1222</span> </span><span class="WHIT">        </span><span class="COMM">//     }</span><span class="WHIT">
+<span class='line'>1223</span> </span><span class="WHIT">        </span><span class="COMM">//     n -= 1;</span><span class="WHIT">
+<span class='line'>1224</span> </span><span class="WHIT">        </span><span class="COMM">// }</span><span class="WHIT">
+<span class='line'>1225</span> </span><span class="WHIT">        </span><span class="COMM">// return x;</span><span class="WHIT">
+<span class='line'>1226</span> 
+<span class='line'>1227</span> </span><span class="WHIT">        </span><span class="COMM">// Implementation from</span><span class="WHIT">
+<span class='line'>1228</span> </span><span class="WHIT">        </span><span class="COMM">// Monte Carlo by George S. Fishman</span><span class="WHIT">
+<span class='line'>1229</span> </span><span class="WHIT">        </span><span class="COMM">// https://link.springer.com/book/10.1007/978-1-4757-2553-7</span><span class="WHIT">
+<span class='line'>1230</span> </span><span class="WHIT">        </span><span class="COMM">// page 218</span><span class="WHIT">
+<span class='line'>1231</span> </span><span class="WHIT">        </span><span class="COMM">//</span><span class="WHIT">
+<span class='line'>1232</span> </span><span class="WHIT">        </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">k</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1233</span> </span><span class="WHIT">        </span><span class="KEYW">while</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1234</span> </span><span class="WHIT">            </span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.random</span><span class="PUNC">(</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1235</span> </span><span class="WHIT">            </span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.floor</span><span class="PUNC">(</span><span class="NAME">u</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">y</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">d1</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">)</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1236</span> </span><span class="WHIT">            </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1237</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1238</span> </span><span class="WHIT">        </span><span class="NAME">x</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">d2</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">y</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1239</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">good</span><span class="WHIT"> </span><span class="PUNC">&lt;=</span><span class="WHIT"> </span><span class="NAME">bad</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1240</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1241</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1242</span> </span><span class="WHIT">            </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1243</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1244</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1245</span> 
+<span class='line'>1246</span> </span><span class="WHIT">    </span><span class="COMM">/**
+<span class='line'>1247</span>      * Compute the histogram of a dataset.
+<span class='line'>1248</span>      * Optional parameters can be supplied through a JavaScript object
+<span class='line'>1249</span>      * with the following default values:
+<span class='line'>1250</span>      * &lt;pre>
+<span class='line'>1251</span>      * {
+<span class='line'>1252</span>      *   bins: 10,          // Number of bins
+<span class='line'>1253</span>      *   range: false,      // false or array. The lower and upper range of the bins.
+<span class='line'>1254</span>      *                      // If not provided, range is simply [min(x), max(x)].
+<span class='line'>1255</span>      *                      // Values outside the range are ignored.
+<span class='line'>1256</span>      *   density: false,    // If true, normalize the counts by dividing by sum(counts)
+<span class='line'>1257</span>      *   cumulative: false
+<span class='line'>1258</span>      * }
+<span class='line'>1259</span>      * &lt;/pre>
+<span class='line'>1260</span>      * The function returns an array containing two arrays. The first array is of length bins+1
+<span class='line'>1261</span>      * containing the start values of the bins. The last entry contains the end values of the last bin.
+<span class='line'>1262</span>      * &lt;p>
+<span class='line'>1263</span>      * The second array contains the counts of each bin.
+<span class='line'>1264</span>      * @param {Array} x
+<span class='line'>1265</span>      * @param {Object} opt Optional parameters
+<span class='line'>1266</span>      * @returns Array [bin, counts] Array bins contains start values of bins, array counts contains
+<span class='line'>1267</span>      * the number of entries of x which are contained in each bin.
+<span class='line'>1268</span>      * @memberof JXG.Math.Statistics
+<span class='line'>1269</span>      *
+<span class='line'>1270</span>      * @example
+<span class='line'>1271</span>      *  let board = JXG.JSXGraph.initBoard('jxgbox',
+<span class='line'>1272</span>      *       { boundingbox: [-1.7, .5, 20, -.03], axis: true});
+<span class='line'>1273</span>      *  let board2 = JXG.JSXGraph.initBoard('jxgbox2',
+<span class='line'>1274</span>      *       { boundingbox: [-1.6, 1.1, 20, -.06], axis: true});
+<span class='line'>1275</span>      *
+<span class='line'>1276</span>      *   let runs = [
+<span class='line'>1277</span>      *       [0.5, 1.0, 'brown'],
+<span class='line'>1278</span>      *       [1.0, 2.0, 'red'],
+<span class='line'>1279</span>      *       [2.0, 2.0, 'orange'],
+<span class='line'>1280</span>      *       [3.0, 2.0, 'yellow'],
+<span class='line'>1281</span>      *       [5.0, 1.0, 'green'],
+<span class='line'>1282</span>      *       [9.0, 0.5, 'black'],
+<span class='line'>1283</span>      *       [7.5, 1.0, 'purple'],
+<span class='line'>1284</span>      *   ]
+<span class='line'>1285</span>      *
+<span class='line'>1286</span>      *   let labelY = .4
+<span class='line'>1287</span>      *   runs.forEach((run,i) => {
+<span class='line'>1288</span>      *       board.create('segment',[[7,labelY-(i/50)],[9,labelY-(i/50)]],{strokeColor:run[2]})
+<span class='line'>1289</span>      *       board.create('text',[10,labelY-(i/50),`k=${run[0]}, &theta;=${run[1]}`])
+<span class='line'>1290</span>      *
+<span class='line'>1291</span>      *       // density
+<span class='line'>1292</span>      *       let x = Array(50000).fill(0).map(() => JXG.Math.Statistics.randomGamma(run[0],run[1]))
+<span class='line'>1293</span>      *       let res = JXG.Math.Statistics.histogram(x, { bins: 50, density: true, cumulative: false, range: [0, 20] });
+<span class='line'>1294</span>      *       board.create('curve', [res[1], res[0]], { strokeColor: run[2], strokeWidth:2});
+<span class='line'>1295</span>      *
+<span class='line'>1296</span>      *       // cumulative density
+<span class='line'>1297</span>      *       res = JXG.Math.Statistics.histogram(x, { bins: 50, density: true, cumulative: true, range: [0, 20] });
+<span class='line'>1298</span>      *       res[0].unshift(0)  // add zero to front so cumulative starts at zero
+<span class='line'>1299</span>      *       res[1].unshift(0)
+<span class='line'>1300</span>      *       board2.create('curve', [res[1], res[0]], { strokeColor: run[2], strokeWidth:2 });
+<span class='line'>1301</span>      *   })
+<span class='line'>1302</span>      *
+<span class='line'>1303</span>      *
+<span class='line'>1304</span>      * &lt;/pre>&lt;div id="JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302" class="jxgbox" style="width: 300px; height: 300px; float:left;">&lt;/div>
+<span class='line'>1305</span>      * &lt;div style='float:left;'>&nbsp;&nbsp;&lt;/div>
+<span class='line'>1306</span>      * &lt;div id="JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302-2" class="jxgbox" style="width: 300px; height: 300px;">&lt;/div>
+<span class='line'>1307</span>      * &lt;script type="text/javascript">
+<span class='line'>1308</span>      * {
+<span class='line'>1309</span>      *  let board = JXG.JSXGraph.initBoard('JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302',
+<span class='line'>1310</span>      *       { boundingbox: [-1.7, .5, 20, -.03], axis: true});
+<span class='line'>1311</span>      *  let board2 = JXG.JSXGraph.initBoard('JXGda56df4d-a5a5-4c87-9ffc-9bbc1b512302-2',
+<span class='line'>1312</span>      *       { boundingbox: [-1.6, 1.1, 20, -.06], axis: true});
+<span class='line'>1313</span>      *
+<span class='line'>1314</span>      *   let runs = [
+<span class='line'>1315</span>      *       [0.5, 1.0, 'brown'],
+<span class='line'>1316</span>      *       [1.0, 2.0, 'red'],
+<span class='line'>1317</span>      *       [2.0, 2.0, 'orange'],
+<span class='line'>1318</span>      *       [3.0, 2.0, 'yellow'],
+<span class='line'>1319</span>      *       [5.0, 1.0, 'green'],
+<span class='line'>1320</span>      *       [9.0, 0.5, 'black'],
+<span class='line'>1321</span>      *       [7.5, 1.0, 'purple'],
+<span class='line'>1322</span>      *   ]
+<span class='line'>1323</span>      *
+<span class='line'>1324</span>      *   let labelY = .4
+<span class='line'>1325</span>      *   runs.forEach((run,i) => {
+<span class='line'>1326</span>      *       board.create('segment',[[7,labelY-(i/50)],[9,labelY-(i/50)]],{strokeColor:run[2]})
+<span class='line'>1327</span>      *       board.create('text',[10,labelY-(i/50),`k=${run[0]}, &theta;=${run[1]}`])
 <span class='line'>1328</span>      *
-<span class='line'>1329</span>      *       // cumulative density
-<span class='line'>1330</span>      *       res = JXG.Math.Statistics.histogram(x, { bins: 50, density: true, cumulative: true, range: [0, 20] });
-<span class='line'>1331</span>      *       res[0].unshift(0)  // add zero to front so cumulative starts at zero
-<span class='line'>1332</span>      *       res[1].unshift(0)
-<span class='line'>1333</span>      *       board2.create('curve', [res[1], res[0]], { strokeColor: run[2], strokeWidth:2 });
-<span class='line'>1334</span>      *   })
-<span class='line'>1335</span>      * }
-<span class='line'>1336</span>      * &lt;/script>&lt;pre>
-<span class='line'>1337</span>      *
-<span class='line'>1338</span>      */</span><span class="WHIT">
-<span class='line'>1339</span> </span><span class="WHIT">    </span><span class="NAME">histogram</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">x</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">opt</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1340</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">le</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">k</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1341</span> </span><span class="WHIT">            </span><span class="NAME">mi</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">ma</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">delta</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1342</span> </span><span class="WHIT">            </span><span class="NAME">range</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1343</span> </span><span class="WHIT">            </span><span class="NAME">s</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1344</span> </span><span class="WHIT">            </span><span class="NAME">counts</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="PUNC">[</span><span class="PUNC">]</span><span class="PUNC">,</span><span class="WHIT">
-<span class='line'>1345</span> </span><span class="WHIT">            </span><span class="NAME">bins</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="PUNC">[</span><span class="PUNC">]</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1346</span> 
-<span class='line'>1347</span> </span><span class="WHIT">        </span><span class="COMM">// Evaluate number of bins</span><span class="WHIT">
-<span class='line'>1348</span> </span><span class="WHIT">        </span><span class="NAME">num_bins</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">opt.bins</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NUMB">10</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1349</span> 
-<span class='line'>1350</span> </span><span class="WHIT">        </span><span class="COMM">// Evaluate range</span><span class="WHIT">
-<span class='line'>1351</span> </span><span class="WHIT">        </span><span class="NAME">range</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">opt.range</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="KEYW">false</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1352</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">range</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="KEYW">false</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1353</span> </span><span class="WHIT">            </span><span class="NAME">mi</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.min.apply</span><span class="PUNC">(</span><span class="KEYW">null</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1354</span> </span><span class="WHIT">            </span><span class="NAME">ma</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.max.apply</span><span class="PUNC">(</span><span class="KEYW">null</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1355</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1356</span> </span><span class="WHIT">            </span><span class="NAME">mi</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">range</span><span class="PUNC">[</span><span class="NUMB">0</span><span class="PUNC">]</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1357</span> </span><span class="WHIT">            </span><span class="NAME">ma</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">range</span><span class="PUNC">[</span><span class="NUMB">1</span><span class="PUNC">]</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1358</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1359</span> 
-<span class='line'>1360</span> </span><span class="WHIT">        </span><span class="COMM">// Set uniform delta</span><span class="WHIT">
-<span class='line'>1361</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">num_bins</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1362</span> </span><span class="WHIT">            </span><span class="NAME">delta</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">ma</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">mi</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">num_bins</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1363</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1364</span> </span><span class="WHIT">            </span><span class="NAME">delta</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1365</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1366</span> 
-<span class='line'>1367</span> </span><span class="WHIT">        </span><span class="COMM">// Set the bins and init the counts array</span><span class="WHIT">
-<span class='line'>1368</span> </span><span class="WHIT">        </span><span class="KEYW">for</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">++</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1369</span> </span><span class="WHIT">            </span><span class="NAME">counts.push</span><span class="PUNC">(</span><span class="NUMB">0</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1370</span> </span><span class="WHIT">            </span><span class="NAME">bins.push</span><span class="PUNC">(</span><span class="NAME">mi</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NAME">delta</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1371</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1372</span> </span><span class="WHIT">        </span><span class="COMM">// bins.push(ma);</span><span class="WHIT">
-<span class='line'>1373</span> 
-<span class='line'>1374</span> </span><span class="WHIT">        </span><span class="COMM">// Determine the counts</span><span class="WHIT">
-<span class='line'>1375</span> </span><span class="WHIT">        </span><span class="NAME">le</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">x.length</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1376</span> </span><span class="WHIT">        </span><span class="KEYW">for</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">le</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">++</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1377</span> </span><span class="WHIT">            </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.floor</span><span class="PUNC">(</span><span class="PUNC">(</span><span class="NAME">x</span><span class="PUNC">[</span><span class="NAME">i</span><span class="PUNC">]</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">mi</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">delta</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1378</span> </span><span class="WHIT">            </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">>=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="WHIT"> </span><span class="PUNC">&&</span><span class="WHIT"> </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1379</span> </span><span class="WHIT">                </span><span class="NAME">counts</span><span class="PUNC">[</span><span class="NAME">k</span><span class="PUNC">]</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1380</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1381</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1382</span> 
-<span class='line'>1383</span> </span><span class="WHIT">        </span><span class="COMM">// Normalize if density===true</span><span class="WHIT">
-<span class='line'>1384</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">opt.density</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1385</span> </span><span class="WHIT">            </span><span class="NAME">s</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">JXG.Math.Statistics.sum</span><span class="PUNC">(</span><span class="NAME">counts</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1386</span> </span><span class="WHIT">            </span><span class="KEYW">for</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">++</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1387</span> </span><span class="WHIT">                </span><span class="NAME">counts</span><span class="PUNC">[</span><span class="NAME">i</span><span class="PUNC">]</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">s</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NAME">delta</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1388</span> </span><span class="WHIT">                </span><span class="COMM">// counts[i] /= s;</span><span class="WHIT">
-<span class='line'>1389</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1390</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1391</span> 
-<span class='line'>1392</span> </span><span class="WHIT">        </span><span class="COMM">// Cumulative counts</span><span class="WHIT">
-<span class='line'>1393</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">opt.cumulative</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1394</span> </span><span class="WHIT">            </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">opt.density</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1395</span> </span><span class="WHIT">                </span><span class="KEYW">for</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">++</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1396</span> </span><span class="WHIT">                    </span><span class="NAME">counts</span><span class="PUNC">[</span><span class="NAME">i</span><span class="PUNC">]</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">delta</span><span class="PUNC">;</span><span class="WHIT">  </span><span class="COMM">// normalize</span><span class="WHIT">
-<span class='line'>1397</span> </span><span class="WHIT">                </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1398</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">for</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">++</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
-<span class='line'>1399</span> </span><span class="WHIT">                </span><span class="NAME">counts</span><span class="PUNC">[</span><span class="NAME">i</span><span class="PUNC">]</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">counts</span><span class="PUNC">[</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">]</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1400</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1401</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1402</span> 
-<span class='line'>1403</span> </span><span class="WHIT">        </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="PUNC">[</span><span class="NAME">counts</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">bins</span><span class="PUNC">]</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1404</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="WHIT">
-<span class='line'>1405</span> </span><span class="PUNC">}</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1406</span> 
-<span class='line'>1407</span> </span><span class="NAME">export</span><span class="WHIT"> </span><span class="KEYW">default</span><span class="WHIT"> </span><span class="NAME">Mat.Statistics</span><span class="PUNC">;</span><span class="WHIT">
-<span class='line'>1408</span> </span></pre></body></html>
+<span class='line'>1329</span>      *       let x = Array(50000).fill(0).map(() => JXG.Math.Statistics.randomGamma(run[0],run[1]))
+<span class='line'>1330</span>      *       let res = JXG.Math.Statistics.histogram(x, { bins: 50, density: true, cumulative: false, range: [0, 20] });
+<span class='line'>1331</span>      *       board.create('curve', [res[1], res[0]], { strokeColor: run[2], strokeWidth:2});
+<span class='line'>1332</span>      *
+<span class='line'>1333</span>      *       // cumulative density
+<span class='line'>1334</span>      *       res = JXG.Math.Statistics.histogram(x, { bins: 50, density: true, cumulative: true, range: [0, 20] });
+<span class='line'>1335</span>      *       res[0].unshift(0)  // add zero to front so cumulative starts at zero
+<span class='line'>1336</span>      *       res[1].unshift(0)
+<span class='line'>1337</span>      *       board2.create('curve', [res[1], res[0]], { strokeColor: run[2], strokeWidth:2 });
+<span class='line'>1338</span>      *   })
+<span class='line'>1339</span>      * }
+<span class='line'>1340</span>      * &lt;/script>&lt;pre>
+<span class='line'>1341</span>      *
+<span class='line'>1342</span>      */</span><span class="WHIT">
+<span class='line'>1343</span> </span><span class="WHIT">    </span><span class="NAME">histogram</span><span class="PUNC">:</span><span class="WHIT"> </span><span class="KEYW">function</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">x</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">opt</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1344</span> </span><span class="WHIT">        </span><span class="KEYW">var</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">le</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">k</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1345</span> </span><span class="WHIT">            </span><span class="NAME">mi</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">ma</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">delta</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1346</span> </span><span class="WHIT">            </span><span class="NAME">range</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1347</span> </span><span class="WHIT">            </span><span class="NAME">s</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1348</span> </span><span class="WHIT">            </span><span class="NAME">counts</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="PUNC">[</span><span class="PUNC">]</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1349</span> </span><span class="WHIT">            </span><span class="NAME">bins</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="PUNC">[</span><span class="PUNC">]</span><span class="PUNC">,</span><span class="WHIT">
+<span class='line'>1350</span> </span><span class="WHIT">            </span><span class="NAME">no_bin</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT">        </span><span class="COMM">// Count of long tail elements not in histogram range</span><span class="WHIT">
+<span class='line'>1351</span> 
+<span class='line'>1352</span> </span><span class="WHIT">        </span><span class="COMM">// Evaluate number of bins</span><span class="WHIT">
+<span class='line'>1353</span> </span><span class="WHIT">        </span><span class="NAME">num_bins</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">opt.bins</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="NUMB">10</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1354</span> 
+<span class='line'>1355</span> </span><span class="WHIT">        </span><span class="COMM">// Evaluate range</span><span class="WHIT">
+<span class='line'>1356</span> </span><span class="WHIT">        </span><span class="NAME">range</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">opt.range</span><span class="WHIT"> </span><span class="PUNC">||</span><span class="WHIT"> </span><span class="KEYW">false</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1357</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">range</span><span class="WHIT"> </span><span class="PUNC">===</span><span class="WHIT"> </span><span class="KEYW">false</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1358</span> </span><span class="WHIT">            </span><span class="NAME">mi</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.min.apply</span><span class="PUNC">(</span><span class="KEYW">null</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1359</span> </span><span class="WHIT">            </span><span class="NAME">ma</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.max.apply</span><span class="PUNC">(</span><span class="KEYW">null</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">x</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1360</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1361</span> </span><span class="WHIT">            </span><span class="NAME">mi</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">range</span><span class="PUNC">[</span><span class="NUMB">0</span><span class="PUNC">]</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1362</span> </span><span class="WHIT">            </span><span class="NAME">ma</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">range</span><span class="PUNC">[</span><span class="NUMB">1</span><span class="PUNC">]</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1363</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1364</span> 
+<span class='line'>1365</span> </span><span class="WHIT">        </span><span class="COMM">// Set uniform delta</span><span class="WHIT">
+<span class='line'>1366</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">num_bins</span><span class="WHIT"> </span><span class="PUNC">></span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1367</span> </span><span class="WHIT">            </span><span class="NAME">delta</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">ma</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">mi</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">num_bins</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1368</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1369</span> </span><span class="WHIT">            </span><span class="NAME">delta</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1370</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1371</span> 
+<span class='line'>1372</span> </span><span class="WHIT">        </span><span class="COMM">// Set the bins and init the counts array</span><span class="WHIT">
+<span class='line'>1373</span> </span><span class="WHIT">        </span><span class="KEYW">for</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">++</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1374</span> </span><span class="WHIT">            </span><span class="NAME">counts.push</span><span class="PUNC">(</span><span class="NUMB">0</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1375</span> </span><span class="WHIT">            </span><span class="NAME">bins.push</span><span class="PUNC">(</span><span class="NAME">mi</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NAME">delta</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1376</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1377</span> </span><span class="WHIT">        </span><span class="COMM">// bins.push(ma);</span><span class="WHIT">
+<span class='line'>1378</span> 
+<span class='line'>1379</span> </span><span class="WHIT">        </span><span class="COMM">// Determine the counts</span><span class="WHIT">
+<span class='line'>1380</span> </span><span class="WHIT">        </span><span class="NAME">le</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">x.length</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1381</span> </span><span class="WHIT">        </span><span class="KEYW">for</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">le</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">++</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1382</span> </span><span class="WHIT">            </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">Math.floor</span><span class="PUNC">(</span><span class="PUNC">(</span><span class="NAME">x</span><span class="PUNC">[</span><span class="NAME">i</span><span class="PUNC">]</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NAME">mi</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="WHIT"> </span><span class="NAME">delta</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1383</span> </span><span class="WHIT">            </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">>=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="WHIT"> </span><span class="PUNC">&&</span><span class="WHIT"> </span><span class="NAME">k</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1384</span> </span><span class="WHIT">                </span><span class="NAME">counts</span><span class="PUNC">[</span><span class="NAME">k</span><span class="PUNC">]</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1385</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">else</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1386</span> </span><span class="WHIT">                </span><span class="NAME">no_bin</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1387</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1388</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1389</span> 
+<span class='line'>1390</span> </span><span class="WHIT">        </span><span class="COMM">// Normalize if density===true</span><span class="WHIT">
+<span class='line'>1391</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">opt.density</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1392</span> </span><span class="WHIT">            </span><span class="NAME">s</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">JXG.Math.Statistics.sum</span><span class="PUNC">(</span><span class="NAME">counts</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="WHIT"> </span><span class="NAME">no_bin</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="COMM">// Normalize including long tail</span><span class="WHIT">
+<span class='line'>1393</span> </span><span class="WHIT">            </span><span class="KEYW">for</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">++</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1394</span> </span><span class="WHIT">                </span><span class="NAME">counts</span><span class="PUNC">[</span><span class="NAME">i</span><span class="PUNC">]</span><span class="WHIT"> </span><span class="PUNC">/</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">s</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="WHIT"> </span><span class="NAME">delta</span><span class="PUNC">)</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1395</span> </span><span class="WHIT">                </span><span class="COMM">// counts[i] /= s;</span><span class="WHIT">
+<span class='line'>1396</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1397</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1398</span> 
+<span class='line'>1399</span> </span><span class="WHIT">        </span><span class="COMM">// Cumulative counts</span><span class="WHIT">
+<span class='line'>1400</span> </span><span class="WHIT">        </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">opt.cumulative</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1401</span> </span><span class="WHIT">            </span><span class="KEYW">if</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">opt.density</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1402</span> </span><span class="WHIT">                </span><span class="KEYW">for</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">0</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">++</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1403</span> </span><span class="WHIT">                    </span><span class="NAME">counts</span><span class="PUNC">[</span><span class="NAME">i</span><span class="PUNC">]</span><span class="WHIT"> </span><span class="PUNC">*</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">delta</span><span class="PUNC">;</span><span class="WHIT">  </span><span class="COMM">// Normalize</span><span class="WHIT">
+<span class='line'>1404</span> </span><span class="WHIT">                </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1405</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT"> </span><span class="KEYW">for</span><span class="WHIT"> </span><span class="PUNC">(</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">&lt;</span><span class="WHIT"> </span><span class="NAME">num_bins</span><span class="PUNC">;</span><span class="WHIT"> </span><span class="NAME">i</span><span class="PUNC">++</span><span class="PUNC">)</span><span class="WHIT"> </span><span class="PUNC">{</span><span class="WHIT">
+<span class='line'>1406</span> </span><span class="WHIT">                </span><span class="NAME">counts</span><span class="PUNC">[</span><span class="NAME">i</span><span class="PUNC">]</span><span class="WHIT"> </span><span class="PUNC">+</span><span class="PUNC">=</span><span class="WHIT"> </span><span class="NAME">counts</span><span class="PUNC">[</span><span class="NAME">i</span><span class="WHIT"> </span><span class="PUNC">-</span><span class="WHIT"> </span><span class="NUMB">1</span><span class="PUNC">]</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1407</span> </span><span class="WHIT">            </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1408</span> </span><span class="WHIT">        </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1409</span> 
+<span class='line'>1410</span> </span><span class="WHIT">        </span><span class="KEYW">return</span><span class="WHIT"> </span><span class="PUNC">[</span><span class="NAME">counts</span><span class="PUNC">,</span><span class="WHIT"> </span><span class="NAME">bins</span><span class="PUNC">]</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1411</span> </span><span class="WHIT">    </span><span class="PUNC">}</span><span class="WHIT">
+<span class='line'>1412</span> </span><span class="PUNC">}</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1413</span> 
+<span class='line'>1414</span> </span><span class="NAME">export</span><span class="WHIT"> </span><span class="KEYW">default</span><span class="WHIT"> </span><span class="NAME">Mat.Statistics</span><span class="PUNC">;</span><span class="WHIT">
+<span class='line'>1415</span> </span></pre></body></html>