diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/README.md b/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/README.md
index a7a07a0ab581..c807c7ade160 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/README.md
@@ -152,9 +152,6 @@ for ( i = 0; i < 10; i++ ) {
 
 <!-- /.examples -->
 
-<!-- Section to include cited references. If references are included, add a horizontal rule *before* the section. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
-
-<section class="references">
 
 <!-- C interface documentation. -->
 
@@ -249,6 +246,18 @@ int main( void ) {
 
 </section>
 
+<!-- /.examples -->
+
+</section>
+
+<!-- /.c -->
+
+<!-- Section to include cited references. If references are included, add a horizontal rule *before* the section. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
+
+<section class="references">
+
+</section>
+
 <!-- /.references -->
 
 <!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. -->
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/benchmark/c/benchmark.c
index 06efcf2aaca6..8e4d9e7bf461 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/benchmark/c/benchmark.c
@@ -93,12 +93,12 @@ static double random_uniform( const double min, const double max ) {
 * @return elapsed time in seconds
 */
 static double benchmark( void ) {
-	double elapsed;
-	double p[ 100 ];
-	double mu[ 100 ];
 	double sigma[ 100 ];
-	double y;
+	double mu[ 100 ];
+	double p[ 100 ];
+	double elapsed;
 	double t;
+	double y;
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/manifest.json b/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/manifest.json
index 749019a51b48..709351d4c22f 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/manifest.json
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/napi/ternary",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/math/base/special/erfinv"
+        "@stdlib/math/base/special/erfinv",
+        "@stdlib/constants/float64/sqrt-two"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/erfinv",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/constants/float64/eps",
+        "@stdlib/constants/float64/sqrt-two"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/erfinv",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/constants/float64/eps",
+        "@stdlib/constants/float64/sqrt-two"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/src/main.c b/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/src/main.c
index 86e8c5b8fd30..039d7d8b12d3 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/src/main.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/quantile/src/main.c
@@ -20,6 +20,7 @@
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/sqrt.h"
 #include "stdlib/math/base/special/erfinv.h"
+#include "stdlib/constants/float64/sqrt_two.h"
 
 /**
 * Evaluates the quantile function for a normal distribution with mean `mu` and standard deviation `sigma` at a probability `p`.
@@ -47,5 +48,5 @@ double stdlib_base_dists_normal_quantile( const double p, const double mu, const
 	if ( sigma == 0.0 ) {
 		return mu;
 	}
-	return mu + ( ( sigma * stdlib_base_sqrt( 2.0 ) ) * stdlib_base_erfinv( ( 2.0 * p ) - 1.0 ) );
+	return mu + ( ( sigma * STDLIB_CONSTANT_FLOAT64_SQRT2 ) * stdlib_base_erfinv( ( 2.0 * p ) - 1.0 ) );
 }