diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/README.md b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/README.md new file mode 100644 index 000000000000..fdb506dbbde6 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/README.md @@ -0,0 +1,138 @@ + + +# Median + +> [Bradford][bradford-distribution] distribution [median][median]. + + + +
+ +The [median][median] for a [Bradford][bradford-distribution] random variable is + + + +```math +\mathop{\mathrm{Median}}\left( X \right) = \frac{\sqrt{1 + c} - 1}{c} +``` + + + + + +where `c` is the shape parameter. + +
+ + + + + +
+ +## Usage + +```javascript +var median = require( '@stdlib/stats/base/dists/bradford/median' ); +``` + +#### median( c ) + +Returns the [median][median] of a [Bradford][bradford-distribution] distribution with shape parameter `c`. + +```javascript +var v = median( 0.1 ); +// returns ~0.488 + +v = median( 10.0 ); +// returns ~0.232 +``` + +If provided a shape parameter `c <= 0`, the function returns `NaN`. + +```javascript +var v = median( 0.0 ); +// returns NaN + +v = median( -1.5 ); +// returns NaN +``` + +
+ + + + + +
+ +
+ + + + + +
+ +## Examples + + + +```javascript +var uniform = require( '@stdlib/random/array/uniform' ); +var median = require( '@stdlib/stats/base/dists/bradford/median' ); + +var c = uniform( 10, 0.1, 10.0 ); + +var v; +var i; +for ( i = 0; i < c.length; i++ ) { + v = median( c[ i ] ); + console.log( 'c: %d, Median(X;c): %d', c[ i ].toFixed( 4 ), v.toFixed( 4 ) ); +} +``` + +
+ + + + + +
+ +
+ + + + + +
+ +[bradford-distribution]: https://en.wikipedia.org/wiki/Bradford%27s_law + +[median]: https://en.wikipedia.org/wiki/Median + +
+ + diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/benchmark/benchmark.js new file mode 100644 index 000000000000..13f37e7676d8 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/benchmark/benchmark.js @@ -0,0 +1,52 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2025 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pkg = require( './../package.json' ).name; +var median = require( './../lib' ); + + +// MAIN // + +bench( pkg, function benchmark( b ) { + var c; + var y; + var i; + + c = uniform( 100, 0.1, 10.0 ); + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = median( c[ i % c.length ] ); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/img/equation_bradford_median.svg b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/img/equation_bradford_median.svg new file mode 100644 index 000000000000..1e1169655201 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/img/equation_bradford_median.svg @@ -0,0 +1,51 @@ + +normal upper M times normal e times normal d times normal i times normal a times normal n left-parenthesis upper X right-parenthesis equals StartFraction StartRoot 1 plus c EndRoot minus 1 Over c EndFraction + + + \ No newline at end of file diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/repl.txt b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/repl.txt new file mode 100644 index 000000000000..9626c52886aa --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/repl.txt @@ -0,0 +1,28 @@ + +{{alias}}( c ) + Returns the median of a Bradford distribution with shape parameter `c`. + + If `c <= 0`, the function returns `NaN`. + + Parameters + ---------- + c: number + Shape parameter. + + Returns + ------- + out: number + Median. + + Examples + -------- + > var v = {{alias}}( 0.1 ) + ~0.488 + > v = {{alias}}( 0.5 ) + ~0.449 + > v = {{alias}}( 10.0 ) + ~0.232 + + See Also + -------- + diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/types/index.d.ts b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/types/index.d.ts new file mode 100644 index 000000000000..da4fe843eb82 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/types/index.d.ts @@ -0,0 +1,60 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2025 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +// TypeScript Version: 4.1 + +/** +* Returns the median of a Bradford distribution. +* +* ## Notes +* +* - If `c <= 0`, the function returns `NaN`. +* +* @param c - shape parameter +* @returns median +* +* @example +* var v = median( 0.1 ); +* // returns ~0.488 +* +* @example +* var v = median( 0.5 ); +* // returns ~0.449 +* +* @example +* var v = median( 10.0 ); +* // returns ~0.232 +* +* @example +* var v = median( 0.0 ); +* // returns NaN +* +* @example +* var v = median( -1.0 ); +* // returns NaN +* +* @example +* var v = median( NaN ); +* // returns NaN +*/ +declare function median( c: number ): number; + + +// EXPORTS // + +export = median; diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/types/test.ts b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/types/test.ts new file mode 100644 index 000000000000..bd0c67d3e16c --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/docs/types/test.ts @@ -0,0 +1,44 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2025 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +import median = require( './index' ); + + +// TESTS // + +// The function returns a number... +{ + median( 0.3 ); // $ExpectType number +} + +// The compiler throws an error if the function is provided a value other than a number... +{ + median( true ); // $ExpectError + median( false ); // $ExpectError + median( null ); // $ExpectError + median( undefined ); // $ExpectError + median( '5' ); // $ExpectError + median( [] ); // $ExpectError + median( {} ); // $ExpectError + median( ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the function is provided insufficient arguments... +{ + median(); // $ExpectError +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/examples/index.js b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/examples/index.js new file mode 100644 index 000000000000..9e968a733db1 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/examples/index.js @@ -0,0 +1,31 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2025 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +var uniform = require( '@stdlib/random/array/uniform' ); +var median = require( './../lib' ); + +var c = uniform( 10, 0.1, 10.0 ); + +var v; +var i; +for ( i = 0; i < c.length; i++ ) { + v = median( c[ i ] ); + console.log( 'c: %d, Median(X;c): %d', c[ i ].toFixed( 4 ), v.toFixed( 4 ) ); +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/lib/index.js b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/lib/index.js new file mode 100644 index 000000000000..fd7148b99859 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/lib/index.js @@ -0,0 +1,43 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2025 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +/** +* Bradford distribution median. +* +* @module @stdlib/stats/base/dists/bradford/median +* +* @example +* var median = require( '@stdlib/stats/base/dists/bradford/median' ); +* +* var v = median( 0.1 ); +* // returns ~0.488 +* +* v = median( 0.5 ); +* // returns ~0.449 +*/ + +// MODULES // + +var median = require( './main.js' ); + + +// EXPORTS // + +module.exports = median; diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/lib/main.js b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/lib/main.js new file mode 100644 index 000000000000..fbb42e0f8505 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/lib/main.js @@ -0,0 +1,72 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2025 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var sqrt = require( '@stdlib/math/base/special/sqrt' ); + + +// MAIN // + +/** +* Returns the median of a Bradford distribution. +* +* @param {PositiveNumber} c - shape parameter +* @returns {PositiveNumber} median +* +* @example +* var v = median( 0.1 ); +* // returns ~0.488 +* +* @example +* var v = median( 0.5 ); +* // returns ~0.449 +* +* @example +* var v = median( 10.0 ); +* // returns ~0.232 +* +* @example +* var v = median( 0.0 ); +* // returns NaN +* +* @example +* var v = median( -1.0 ); +* // returns NaN +* +* @example +* var v = median( NaN ); +* // returns NaN +*/ +function median( c ) { + if ( + isnan( c ) || + c <= 0.0 + ) { + return NaN; + } + return ( sqrt( 1.0 + c ) - 1.0 ) / c; +} + + +// EXPORTS // + +module.exports = median; diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/package.json b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/package.json new file mode 100644 index 000000000000..0928a56fd902 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/package.json @@ -0,0 +1,69 @@ +{ + "name": "@stdlib/stats/base/dists/bradford/median", + "version": "0.0.0", + "description": "Bradford distribution median.", + "license": "Apache-2.0", + "author": { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + }, + "contributors": [ + { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + } + ], + "main": "./lib", + "directories": { + "benchmark": "./benchmark", + "doc": "./docs", + "example": "./examples", + "lib": "./lib", + "test": "./test" + }, + "types": "./docs/types", + "scripts": {}, + "homepage": "https://github.com/stdlib-js/stdlib", + "repository": { + "type": "git", + "url": "git://github.com/stdlib-js/stdlib.git" + }, + "bugs": { + "url": "https://github.com/stdlib-js/stdlib/issues" + }, + "dependencies": {}, + "devDependencies": {}, + "engines": { + "node": ">=0.10.0", + "npm": ">2.7.0" + }, + "os": [ + "aix", + "darwin", + "freebsd", + "linux", + "macos", + "openbsd", + "sunos", + "win32", + "windows" + ], + "keywords": [ + "stdlib", + "stdmath", + "statistics", + "stats", + "distribution", + "dist", + "bradford", + "bradford-distribution", + "parameter", + "shape-parameter", + "continuous", + "skewed", + "median", + "quantile", + "location", + "percentile" + ] +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/test/fixtures/python/data.json b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/test/fixtures/python/data.json new file mode 100644 index 000000000000..0e98ebcf9f0f --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/test/fixtures/python/data.json @@ -0,0 +1 @@ +{"x": [0.1, 0.10990990990990991, 0.11981981981981983, 0.12972972972972974, 0.13963963963963966, 0.14954954954954958, 0.15945945945945947, 0.1693693693693694, 0.17927927927927928, 0.1891891891891892, 0.19909909909909912, 0.20900900900900904, 0.21891891891891893, 0.22882882882882885, 0.23873873873873877, 0.24864864864864866, 0.2585585585585586, 0.2684684684684685, 0.2783783783783784, 0.28828828828828834, 0.29819819819819826, 0.3081081081081081, 0.31801801801801804, 0.32792792792792796, 0.33783783783783783, 0.34774774774774775, 0.35765765765765767, 0.3675675675675676, 0.3774774774774775, 0.3873873873873874, 0.39729729729729735, 0.40720720720720727, 0.4171171171171172, 0.4270270270270271, 0.436936936936937, 0.44684684684684695, 0.45675675675675675, 0.4666666666666667, 0.4765765765765766, 0.4864864864864865, 0.49639639639639643, 0.5063063063063064, 0.5162162162162163, 0.5261261261261262, 0.5360360360360361, 0.545945945945946, 0.555855855855856, 0.5657657657657659, 0.5756756756756757, 0.5855855855855856, 0.5954954954954955, 0.6054054054054054, 0.6153153153153154, 0.6252252252252253, 0.6351351351351352, 0.6450450450450451, 0.654954954954955, 0.664864864864865, 0.6747747747747748, 0.6846846846846847, 0.6945945945945946, 0.7045045045045045, 0.7144144144144144, 0.7243243243243244, 0.7342342342342343, 0.7441441441441442, 0.7540540540540541, 0.763963963963964, 0.773873873873874, 0.7837837837837839, 0.7936936936936938, 0.8036036036036036, 0.8135135135135135, 0.8234234234234235, 0.8333333333333334, 0.8432432432432433, 0.8531531531531532, 0.8630630630630631, 0.872972972972973, 0.882882882882883, 0.8927927927927929, 0.9027027027027028, 0.9126126126126127, 0.9225225225225226, 0.9324324324324325, 0.9423423423423424, 0.9522522522522523, 0.9621621621621622, 0.9720720720720721, 0.981981981981982, 0.991891891891892, 1.001801801801802, 1.0117117117117118, 1.0216216216216218, 1.0315315315315317, 1.0414414414414417, 1.0513513513513515, 1.0612612612612613, 1.0711711711711713, 1.0810810810810811, 1.0909909909909912, 1.1009009009009012, 1.110810810810811, 1.1207207207207208, 1.1306306306306309, 1.1405405405405407, 1.1504504504504507, 1.1603603603603605, 1.1702702702702705, 1.1801801801801803, 1.1900900900900904, 1.2000000000000002, 1.2099099099099102, 1.21981981981982, 1.22972972972973, 1.2396396396396399, 1.2495495495495497, 1.2594594594594597, 1.2693693693693695, 1.2792792792792795, 1.2891891891891893, 1.2990990990990994, 1.3090090090090092, 1.3189189189189192, 1.328828828828829, 1.338738738738739, 1.3486486486486489, 1.358558558558559, 1.3684684684684687, 1.3783783783783785, 1.3882882882882885, 1.3981981981981983, 1.4081081081081084, 1.4180180180180182, 1.4279279279279282, 1.437837837837838, 1.447747747747748, 1.4576576576576579, 1.467567567567568, 1.4774774774774777, 1.4873873873873877, 1.4972972972972975, 1.5072072072072074, 1.5171171171171174, 1.5270270270270272, 1.5369369369369372, 1.546846846846847, 1.556756756756757, 1.5666666666666669, 1.576576576576577, 1.5864864864864867, 1.5963963963963967, 1.6063063063063066, 1.6162162162162166, 1.6261261261261264, 1.6360360360360362, 1.6459459459459462, 1.655855855855856, 1.665765765765766, 1.6756756756756759, 1.685585585585586, 1.6954954954954957, 1.7054054054054057, 1.7153153153153156, 1.7252252252252256, 1.7351351351351354, 1.7450450450450454, 1.7549549549549552, 1.764864864864865, 1.774774774774775, 1.7846846846846849, 1.794594594594595, 1.8045045045045047, 1.8144144144144148, 1.8243243243243246, 1.8342342342342346, 1.8441441441441444, 1.8540540540540544, 1.8639639639639642, 1.8738738738738743, 1.883783783783784, 1.893693693693694, 1.903603603603604, 1.9135135135135137, 1.9234234234234238, 1.9333333333333336, 1.9432432432432436, 1.9531531531531534, 1.9630630630630634, 1.9729729729729732, 1.9828828828828833, 1.992792792792793, 2.002702702702703, 2.012612612612613, 2.0225225225225225, 2.0324324324324325, 2.0423423423423426, 2.0522522522522526, 2.062162162162162, 2.0720720720720722, 2.0819819819819823, 2.0918918918918923, 2.1018018018018023, 2.111711711711712, 2.121621621621622, 2.131531531531532, 2.1414414414414416, 2.1513513513513516, 2.1612612612612616, 2.1711711711711716, 2.1810810810810812, 2.1909909909909913, 2.2009009009009013, 2.2108108108108113, 2.220720720720721, 2.230630630630631, 2.240540540540541, 2.250450450450451, 2.2603603603603606, 2.2702702702702706, 2.2801801801801806, 2.2900900900900902, 2.3000000000000003, 2.3099099099099103, 2.3198198198198203, 2.32972972972973, 2.33963963963964, 2.34954954954955, 2.35945945945946, 2.3693693693693696, 2.3792792792792796, 2.3891891891891897, 2.3990990990990992, 2.4090090090090093, 2.4189189189189193, 2.4288288288288293, 2.438738738738739, 2.448648648648649, 2.458558558558559, 2.468468468468469, 2.4783783783783786, 2.4882882882882886, 2.4981981981981987, 2.5081081081081087, 2.5180180180180183, 2.5279279279279283, 2.5378378378378383, 2.547747747747748, 2.557657657657658, 2.567567567567568, 2.577477477477478, 2.5873873873873876, 2.5972972972972976, 2.6072072072072077, 2.6171171171171177, 2.6270270270270273, 2.6369369369369373, 2.6468468468468473, 2.656756756756757, 2.666666666666667, 2.676576576576577, 2.686486486486487, 2.6963963963963966, 2.7063063063063066, 2.7162162162162167, 2.7261261261261267, 2.7360360360360363, 2.7459459459459463, 2.7558558558558564, 2.765765765765766, 2.775675675675676, 2.785585585585586, 2.795495495495496, 2.8054054054054056, 2.8153153153153156, 2.8252252252252257, 2.8351351351351357, 2.8450450450450453, 2.8549549549549553, 2.8648648648648654, 2.8747747747747754, 2.884684684684685, 2.894594594594595, 2.904504504504505, 2.9144144144144146, 2.9243243243243247, 2.9342342342342347, 2.9441441441441447, 2.9540540540540543, 2.9639639639639643, 2.9738738738738744, 2.9837837837837844, 2.993693693693694, 3.003603603603604, 3.013513513513514, 3.0234234234234236, 3.0333333333333337, 3.0432432432432437, 3.0531531531531537, 3.0630630630630633, 3.0729729729729733, 3.0828828828828834, 3.0927927927927934, 3.102702702702703, 3.112612612612613, 3.122522522522523, 3.132432432432433, 3.1423423423423427, 3.1522522522522527, 3.1621621621621627, 3.1720720720720723, 3.1819819819819823, 3.1918918918918924, 3.2018018018018024, 3.211711711711712, 3.221621621621622, 3.231531531531532, 3.241441441441442, 3.2513513513513517, 3.2612612612612617, 3.2711711711711717, 3.2810810810810813, 3.2909909909909913, 3.3009009009009014, 3.3108108108108114, 3.320720720720721, 3.330630630630631, 3.340540540540541, 3.350450450450451, 3.3603603603603607, 3.3702702702702707, 3.3801801801801807, 3.3900900900900908, 3.4000000000000004, 3.4099099099099104, 3.4198198198198204, 3.42972972972973, 3.43963963963964, 3.44954954954955, 3.45945945945946, 3.4693693693693697, 3.4792792792792797, 3.4891891891891897, 3.4990990990990998, 3.5090090090090094, 3.5189189189189194, 3.5288288288288294, 3.538738738738739, 3.548648648648649, 3.558558558558559, 3.568468468468469, 3.5783783783783787, 3.5882882882882887, 3.5981981981981987, 3.608108108108109, 3.6180180180180184, 3.6279279279279284, 3.6378378378378384, 3.6477477477477485, 3.657657657657658, 3.667567567567568, 3.677477477477478, 3.6873873873873877, 3.6972972972972977, 3.7072072072072078, 3.717117117117118, 3.7270270270270274, 3.7369369369369374, 3.7468468468468474, 3.7567567567567575, 3.766666666666667, 3.776576576576577, 3.786486486486487, 3.7963963963963967, 3.8063063063063067, 3.8162162162162168, 3.826126126126127, 3.8360360360360364, 3.8459459459459464, 3.8558558558558564, 3.8657657657657665, 3.875675675675676, 3.885585585585586, 3.895495495495496, 3.9054054054054057, 3.9153153153153157, 3.9252252252252258, 3.935135135135136, 3.9450450450450454, 3.9549549549549554, 3.9648648648648654, 3.9747747747747755, 3.984684684684685, 3.994594594594595, 4.004504504504505, 4.014414414414415, 4.024324324324325, 4.034234234234234, 4.044144144144145, 4.054054054054054, 4.063963963963964, 4.0738738738738745, 4.083783783783784, 4.093693693693694, 4.103603603603604, 4.113513513513514, 4.123423423423423, 4.133333333333334, 4.143243243243243, 4.153153153153153, 4.163063063063063, 4.172972972972973, 4.182882882882883, 4.192792792792793, 4.202702702702703, 4.212612612612613, 4.222522522522523, 4.232432432432432, 4.242342342342343, 4.252252252252252, 4.262162162162162, 4.272072072072072, 4.281981981981982, 4.291891891891892, 4.301801801801802, 4.311711711711712, 4.321621621621622, 4.331531531531532, 4.341441441441441, 4.351351351351352, 4.361261261261261, 4.371171171171171, 4.381081081081081, 4.390990990990991, 4.4009009009009015, 4.410810810810811, 4.420720720720721, 4.430630630630631, 4.440540540540541, 4.45045045045045, 4.460360360360361, 4.47027027027027, 4.48018018018018, 4.49009009009009, 4.5, 4.5099099099099105, 4.51981981981982, 4.52972972972973, 4.53963963963964, 4.54954954954955, 4.559459459459459, 4.56936936936937, 4.579279279279279, 4.589189189189189, 4.599099099099099, 4.609009009009009, 4.6189189189189195, 4.628828828828829, 4.638738738738739, 4.648648648648649, 4.658558558558559, 4.668468468468468, 4.678378378378379, 4.688288288288288, 4.698198198198198, 4.708108108108108, 4.718018018018018, 4.7279279279279285, 4.737837837837838, 4.747747747747748, 4.757657657657658, 4.767567567567568, 4.777477477477477, 4.787387387387388, 4.797297297297297, 4.807207207207207, 4.817117117117117, 4.827027027027027, 4.8369369369369375, 4.846846846846847, 4.856756756756757, 4.866666666666667, 4.876576576576577, 4.886486486486486, 4.896396396396397, 4.906306306306306, 4.916216216216217, 4.926126126126126, 4.936036036036036, 4.9459459459459465, 4.955855855855856, 4.965765765765766, 4.975675675675676, 4.985585585585586, 4.995495495495495, 5.005405405405406, 5.015315315315315, 5.025225225225226, 5.035135135135135, 5.045045045045045, 5.0549549549549555, 5.064864864864865, 5.074774774774775, 5.084684684684685, 5.094594594594595, 5.104504504504504, 5.114414414414415, 5.124324324324324, 5.134234234234235, 5.1441441441441444, 5.154054054054054, 5.1639639639639645, 5.173873873873874, 5.183783783783784, 5.193693693693694, 5.203603603603604, 5.213513513513513, 5.223423423423424, 5.233333333333333, 5.243243243243244, 5.2531531531531535, 5.263063063063063, 5.2729729729729735, 5.282882882882883, 5.292792792792793, 5.302702702702703, 5.312612612612613, 5.322522522522522, 5.332432432432433, 5.342342342342342, 5.352252252252253, 5.3621621621621625, 5.372072072072072, 5.3819819819819825, 5.391891891891892, 5.401801801801802, 5.411711711711712, 5.421621621621622, 5.431531531531531, 5.441441441441442, 5.451351351351351, 5.461261261261262, 5.4711711711711715, 5.481081081081081, 5.4909909909909915, 5.500900900900901, 5.510810810810811, 5.520720720720721, 5.530630630630631, 5.540540540540541, 5.550450450450451, 5.56036036036036, 5.570270270270271, 5.5801801801801805, 5.59009009009009, 5.6000000000000005, 5.60990990990991, 5.61981981981982, 5.62972972972973, 5.63963963963964, 5.64954954954955, 5.65945945945946, 5.669369369369369, 5.67927927927928, 5.6891891891891895, 5.699099099099099, 5.7090090090090095, 5.718918918918919, 5.728828828828829, 5.738738738738739, 5.748648648648649, 5.758558558558559, 5.768468468468469, 5.778378378378378, 5.788288288288289, 5.7981981981981985, 5.808108108108108, 5.8180180180180185, 5.827927927927928, 5.837837837837838, 5.847747747747748, 5.857657657657658, 5.867567567567568, 5.877477477477478, 5.887387387387387, 5.897297297297298, 5.9072072072072075, 5.917117117117117, 5.9270270270270276, 5.936936936936937, 5.946846846846847, 5.956756756756757, 5.966666666666667, 5.976576576576577, 5.986486486486487, 5.996396396396396, 6.006306306306307, 6.0162162162162165, 6.026126126126126, 6.036036036036037, 6.045945945945946, 6.055855855855857, 6.065765765765766, 6.075675675675676, 6.085585585585586, 6.095495495495496, 6.105405405405405, 6.115315315315316, 6.1252252252252255, 6.135135135135135, 6.145045045045046, 6.154954954954955, 6.164864864864866, 6.174774774774775, 6.184684684684685, 6.194594594594595, 6.204504504504505, 6.2144144144144144, 6.224324324324325, 6.2342342342342345, 6.244144144144144, 6.254054054054055, 6.263963963963964, 6.273873873873875, 6.283783783783784, 6.293693693693694, 6.303603603603604, 6.313513513513514, 6.3234234234234235, 6.333333333333334, 6.3432432432432435, 6.353153153153153, 6.363063063063064, 6.372972972972973, 6.382882882882884, 6.392792792792793, 6.402702702702703, 6.412612612612613, 6.422522522522523, 6.4324324324324325, 6.442342342342343, 6.4522522522522525, 6.462162162162162, 6.472072072072073, 6.481981981981982, 6.491891891891893, 6.501801801801802, 6.511711711711712, 6.521621621621622, 6.531531531531532, 6.5414414414414415, 6.551351351351352, 6.5612612612612615, 6.571171171171171, 6.581081081081082, 6.590990990990991, 6.600900900900902, 6.610810810810811, 6.620720720720721, 6.630630630630631, 6.640540540540541, 6.6504504504504505, 6.660360360360361, 6.6702702702702705, 6.680180180180181, 6.690090090090091, 6.7, 6.709909909909911, 6.71981981981982, 6.72972972972973, 6.73963963963964, 6.74954954954955, 6.7594594594594595, 6.76936936936937, 6.7792792792792795, 6.78918918918919, 6.7990990990991, 6.809009009009009, 6.81891891891892, 6.828828828828829, 6.838738738738739, 6.848648648648649, 6.858558558558559, 6.8684684684684685, 6.878378378378379, 6.8882882882882885, 6.898198198198199, 6.908108108108109, 6.918018018018018, 6.927927927927929, 6.937837837837838, 6.947747747747748, 6.957657657657658, 6.967567567567568, 6.9774774774774775, 6.987387387387388, 6.9972972972972975, 7.007207207207208, 7.017117117117118, 7.027027027027027, 7.036936936936938, 7.046846846846847, 7.056756756756757, 7.066666666666667, 7.076576576576577, 7.0864864864864865, 7.096396396396397, 7.1063063063063066, 7.116216216216217, 7.126126126126127, 7.136036036036036, 7.145945945945947, 7.155855855855856, 7.165765765765766, 7.175675675675676, 7.185585585585586, 7.195495495495496, 7.205405405405406, 7.215315315315316, 7.225225225225226, 7.235135135135136, 7.245045045045045, 7.254954954954956, 7.264864864864865, 7.274774774774775, 7.284684684684685, 7.294594594594595, 7.304504504504505, 7.314414414414415, 7.324324324324325, 7.334234234234235, 7.344144144144145, 7.354054054054054, 7.363963963963965, 7.373873873873874, 7.383783783783784, 7.393693693693694, 7.403603603603604, 7.413513513513514, 7.423423423423424, 7.433333333333334, 7.443243243243244, 7.453153153153154, 7.463063063063063, 7.472972972972974, 7.482882882882883, 7.492792792792793, 7.502702702702703, 7.512612612612613, 7.522522522522523, 7.532432432432433, 7.542342342342343, 7.552252252252253, 7.562162162162163, 7.572072072072072, 7.581981981981983, 7.591891891891892, 7.601801801801802, 7.611711711711712, 7.621621621621622, 7.631531531531532, 7.641441441441442, 7.651351351351352, 7.661261261261262, 7.671171171171172, 7.681081081081081, 7.690990990990992, 7.700900900900901, 7.710810810810811, 7.720720720720721, 7.730630630630631, 7.740540540540541, 7.750450450450451, 7.760360360360361, 7.770270270270271, 7.780180180180181, 7.79009009009009, 7.800000000000001, 7.80990990990991, 7.819819819819821, 7.82972972972973, 7.83963963963964, 7.84954954954955, 7.85945945945946, 7.86936936936937, 7.87927927927928, 7.88918918918919, 7.899099099099099, 7.90900900900901, 7.918918918918919, 7.92882882882883, 7.938738738738739, 7.948648648648649, 7.958558558558559, 7.968468468468469, 7.978378378378379, 7.988288288288289, 7.998198198198199, 8.00810810810811, 8.018018018018019, 8.027927927927928, 8.03783783783784, 8.04774774774775, 8.057657657657659, 8.067567567567568, 8.077477477477478, 8.087387387387388, 8.097297297297299, 8.107207207207209, 8.117117117117118, 8.127027027027028, 8.136936936936937, 8.146846846846847, 8.156756756756756, 8.166666666666668, 8.176576576576577, 8.186486486486487, 8.196396396396397, 8.206306306306306, 8.216216216216218, 8.226126126126127, 8.236036036036037, 8.245945945945946, 8.255855855855856, 8.265765765765765, 8.275675675675677, 8.285585585585586, 8.295495495495496, 8.305405405405406, 8.315315315315315, 8.325225225225227, 8.335135135135136, 8.345045045045046, 8.354954954954955, 8.364864864864865, 8.374774774774774, 8.384684684684686, 8.394594594594595, 8.404504504504505, 8.414414414414415, 8.424324324324324, 8.434234234234236, 8.444144144144145, 8.454054054054055, 8.463963963963964, 8.473873873873874, 8.483783783783784, 8.493693693693695, 8.503603603603604, 8.513513513513514, 8.523423423423424, 8.533333333333333, 8.543243243243245, 8.553153153153154, 8.563063063063064, 8.572972972972973, 8.582882882882883, 8.592792792792794, 8.602702702702704, 8.612612612612613, 8.622522522522523, 8.632432432432433, 8.642342342342342, 8.652252252252254, 8.662162162162163, 8.672072072072073, 8.681981981981982, 8.691891891891892, 8.701801801801803, 8.711711711711713, 8.721621621621622, 8.731531531531532, 8.741441441441442, 8.751351351351351, 8.761261261261263, 8.771171171171172, 8.781081081081082, 8.790990990990991, 8.800900900900901, 8.810810810810812, 8.820720720720722, 8.830630630630631, 8.840540540540541, 8.85045045045045, 8.86036036036036, 8.870270270270272, 8.880180180180181, 8.89009009009009, 8.9, 8.90990990990991, 8.919819819819821, 8.929729729729731, 8.93963963963964, 8.94954954954955, 8.95945945945946, 8.96936936936937, 8.97927927927928, 8.98918918918919, 8.9990990990991, 9.00900900900901, 9.018918918918919, 9.02882882882883, 9.03873873873874, 9.04864864864865, 9.058558558558559, 9.068468468468469, 9.078378378378378, 9.08828828828829, 9.0981981981982, 9.108108108108109, 9.118018018018018, 9.127927927927928, 9.13783783783784, 9.147747747747749, 9.157657657657658, 9.167567567567568, 9.177477477477478, 9.187387387387387, 9.197297297297299, 9.207207207207208, 9.217117117117118, 9.227027027027027, 9.236936936936937, 9.246846846846848, 9.256756756756758, 9.266666666666667, 9.276576576576577, 9.286486486486487, 9.296396396396396, 9.306306306306308, 9.316216216216217, 9.326126126126127, 9.336036036036036, 9.345945945945946, 9.355855855855857, 9.365765765765767, 9.375675675675677, 9.385585585585586, 9.395495495495496, 9.405405405405405, 9.415315315315317, 9.425225225225226, 9.435135135135136, 9.445045045045045, 9.454954954954955, 9.464864864864866, 9.474774774774776, 9.484684684684686, 9.494594594594595, 9.504504504504505, 9.514414414414414, 9.524324324324326, 9.534234234234235, 9.544144144144145, 9.554054054054054, 9.563963963963964, 9.573873873873875, 9.583783783783785, 9.593693693693695, 9.603603603603604, 9.613513513513514, 9.623423423423423, 9.633333333333335, 9.643243243243244, 9.653153153153154, 9.663063063063063, 9.672972972972973, 9.682882882882884, 9.692792792792794, 9.702702702702704, 9.712612612612613, 9.722522522522523, 9.732432432432434, 9.742342342342344, 9.752252252252253, 9.762162162162163, 9.772072072072072, 9.781981981981982, 9.791891891891893, 9.801801801801803, 9.811711711711713, 9.821621621621622, 9.831531531531532, 9.841441441441443, 9.851351351351353, 9.861261261261262, 9.871171171171172, 9.881081081081081, 9.890990990990991, 9.900900900900902, 9.910810810810812, 9.920720720720722, 9.930630630630631, 9.94054054054054, 9.950450450450452, 9.960360360360362, 9.970270270270271, 9.98018018018018, 9.99009009009009, 10.0], "expected": [0.48808848170151553, 0.4869680959334418, 0.48585779844963695, 0.4847574211445585, 0.48366680005044566, 0.48258577520226276, 0.4815141905081331, 0.48045189362499346, 0.47939873583922304, 0.47835457195200415, 0.47731926016919324, 0.47629266199548936, 0.47527464213269793, 0.47426506838190435, 0.47326381154937236, 0.472270745356002, 0.4712857463501854, 0.4703086938239047, 0.4693394697319325, 0.46837795861399456, 0.4674240475197656, 0.4664776259365753, 0.4655385857197078, 0.4646068210251813, 0.4636822282449053, 0.4627647059441123, 0.46185415480097003, 0.4609504775482829, 0.4600535789171963, 0.4591633655828209, 0.4582797461117001, 0.45740263091104344, 0.45653193217965865, 0.45566756386051027, 0.45480944159484293, 0.4539574826778068, 0.45311160601552825, 0.4522717320835644, 0.4514377828866942, 0.45060968191999123, 0.44978735413112997, 0.44897072588387793, 0.4481597249227319, 0.4473542803386525, 0.446554322535857, 0.4457597831996314, 0.4449705952651259, 0.4441866928870954, 0.44340801141055175, 0.4426344873422951, 0.44186605832329423, 0.44110266310188107, 0.44034424150773716, 0.439590734426638, 0.43884208377593387, 0.4380982324807363, 0.43735912445079134, 0.43662470455801183, 0.4358949186146489, 0.4351697133520781, 0.43444903640018234, 0.43373283626731013, 0.43302106232079, 0.4323136647679824, 0.4316105946378526, 0.4309118037630458, 0.4302172447624497, 0.4295268710242259, 0.42884063668929917, 0.4281584966352877, 0.42748040646085794, 0.4268063224704944, 0.42613620165966987, 0.42547000170040356, 0.4248076809271922, 0.42414919832330905, 0.4234945135074525, 0.42284358672073774, 0.4221963788140194, 0.4215528512355349, 0.42091296601886075, 0.4202766857711694, 0.4196439736617804, 0.4190147934109944, 0.41838910927920425, 0.41776688605627366, 0.41714808905117473, 0.4165326840818787, 0.415920637465492, 0.415311916008628, 0.4147064869980119, 0.41410431819130905, 0.413505377808171, 0.4129096345214938, 0.4123170574488815, 0.4117276161443095, 0.411141280589984, 0.4105580211883867, 0.4099778087545071, 0.4094006145082517, 0.40882641006702614, 0.4082551674384875, 0.40768685901346047, 0.40712145755901363, 0.40655893621169115, 0.4059992684708961, 0.4054424281924197, 0.40488838958211726, 0.4043371271897198, 0.40378861590278425, 0.40324283094077573, 0.4026997478492772, 0.40215934249432655, 0.401621591056875, 0.40108647002736425, 0.40055395620042117, 0.4000240266696634, 0.39949665882261787, 0.3989718303357441, 0.3984495191695638, 0.39792970356389257, 0.39741236203317043, 0.3968974733618887, 0.3963850166001141, 0.3958749710591028, 0.395367316307004, 0.39486203216465493, 0.39435909870145647, 0.3938584962313365, 0.39336020530879257, 0.3928642067250154, 0.3923704815040881, 0.3918790108992648, 0.3913897763893178, 0.3909027596749634, 0.39041794267535285, 0.3899353075246348, 0.389454836568585, 0.3889765123613009, 0.38850031766196025, 0.3880262354316437, 0.3875542488302166, 0.3870843412132726, 0.38661649612913473, 0.38615069731591267, 0.3856869286986181, 0.3852251743863321, 0.38476541866942776, 0.38430764601684386, 0.3838518410734096, 0.3833979886572193, 0.38294607375705647, 0.3824960815298635, 0.38204799729826, 0.38160180654810616, 0.381157494926109, 0.3807150482374756, 0.38027445244360447, 0.37983569365982267, 0.379398758153162, 0.37896363234017394, 0.37853030278478594, 0.37809875619619465, 0.37766897942679595, 0.377240959470154, 0.37681468345900276, 0.3763901386632851, 0.37596731248822596, 0.3755461924724378, 0.37512676628606034, 0.37470902172893217, 0.3742929467287932, 0.37387852933951937, 0.3734657577393865, 0.3730546202293641, 0.3726451052314387, 0.37223720128696525, 0.3718308970550464, 0.3714261813109389, 0.3710230429444877, 0.3706214709585849, 0.37022145446765586, 0.369822982696169, 0.369426044977172, 0.3690306307508508, 0.36863672956311233, 0.3682443310641924, 0.36785342500728385, 0.36746400124718837, 0.36707604973899066, 0.366689560536753, 0.3663045237922316, 0.3659209297536135, 0.3655387687642737, 0.3651580312615521, 0.3647787077755501, 0.3644007889279463, 0.36402426543083105, 0.3636491280855589, 0.3632753677816199, 0.3629029754955283, 0.36253194228972835, 0.36216225931151774, 0.36179391779198705, 0.36142690904497643, 0.361061224466048, 0.36069685553147424, 0.3603337937972416, 0.35997203089806973, 0.3596115585464455, 0.3592523685316723, 0.35889445271893233, 0.3585378030483649, 0.3581824115341576, 0.35782827026365077, 0.35747537139645674, 0.357123707163591, 0.3567732698666165, 0.356424051876802, 0.3560760456342906, 0.3557292436472829, 0.3553836384912303, 0.35503922280804123, 0.35469598930529933, 0.35435393075549143, 0.35401303999524864, 0.35367330992459733, 0.3533347335062216, 0.35299730376473537, 0.35266101378596637, 0.35232585671624944, 0.3519918257617304, 0.35165891418767925, 0.351327115317814, 0.35099642253363367, 0.35066682927376025, 0.3503383290332913, 0.35001091536316004, 0.3496845818695055, 0.3493593222130515, 0.3490351301084932, 0.3487119993238941, 0.34838992368008914, 0.3480688970500981, 0.3477489133585459, 0.34742996658109104, 0.34711205074386264, 0.34679515992290383, 0.3464792882436242, 0.3461644298802588, 0.34585057905533484, 0.3455377300391451, 0.3452258771492294, 0.34491501474986264, 0.3446051372515489, 0.3442962391105235, 0.34398831482826153, 0.3436813589509921, 0.3433753660692202, 0.3430703308172535, 0.34276624787273746, 0.3424631119561941, 0.3421609178305688, 0.3418596603007824, 0.3415593342132888, 0.3412599344556388, 0.3409614559560493, 0.3406638936829787, 0.3403672426447069, 0.3400714978889217, 0.3397766545023096, 0.3394827076101522, 0.33918965237592896, 0.3388974840009224, 0.3386061977238312, 0.3383157888203862, 0.33802625260297253, 0.33773758442025564, 0.33744977965681244, 0.3371628337327673, 0.33687674210343244, 0.33659150025895257, 0.33630710372395456, 0.3360235480572007, 0.335740828851247, 0.3354589417321054, 0.33517788235891016, 0.3348976464235887, 0.33461822965053584, 0.3343396277962927, 0.33406183664922917, 0.3337848520292304, 0.33350866978738725, 0.3332332858056901, 0.3329586959967266, 0.3326848963033838, 0.33241188269855243, 0.3321396511848358, 0.33186819779426185, 0.3315975185879991, 0.33132760965607516, 0.3310584671170996, 0.33079008711798924, 0.33052246583369715, 0.33025559946694505, 0.32998948424795815, 0.3297241164342041, 0.3294594923101344, 0.32919560818692895, 0.32893246040224333, 0.32867004531995986, 0.3284083593299408, 0.3281473988477843, 0.3278871603145846, 0.3276276401966925, 0.3273688349854815, 0.32711074119711414, 0.32685335537231275, 0.32659667407613197, 0.3263406938977343, 0.32608541145016723, 0.3258308233701449, 0.32557692631783014, 0.32532371697662044, 0.32507119205293555, 0.3248193482760075, 0.3245681823976737, 0.3243176911921717, 0.32406787145593646, 0.3238187200074001, 0.32357023368679316, 0.32332240935594914, 0.32307524389811015, 0.3228287342177362, 0.32258287724031454, 0.32233766991217283, 0.3220931092002939, 0.3218491920921319, 0.32160591559543167, 0.3213632767380485, 0.3211212725677723, 0.32087990015215, 0.32063915657831443, 0.32039903895281074, 0.32015954440142785, 0.3199206700690304, 0.319682413119392, 0.31944477073503136, 0.31920774011705033, 0.3189713184849717, 0.31873550307658166, 0.3185002911477715, 0.3182656799723824, 0.3180316668420509, 0.3177982490660569, 0.3175654239711731, 0.31733318890151524, 0.31710154121839507, 0.3168704783001737, 0.3166399975421179, 0.31641009635625644, 0.31618077217123897, 0.31595202243219567, 0.3157238446005994, 0.3154962361541274, 0.3152691945865271, 0.3150427174074806, 0.31481680214247265, 0.31459144633265845, 0.31436664753473426, 0.3141424033208077, 0.3139187112782708, 0.3136955690096743, 0.31347297413260156, 0.3132509242795461, 0.31302941709778836, 0.3128084502492751, 0.3125880214104996, 0.3123681282723823, 0.3121487685401544, 0.3119299399332405, 0.31171164018514447, 0.3114938670433343, 0.31127661826913067, 0.3110598916375941, 0.31084368493741543, 0.3106279959708049, 0.3104128225533854, 0.31019816251408416, 0.30998401369502676, 0.3097703739514317, 0.30955724115150707, 0.3093446131763465, 0.30913248791982706, 0.30892086328850926, 0.30870973720153544, 0.308499107590532, 0.30828897239951014, 0.30807932958476925, 0.3078701771148003, 0.30766151297019106, 0.307453335143531, 0.3072456416393183, 0.30703843047386686, 0.306831699675215, 0.3066254472830343, 0.3064196713485397, 0.30621436993440093, 0.3060095411146532, 0.30580518297461096, 0.3056012936107804, 0.3053978711307745, 0.3051949136532274, 0.30499241930771054, 0.30479038623464977, 0.3045888125852422, 0.3043876965213745, 0.3041870362155423, 0.3039868298507695, 0.30378707562052915, 0.3035877717286645, 0.3033889163893112, 0.30319050782681944, 0.30299254427567807, 0.3027950239804386, 0.3025979451956398, 0.3024013061857333, 0.30220510522501004, 0.30200934059752704, 0.30181401059703494, 0.30161911352690635, 0.30142464770006466, 0.30123061143891366, 0.30103700307526765, 0.30084382095028256, 0.3006510634143872, 0.3004587288272151, 0.30026681555753765, 0.30007532198319736, 0.2998842464910413, 0.29969358747685654, 0.29950334334530404, 0.29931351250985555, 0.29912409339272894, 0.2989350844248255, 0.2987464840456673, 0.298558290703335, 0.29837050285440647, 0.29818311896389593, 0.29799613750519327, 0.29780955696000444, 0.29762337581829196, 0.29743759257821645, 0.2972522057460776, 0.2970672138362573, 0.2968826153711616, 0.2966984088811644, 0.29651459290455096, 0.2963311659874625, 0.29614812668384016, 0.2959654735553708, 0.2957832051714327, 0.29560132010904105, 0.29541981695279546, 0.2952386942948262, 0.29505795073474234, 0.2948775848795794, 0.2946975953437477, 0.29451798074898133, 0.2943387397242879, 0.29415987090589746, 0.29398137293721327, 0.29380324446876205, 0.2936254841581454, 0.2934480906699906, 0.2932710626759034, 0.2930943988544192, 0.29291809789095663, 0.29274215847776985, 0.2925665793139025, 0.2923913591051415, 0.29221649656397125, 0.29204199040952783, 0.29186783936755495, 0.2916940421703581, 0.29152059755676135, 0.2913475042720631, 0.2911747610679923, 0.29100236670266566, 0.2908303199405447, 0.2906586195523934, 0.29048726431523597, 0.29031625301231523, 0.2901455844330512, 0.28997525737299984, 0.2898052706338128, 0.28963562302319645, 0.2894663133548725, 0.28929734044853755, 0.2891287031298242, 0.2889604002302618, 0.2887924305872378, 0.288624793043959, 0.288457486449414, 0.2882905096583346, 0.2881238615311589, 0.28795754093399395, 0.2877915467385786, 0.2876258778222472, 0.2874605330678931, 0.2872955113639327, 0.2871308116042694, 0.2869664326882589, 0.2868023735206731, 0.28663863301166587, 0.28647521007673815, 0.28631210363670384, 0.28614931261765525, 0.28598683595093005, 0.285824672573077, 0.2856628214258233, 0.28550128145604103, 0.28534005161571496, 0.2851791308619097, 0.28501851815673757, 0.28485821246732673, 0.2846982127657892, 0.2845385180291898, 0.28437912723951453, 0.2842200393836397, 0.2840612534533015, 0.28390276844506496, 0.283744583360294, 0.2835866972051213, 0.28342910899041895, 0.28327181773176785, 0.2831148224494296, 0.28295812216831645, 0.2828017159179627, 0.28264560273249645, 0.2824897816506104, 0.2823342517155343, 0.2821790119750064, 0.2820240614812461, 0.28186939929092597, 0.28171502446514474, 0.2815609360693998, 0.28140713317356053, 0.28125361485184136, 0.28110038018277517, 0.2809474282491869, 0.28079475813816784, 0.28064236894104927, 0.28049025975337655, 0.2803384296748839, 0.2801868778094689, 0.2800356032651673, 0.27988460515412794, 0.2797338825925882, 0.27958343470084895, 0.27943326060325047, 0.27928335942814797, 0.27913373030788774, 0.27898437237878293, 0.27883528478109015, 0.27868646665898583, 0.2785379171605424, 0.27838963543770584, 0.2782416206462722, 0.2780938719458648, 0.2779463884999114, 0.27779916947562194, 0.27765221404396606, 0.277505521379651, 0.27735909066109904, 0.2772129210704265, 0.27706701179342147, 0.27692136201952194, 0.2767759709417954, 0.2766308377569165, 0.27648596166514655, 0.2763413418703127, 0.2761969775797865, 0.2760528680044642, 0.27590901235874515, 0.2757654098605126, 0.27562205973111253, 0.27547896119533416, 0.27533611348138975, 0.2751935158208953, 0.2750511674488501, 0.27490906760361833, 0.27476721552690847, 0.2746256104637552, 0.2744842516624996, 0.2743431383747708, 0.2742022698554664, 0.2740616453627345, 0.27392126415795487, 0.2737811255057206, 0.2736412286738198, 0.2735015729332173, 0.27336215755803706, 0.2732229818255436, 0.27308404501612504, 0.2729453464132746, 0.2728068853035739, 0.27266866097667486, 0.272530672725283, 0.27239291984513964, 0.27225540163500556, 0.27211811739664354, 0.2719810664348018, 0.2718442480571975, 0.2717076615744995, 0.2715713063003127, 0.271435181551161, 0.27129928664647174, 0.27116362090855894, 0.271028183662608, 0.2708929742366588, 0.2707579919615909, 0.27062323617110734, 0.2704887062017192, 0.2703544013927301, 0.2702203210862206, 0.2700864646270337, 0.26995283136275855, 0.2698194206437165, 0.2696862318229453, 0.26955326425618475, 0.26942051730186195, 0.2692879903210761, 0.2691556826775848, 0.26902359373778867, 0.26889172287071805, 0.2687600694480176, 0.268628632843933, 0.2684974124352964, 0.26836640760151287, 0.26823561772454607, 0.26810504218890446, 0.2679746803816282, 0.2678445316922748, 0.2677145955129064, 0.2675848712380751, 0.2674553582648112, 0.26732605599260884, 0.26719696382341296, 0.2670680811616066, 0.2669394074139979, 0.26681094198980676, 0.26668268430065234, 0.2665546337605404, 0.2664267897858504, 0.2662991517953231, 0.26617171921004795, 0.26604449145345105, 0.26591746795128207, 0.26579064813160297, 0.26566403142477496, 0.26553761726344716, 0.2654114050825439, 0.26528539431925335, 0.26515958441301535, 0.2650339748055099, 0.264908564940645, 0.2647833542645457, 0.2646583422255417, 0.264533528274157, 0.26440891186309734, 0.2642844924472395, 0.2641602694836203, 0.26403624243142454, 0.26391241075197475, 0.26378877390871974, 0.2636653313672237, 0.26354208259515516, 0.2634190270622767, 0.2632961642404334, 0.2631734936035428, 0.26305101462758396, 0.2629287267905869, 0.2628066295726222, 0.2626847224557907, 0.2625630049242129, 0.2624414764640188, 0.2623201365633375, 0.2621989847122876, 0.262078020402966, 0.26195724312943935, 0.26183665238773257, 0.26171624767582014, 0.2615960284936154, 0.2614759943429616, 0.261356144727621, 0.2612364791532665, 0.2611169971274712, 0.2609976981596991, 0.2608785817612958, 0.2607596474454783, 0.26064089472732704, 0.2605223231237751, 0.2604039321535998, 0.2602857213374137, 0.2601676901976542, 0.2600498382585762, 0.25993216504624184, 0.25981467008851167, 0.2596973529150365, 0.2595802130572474, 0.2594632500483479, 0.2593464634233046, 0.2592298527188388, 0.25911341747341776, 0.2589971572272457, 0.2588810715222564, 0.2587651599021033, 0.2586494219121519, 0.25853385709947135, 0.2584184650128256, 0.2583032452026659, 0.25818819722112163, 0.258073320621993, 0.2579586149607425, 0.2578440797944865, 0.25772971468198796, 0.25761551918364756, 0.2575014928614967, 0.25738763527918873, 0.25727394600199166, 0.2571604245967804, 0.25704707063202814, 0.2569338836778, 0.2568208633057443, 0.25670800908908536, 0.25659532060261603, 0.25648279742268953, 0.2563704391272131, 0.25625824529563934, 0.25614621550895966, 0.2560343493496966, 0.25592264640189644, 0.2558111062511224, 0.25569972848444666, 0.25558851269044375, 0.25547745845918346, 0.2553665653822233, 0.255255833052602, 0.2551452610648317, 0.255034849014892, 0.25492459650022214, 0.2548145031197147, 0.2547045684737085, 0.25459479216398145, 0.2544851737937448, 0.2543757129676353, 0.25426640929170896, 0.2541572623734346, 0.2540482718216869, 0.25393943724674006, 0.253830758260261, 0.2537222344753033, 0.25361386550630016, 0.25350565096905836, 0.253397590480752, 0.25328968365991567, 0.2531819301264384, 0.2530743295015576, 0.25296688140785234, 0.2528595854692376, 0.2527524413109578, 0.25264544855958077, 0.2525386068429915, 0.25243191579038654, 0.2523253750322675, 0.2522189842004348, 0.2521127429279827, 0.25200665084929247, 0.2519007076000266, 0.2517949128171236, 0.2516892661387912, 0.25158376720450143, 0.2514784156549842, 0.2513732111322221, 0.251268153279444, 0.25116324174112037, 0.2510584761629568, 0.2509538561918887, 0.25084938147607555, 0.2507450516648961, 0.25064086640894156, 0.2505368253600114, 0.25043292817110696, 0.2503291744964268, 0.25022556399136053, 0.2501220963124841, 0.2500187711175543, 0.24991558806550296, 0.24981254681643258, 0.24970964703161036, 0.2496068883734633, 0.2495042705055729, 0.2494017930926701, 0.24929945580063015, 0.2491972582964674, 0.24909520024833032, 0.24899328132549653, 0.2488915011983675, 0.24878985953846405, 0.24868835601842085, 0.24858699031198186, 0.24848576209399537, 0.24838467104040873, 0.2482837168282643, 0.2481828991356938, 0.24808221764191396, 0.24798167202722163, 0.24788126197298926, 0.2477809871616596, 0.24768084727674164, 0.24758084200280553, 0.24748097102547822, 0.24738123403143847, 0.24728163070841283, 0.24718216074517027, 0.2470828238315185, 0.24698361965829888, 0.24688454791738204, 0.24678560830166332, 0.24668680050505887, 0.24658812422250034, 0.2464895791499312, 0.24639116498430189, 0.24629288142356612, 0.24619472816667554, 0.24609670491357638, 0.24599881136520466, 0.24590104722348194, 0.24580341219131133, 0.24570590597257322, 0.24560852827212054, 0.24551127879577567, 0.2454141572503251, 0.2453171633435162, 0.24522029678405277, 0.24512355728159071, 0.24502694454673443, 0.24493045829103266, 0.24483409822697416, 0.24473786406798384, 0.24464175552841913, 0.2445457723235655, 0.24444991416963313, 0.24435418078375223, 0.24425857188396985, 0.2441630871892454, 0.24406772641944746, 0.24397248929534943, 0.24387737553862562, 0.24378238487184828, 0.2436875170184827, 0.24359277170288438, 0.2434981486502946, 0.2434036475868374, 0.24330926823951515, 0.24321501033620543, 0.24312087360565715, 0.2430268577774868, 0.2429329625821752, 0.24283918775106342, 0.24274553301634955, 0.242651998111085, 0.24255858276917092, 0.24246528672535475, 0.24237210971522663, 0.24227905147521606, 0.24218611174258828, 0.24209329025544069, 0.24200058675269978, 0.24190800097411735, 0.24181553266026742, 0.24172318155254244, 0.2416309473931503, 0.24153882992511091, 0.24144682889225252, 0.2413549440392089, 0.24126317511141565, 0.2411715218551073, 0.24107998401731343, 0.24098856134585608, 0.24089725358934616, 0.24080606049718034, 0.24071498181953765, 0.24062401730737656, 0.24053316671243172, 0.2404424297872107, 0.24035180628499092, 0.24026129595981668, 0.2401708985664957, 0.24008061386059634, 0.2399904415984444, 0.23990038153712018, 0.23981043343445518, 0.2397205970490292, 0.23963087214016743, 0.23954125846793745, 0.23945175579314576, 0.2393623638773357, 0.23927308248278376, 0.23918391137249675, 0.2390948503102093, 0.2390058990603802, 0.23891705738819038, 0.23882832505953933, 0.2387397018410426, 0.23865118750002878, 0.2385627818045367, 0.2384744845233125, 0.2383862954258073, 0.23829821428217368, 0.23821024086326334, 0.23812237494062413, 0.23803461628649752, 0.23794696467381576, 0.23785941987619896, 0.23777198166795266, 0.23768464982406487, 0.23759742412020354, 0.2375103043327138, 0.2374232902386154, 0.2373363816155999, 0.23724957824202808, 0.23716287989692764, 0.23707628635998979, 0.23698979741156753, 0.23690341283267266, 0.23681713240497296, 0.23673095591079016, 0.23664488313309692, 0.23655891385551456, 0.23647304786231055, 0.23638728493839573, 0.23630162486932205, 0.23621606744128007, 0.2361306124410962, 0.23604525965623077, 0.23596000887477495, 0.23587485988544873, 0.23578981247759842, 0.23570486644119396, 0.23562002156682713, 0.23553527764570825, 0.2354506344696647, 0.2353660918311379, 0.2352816495231812, 0.2351973073394576, 0.23511306507423715, 0.2350289225223951, 0.23494487947940892, 0.2348609357413566, 0.23477709110491388, 0.2346933453673522, 0.2346096983265366, 0.2345261497809231, 0.2344426995295566, 0.2343593473720687, 0.23427609310867528, 0.23419293654017448, 0.23410987746794454, 0.23402691569394116, 0.2339440510206959, 0.23386128325131356, 0.2337786121894702, 0.2336960376394109, 0.23361355940594786, 0.2335311772944577, 0.2334488911108799, 0.2333667006617145, 0.2332846057540198, 0.2332026061954105, 0.23312070179405542, 0.23303889235867573, 0.2329571776985423, 0.23287555762347437, 0.23279403194383694, 0.23271260047053885, 0.23263126301503104, 0.23255001938930397, 0.2324688694058863, 0.23238781287784205, 0.2323068496187696, 0.2322259794427987, 0.2321452021645893, 0.23206451759932906, 0.23198392556273142, 0.23190342587103427, 0.23182301834099714, 0.23174270278989972, 0.23166247903554002]} diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/test/fixtures/python/runner.py b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/test/fixtures/python/runner.py new file mode 100644 index 000000000000..fc231f83adda --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/test/fixtures/python/runner.py @@ -0,0 +1,71 @@ +#!/usr/bin/env python +# +# @license Apache-2.0 +# +# Copyright (c) 2025 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. + +"""Generate fixtures.""" + +import os +import json +import numpy as np +from scipy.stats import bradford + +# Get the file path: +FILE = os.path.realpath(__file__) + +# Extract the directory in which this file resides: +DIR = os.path.dirname(FILE) + + +def gen(x, name): + """Generate fixture data and write to file. + + # Arguments + + * `x`: domain + * `name::str`: output filename + + # Examples + + ``` python + python> x = linspace(0.1, 100, 2001) + python> gen(x, './data.json') + ``` + """ + y = bradford.median(x) + + # Store data to be written to file as a dictionary: + data = { + "x": x.tolist(), + "expected": y.tolist() + } + + # Based on the script directory, create an output filepath: + filepath = os.path.join(DIR, name) + + # Write the data to the output filepath as JSON: + with open(filepath, "w", encoding="utf-8") as outfile: + json.dump(data, outfile) + + +def main(): + """Generate fixture data.""" + c = np.linspace(0.1, 10, 1000) + gen(c, "data.json") + + +if __name__ == "__main__": + main() diff --git a/lib/node_modules/@stdlib/stats/base/dists/bradford/median/test/test.js b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/test/test.js new file mode 100644 index 000000000000..ade2cc592b55 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/bradford/median/test/test.js @@ -0,0 +1,86 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2025 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var tape = require( 'tape' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var abs = require( '@stdlib/math/base/special/abs' ); +var NINF = require( '@stdlib/constants/float64/ninf' ); +var EPS = require( '@stdlib/constants/float64/eps' ); +var median = require( './../lib' ); + + +// FIXTURES // + +var data = require( './fixtures/python/data.json' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof median, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'if provided `NaN` for `c`, the function returns `NaN`', function test( t ) { + var v = median( NaN ); + t.equal( isnan( v ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'if provided `c <= 0`, the function returns `NaN`', function test( t ) { + var v; + + v = median( 0.0 ); + t.equal( isnan( v ), true, 'returns expected value' ); + + v = median( -1.0 ); + t.equal( isnan( v ), true, 'returns expected value' ); + + v = median( NINF ); + t.equal( isnan( v ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns the median of a Bradford distribution', function test( t ) { + var expected; + var delta; + var tol; + var i; + var c; + var y; + + expected = data.expected; + c = data.x; + for ( i = 0; i < expected.length; i++ ) { + y = median( c[i] ); + if ( y === expected[i] ) { + t.equal( y, expected[i], 'c: '+c[i]+', y: '+y+', expected: '+expected[i] ); + } else { + delta = abs( y - expected[ i ] ); + tol = 7.2 * EPS * abs( expected[ i ] ); + t.ok( delta <= tol, 'within tolerance. c: '+c[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' ); + } + } + t.end(); +});