d9a8899b-48d6-4974-ac9e-9b4cfb2c7329 9 6 5 fraction

&[DATE] &[TIME] - &[FILENAME]

This algorithm is sourced from the Mathsoft Collaboratory. The programmed version reproduces the original from author Mircea Andreacut [Roumanie]... from recollection. Plots fast Mathcad 11, almost impossible Smath ?

V n place n  decimals in V i 1 V rows range r i el V i el n round : for r 31 line :

t n Φ 1 t n - sign + 2 / :

M Min Max scale Max Min - (M M min - M max M min - / * Min + :

n 200 :

k 4 :

k 4, 8, 16 ... ≡

Leave 'n' minimal, otherwise Smath issues: Exception of Type "SystemOutofMemoryException" was thrown

X N k n * : u 1 el 1 Random 1 el : v 1 el 1 Random 1 el : i 1 N range u i 1 + el 4 u i el * 1 u i el - (* : v i 1 + el 4 v i el * 1 v i el - (* : b i el u i el v i el Φ : 31 line for m 1 N 1 - k / 1 - range w m el b m k * j + el 2 j^* (j 0 k sum eval : for w 1 - 1 scale 61 line :

Y N k n * : u 1 el 1 Random 1 el : v 1 el 1 Random 1 el : i 1 N range u i 1 + el 4 u i el * 1 u i el - (* : v i 1 + el 4 v i el * 1 v i el - (* : b i el u i el v i el Φ : 31 line for m 1 N 1 - k / 1 - range w m el b m k * j + el 2 j^* (j 0 k sum eval : for w 1 - 1 scale 61 line :

t0 1 time :

X Y augment

1 time t0 - 39 sec *

Given as is: I don't see better distribution. Curiosity for [n=500, k=4] Timing for [n=500, k=8] => 700 sec.