Skip to content

Instantly share code, notes, and snippets.

@d--b

d--b/brownian_motion_simulation.sbd

Last active August 29, 2015 14:04
Show Gist options
  • Save d--b/125260dcdc636f0191ac to your computer and use it in GitHub Desktop.
Save d--b/125260dcdc636f0191ac to your computer and use it in GitHub Desktop.
Download ZIP
Updated through Spreadboard
Raw
brownian_motion_simulation.sbd
{
"tablenames": {
"0": "table_0",
"1": "table_1",
"2": "table_2"
},
"calctree": {
"0": [
["=function () {\n hist = Array.apply(null, Array(1000)).map(function (x) {\n return 0;\n });\n setInterval(function () {\n spreadboard.SetFormulaNoRecalc(0, 0, 1, '=Math.random()');\n spreadboard.CalculateAll();\n }, 1);\n return 'init done.';\n}()", "=Math.random()", "=function NormSInv(p) {\n var a1 = -39.6968302866538, a2 = 220.946098424521, a3 = -275.928510446969;\n var a4 = 138.357751867269, a5 = -30.6647980661472, a6 = 2.50662827745924;\n var b1 = -54.4760987982241, b2 = 161.585836858041, b3 = -155.698979859887;\n var b4 = 66.8013118877197, b5 = -13.2806815528857, c1 = -0.00778489400243029;\n var c2 = -0.322396458041136, c3 = -2.40075827716184, c4 = -2.54973253934373;\n var c5 = 4.37466414146497, c6 = 2.93816398269878, d1 = 0.00778469570904146;\n var d2 = 0.32246712907004, d3 = 2.445134137143, d4 = 3.75440866190742;\n var p_low = 0.02425, p_high = 1 - p_low;\n var q, r;\n var retVal;\n if (p < 0 || p > 1) {\n alert('NormSInv: Argument out of range.');\n retVal = 0;\n } else if (p < p_low) {\n q = Math.sqrt(-2 * Math.log(p));\n retVal = (((((c1 * q + c2) * q + c3) * q + c4) * q + c5) * q + c6) / ((((d1 * q + d2) * q + d3) * q + d4) * q + 1);\n } else if (p <= p_high) {\n q = p - 0.5;\n r = q * q;\n retVal = (((((a1 * r + a2) * r + a3) * r + a4) * r + a5) * r + a6) * q / (((((b1 * r + b2) * r + b3) * r + b4) * r + b5) * r + 1);\n } else {\n q = Math.sqrt(-2 * Math.log(1 - p));\n retVal = -(((((c1 * q + c2) * q + c3) * q + c4) * q + c5) * q + c6) / ((((d1 * q + d2) * q + d3) * q + d4) * q + 1);\n }\n return retVal;\n}", "", ""],
["hist", "=function () {\n hist.splice(0, 1);\n hist.push(C1(B1) + hist[hist.length - 1]);\n return hist;\n}()", "=Chart.sparkline(B2)", "", ""],
["", "", "", "", ""],
["last", "=hist[hist.length - 1] + 0 * $B$1", "", "", ""],
["", "", "", "", ""],
["mvavg", "=function (arr, per) {\n var m = 0;\n for (var i = 0; i < per; i++) {\n m = m + arr[i] / per;\n }\n var result = new Array(arr.length);\n for (var i = 0; i < per; i++) {\n result[i] = m;\n }\n for (var i = 1; i < arr.length - per + 1; i++) {\n result[i + per - 1] = result[i + per - 2] + (-1 * arr[i - 1] + arr[i + per - 1]) / per;\n }\n return result;\n}", "", "", ""],
[100, "=B6(B2, A7)", "=Chart.sparkline(B7)", "", ""],
[300, "=B6(B2, A8)", "=Chart.sparkline(B8)", "", ""],
["", "", "", "", ""],
["deriv", "=function (arr, n) {\n var result = new Array(arr.length - n);\n for (var i = 0; i < arr.length - n; i++) {\n result[i] = arr[i + n] - arr[i];\n }\n return result;\n}", "", "", ""],
[40, "=B10(B2, A11)", "=Chart.sparkline(B11)", "", ""],
["", "=B10(B8, A11)", "=Chart.sparkline(B12)", "", ""],
["", "", "", "", ""],
["", "", "", "", ""],
["", "", "", "", ""],
["", "", "", "", ""],
["", "", "", "", ""]
],
"1": [
["=Chart.sparkline([\n table_0.$B$2,\n table_0.B7,\n table_0.B8\n])"]
],
"2": [
["=Chart.sparkline(table_0.B11_12)"]
]
},
"view": {
"TableDetails": {
"0": {
"left": 120,
"top": 64,
"nrows": 17,
"ncols": 5,
"colwidths": [63, 63, 63, 63, 63],
"rowheights": [15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15]
},
"1": {
"left": 120,
"top": 360,
"nrows": 1,
"ncols": 1,
"colwidths": [319],
"rowheights": [87]
},
"2": {
"left": 120,
"top": 456,
"nrows": 1,
"ncols": 1,
"colwidths": [319],
"rowheights": [79]
}
}
}
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment