diff --git a/lib/node_modules/@stdlib/blas/base/dsyr/README.md b/lib/node_modules/@stdlib/blas/base/dsyr/README.md
index 2c37c5482c79..e20d2152f1d1 100644
--- a/lib/node_modules/@stdlib/blas/base/dsyr/README.md
+++ b/lib/node_modules/@stdlib/blas/base/dsyr/README.md
@@ -307,7 +307,7 @@ int main( void ) {
 
     // Print the result:
     for ( int i = 0; i < N*N; i++ ) {
-        printf( "A2[ %i ] = %lf\n", i, A[ i ] );
+        printf( "A2[ %i ] = %lf\n", i, A2[ i ] );
     }
 }
 ```
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/README.md b/lib/node_modules/@stdlib/blas/base/sdsdot/README.md
index d6b8f9131d58..c1eb8301d180 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/README.md
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/README.md
@@ -213,7 +213,7 @@ const float x[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
 const float y[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
 
 float v = c_sdsdot( 5, 0.0f, x, 1, y, -1 );
-// returns -120.0f
+// returns 35.0f
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/blas/base/sger/README.md b/lib/node_modules/@stdlib/blas/base/sger/README.md
index 3220d4a0614c..737fe0646aee 100644
--- a/lib/node_modules/@stdlib/blas/base/sger/README.md
+++ b/lib/node_modules/@stdlib/blas/base/sger/README.md
@@ -335,7 +335,7 @@ int main( void ) {
    }
 
    // Perform operation using alternative indexing semantics:
-   c_sger( CblasRowMajor, M, N, 1.0f, x, strideX, 0, y, 0, strideY, A, N, 1, 0 );
+   c_sger_ndarray( M, N, 1.0f, x, strideX, 0, y, strideY, 0, A, N, 1, 0 );
 
    // Print the result:
    for ( int i = 0; i < M; i++ ) {
diff --git a/lib/node_modules/@stdlib/blas/base/sspr/README.md b/lib/node_modules/@stdlib/blas/base/sspr/README.md
index 0e5f73306609..85bd2f57a9e7 100644
--- a/lib/node_modules/@stdlib/blas/base/sspr/README.md
+++ b/lib/node_modules/@stdlib/blas/base/sspr/README.md
@@ -225,7 +225,7 @@ Performs the symmetric rank 1 operation `A = α*x*x^T + A` where `α` is a scala
 float AP[] = { 1.0f, 2.0f, 3.0f, 1.0f, 2.0f, 1.0f };
 const float x[] = { 1.0f, 2.0f, 3.0f };
 
-c_sspr_ndarray( CblasColMajor, CblasUpper, 3, 1.0f, x, 1, AP, 1, 0 );
+c_sspr_ndarray( CblasColMajor, CblasUpper, 3, 1.0f, x, 1, 0, AP, 1, 0 );
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/blas/base/ssyr/README.md b/lib/node_modules/@stdlib/blas/base/ssyr/README.md
index 0526eef6cb2a..1ed5fad234b1 100644
--- a/lib/node_modules/@stdlib/blas/base/ssyr/README.md
+++ b/lib/node_modules/@stdlib/blas/base/ssyr/README.md
@@ -307,7 +307,7 @@ int main( void ) {
 
     // Print the result:
     for ( int i = 0; i < N*N; i++ ) {
-        printf( "A2[ %i ] = %f\n", i, A[ i ] );
+        printf( "A2[ %i ] = %f\n", i, A2[ i ] );
     }
 }
 ```
diff --git a/lib/node_modules/@stdlib/blas/ext/base/cfill/README.md b/lib/node_modules/@stdlib/blas/ext/base/cfill/README.md
index 3b560d4c1002..d56850103ad5 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/cfill/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/cfill/README.md
@@ -221,7 +221,7 @@ Fills a single-precision complex floating-point strided array `X` with a specifi
 float x[] = { 1.0f, 2.0f, 3.0f, 4.0f };
 const stdlib_complex64_t alpha = stdlib_complex64( 2.0f, 2.0f );
 
-stdlib_strided_cfill_ndarray( 4, alpha, (stdlib_complex64_t *x), 1, 0 );
+stdlib_strided_cfill_ndarray( 4, alpha, (stdlib_complex64_t *)x, 1, 0 );
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/blas/ext/base/dasumpw/README.md b/lib/node_modules/@stdlib/blas/ext/base/dasumpw/README.md
index e901a320238e..d9ce0e8a4442 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/dasumpw/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/dasumpw/README.md
@@ -192,7 +192,7 @@ console.log( v );
 Computes the sum of absolute values ([_L1_ norm][l1norm]) of double-precision floating-point strided array elements using pairwise summation.
 
 ```c
-const double x[] = { 1.0, 2.0, 3.0, 4.0 }
+const double x[] = { 1.0, 2.0, 3.0, 4.0 };
 
 double v = stdlib_strided_dasumpw( 4, x, 1 );
 // returns 10.0
@@ -213,7 +213,7 @@ double stdlib_strided_dasumpw( const CBLAS_INT N, const double *X, const CBLAS_I
 Computes the sum of absolute values ([_L1_ norm][l1norm]) of double-precision floating-point strided array elements using pairwise summation and alternative indexing semantics.
 
 ```c
-const double x[] = { 1.0, 2.0, 3.0, 4.0 }
+const double x[] = { 1.0, 2.0, 3.0, 4.0 };
 
 double v = stdlib_strided_dasumpw_ndarray( 4, x, 1, 0 );
 // returns 10.0
diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcusumkbn2/README.md b/lib/node_modules/@stdlib/blas/ext/base/dcusumkbn2/README.md
index 2c80cfad2104..131a16f081b3 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/dcusumkbn2/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/dcusumkbn2/README.md
@@ -199,8 +199,8 @@ console.log( y );
 Computes the cumulative sum of double-precision floating-point strided array elements using a second-order iterative Kahan–Babuška algorithm.
 
 ```c
-const double x[] = { 1.0, 2.0, 3.0, 4.0 }
-double y[] = { 0.0, 0.0, 0.0, 0.0 }
+const double x[] = { 1.0, 2.0, 3.0, 4.0 };
+double y[] = { 0.0, 0.0, 0.0, 0.0 };
 
 stdlib_strided_dcusumkbn2( 4, 0.0, x, 1, y, 1 );
 ```
@@ -227,8 +227,8 @@ void stdlib_strided_dcusumkbn2( const CBLAS_INT N, const double sum, const doubl
 Computes the cumulative sum of double-precision floating-point strided array elements using a second-order iterative Kahan–Babuška algorithm and alternative indexing semantics.
 
 ```c
-const double x[] = { 1.0, 2.0, 3.0, 4.0 }
-double y[] = { 0.0, 0.0, 0.0, 0.0 }
+const double x[] = { 1.0, 2.0, 3.0, 4.0 };
+double y[] = { 0.0, 0.0, 0.0, 0.0 };
 
 stdlib_strided_dcusumkbn2_ndarray( 4, 0.0, x, 1, 0, y, 1, 0 );
 ```
diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcusumors/README.md b/lib/node_modules/@stdlib/blas/ext/base/dcusumors/README.md
index 74ae9a85ad86..61b7fc5618a1 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/dcusumors/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/dcusumors/README.md
@@ -228,8 +228,8 @@ void stdlib_strided_dcusumors( const CBLAS_INT N, double sum, const double *X, c
 Computes the cumulative sum of double-precision floating-point strided array elements using ordinary recursive summation and alternative indexing semantics.
 
 ```c
-const double x[] = { 1.0, 2.0, 3.0, 4.0 }
-double y[] = { 0.0, 0.0, 0.0, 0.0 }
+const double x[] = { 1.0, 2.0, 3.0, 4.0 };
+double y[] = { 0.0, 0.0, 0.0, 0.0 };
 
 stdlib_strided_dcusumors_ndarray( 4, 0.0, x, 1, 0, y, 1, 0 );
 ```
diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcusumpw/README.md b/lib/node_modules/@stdlib/blas/ext/base/dcusumpw/README.md
index a2fddff55051..1b0c9e3cef22 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/dcusumpw/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/dcusumpw/README.md
@@ -228,8 +228,8 @@ void stdlib_strided_dcusumpw( const CBLAS_INT N, const double sum, const double
 Computes the cumulative sum of double-precision floating-point strided array elements using pairwise summation and alternative indexing semantics.
 
 ```c
-const double x[] = { 1.0, 2.0, 3.0, 4.0 }
-double y[] = { 0.0, 0.0, 0.0, 0.0 }
+const double x[] = { 1.0, 2.0, 3.0, 4.0 };
+double y[] = { 0.0, 0.0, 0.0, 0.0 };
 
 stdlib_strided_dcusumpw_ndarray( 4, 0.0, x, 1, 0, y, 1, 0 );
 ```
diff --git a/lib/node_modules/@stdlib/blas/ext/base/dindex-of/README.md b/lib/node_modules/@stdlib/blas/ext/base/dindex-of/README.md
index 31712a75019b..6538aced11c0 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/dindex-of/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/dindex-of/README.md
@@ -222,7 +222,7 @@ Returns the first index of a specified search element in a double-precision floa
 ```c
 double x[] = { 1.0, 2.0, 3.0, 4.0 };
 
-int idx = stdlib_strided_dindex_of( 4, 3.0, x, 1, 0 );
+int idx = stdlib_strided_dindex_of_ndarray( 4, 3.0, x, 1, 0 );
 // returns 2
 ```
 
diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnannsumors/README.md b/lib/node_modules/@stdlib/blas/ext/base/dnannsumors/README.md
index 65de33ac5618..e1952e1bd894 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/dnannsumors/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/dnannsumors/README.md
@@ -259,7 +259,7 @@ double stdlib_strided_dnannsumors_ndarray( const CBLAS_INT N, const double *X, c
 
 ```c
 #include "stdlib/blas/ext/base/dnannsumors.h"
-#include "stdlib/blase/base/shared.h"
+#include "stdlib/blas/base/shared.h"
 #include <stdio.h>
 
 int main( void ) {
diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnannsumpw/README.md b/lib/node_modules/@stdlib/blas/ext/base/dnannsumpw/README.md
index 9046d473826b..01d48459a332 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/dnannsumpw/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/dnannsumpw/README.md
@@ -262,7 +262,7 @@ double stdlib_strided_dnannsumpw_ndarray( const CBLAS_INT N, const double *X, co
 
 ```c
 #include "stdlib/blas/ext/base/dnannsumpw.h"
-#include "stdlib/blase/base/shared.h"
+#include "stdlib/blas/base/shared.h"
 #include <stdio.h>
 
 int main( void ) {
diff --git a/lib/node_modules/@stdlib/blas/ext/base/drrss/README.md b/lib/node_modules/@stdlib/blas/ext/base/drrss/README.md
index 1382897c81c8..660d891429a2 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/drrss/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/drrss/README.md
@@ -238,7 +238,7 @@ Computes the square root of the [residual sum of squares][wikipedia-residual-sum
 const double x[] = { 1.0, -2.0, 2.0 };
 const double y[] = { 1.0, 1.0, -4.0 };
 
-double v = stdlib_strided_drrss_ndarray( 3, x, 1, 0, 1, 0 );
+double v = stdlib_strided_drrss_ndarray( 3, x, 1, 0, y, 1, 0 );
 // returns ~6.7
 ```
 
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sindex-of/README.md b/lib/node_modules/@stdlib/blas/ext/base/sindex-of/README.md
index f2cb92b0e5c1..dc945842cf96 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sindex-of/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/sindex-of/README.md
@@ -222,7 +222,7 @@ Returns the first index of a specified search element in a single-precision floa
 ```c
 float x[] = { 1.0f, 2.0f, 3.0f, 4.0f };
 
-int idx = stdlib_strided_sindex_of( 4, 3.0f, x, 1, 0 );
+int idx = stdlib_strided_sindex_of_ndarray( 4, 3.0f, x, 1, 0 );
 // returns 2
 ```
 
diff --git a/lib/node_modules/@stdlib/complex/float32/base/mul-add/README.md b/lib/node_modules/@stdlib/complex/float32/base/mul-add/README.md
index ac0da3f66d3c..274afb7653af 100644
--- a/lib/node_modules/@stdlib/complex/float32/base/mul-add/README.md
+++ b/lib/node_modules/@stdlib/complex/float32/base/mul-add/README.md
@@ -183,9 +183,9 @@ Performs a multiply-add operation involving three single-precision complex float
 #include "stdlib/complex/float32/real.h"
 #include "stdlib/complex/float32/imag.h"
 
-stdlib_complex64_t z1 = sstdlib_complex64( 5.0f, 3.0f );
-stdlib_complex64_t z2 = sstdlib_complex64( -2.0f, 1.0f );
-stdlib_complex64_t z3 = sstdlib_complex64( 7.0f, -8.0f );
+stdlib_complex64_t z1 = stdlib_complex64( 5.0f, 3.0f );
+stdlib_complex64_t z2 = stdlib_complex64( -2.0f, 1.0f );
+stdlib_complex64_t z3 = stdlib_complex64( 7.0f, -8.0f );
 
 stdlib_complex64_t out = stdlib_base_complex64_muladd( z1, z2, z3 );
 
diff --git a/lib/node_modules/@stdlib/math/base/special/acotd/README.md b/lib/node_modules/@stdlib/math/base/special/acotd/README.md
index 829092b881e1..4760db9ef127 100644
--- a/lib/node_modules/@stdlib/math/base/special/acotd/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/acotd/README.md
@@ -114,7 +114,7 @@ double out = stdlib_base_acotd( 0.0 );
 // returns 90.0
 
 out = stdlib_base_acotd( 3.0 );
-// returns ~30.0
+// returns ~18.43495
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/math/base/special/acotdf/README.md b/lib/node_modules/@stdlib/math/base/special/acotdf/README.md
index e0ea3bd35495..2f76954bfb03 100644
--- a/lib/node_modules/@stdlib/math/base/special/acotdf/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/acotdf/README.md
@@ -114,7 +114,7 @@ float out = stdlib_base_acotdf( 0.0f );
 // returns 90.0f
 
 out = stdlib_base_acotdf( 3.0f );
-// returns ~30.0f
+// returns 18.43495f
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/math/base/special/ceiln/README.md b/lib/node_modules/@stdlib/math/base/special/ceiln/README.md
index 1e92fe9ad0f4..2f5a5a0cdce9 100644
--- a/lib/node_modules/@stdlib/math/base/special/ceiln/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/ceiln/README.md
@@ -134,11 +134,11 @@ double y = stdlib_base_ceiln( 3.141592653589793, -2 );
 // returns 3.15
 
 // If n = 0, `ceiln` behaves like `ceil`:
-double y = stdlib_base_ceiln( 3.141592653589793, 0 );
+y = stdlib_base_ceiln( 3.141592653589793, 0 );
 // returns 4.0
 
 // Round a value to the nearest thousand:
-double y = stdlib_base_ceiln( 12368.0, 3 );
+y = stdlib_base_ceiln( 12368.0, 3 );
 // returns 13000.0
 ```
 
diff --git a/lib/node_modules/@stdlib/math/base/special/cflipsign/README.md b/lib/node_modules/@stdlib/math/base/special/cflipsign/README.md
index 65fe10322c4f..8fd0610d4ddb 100644
--- a/lib/node_modules/@stdlib/math/base/special/cflipsign/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/cflipsign/README.md
@@ -193,6 +193,7 @@ int main( void ) {
         stdlib_complex128_reim( v, &re1, &im1 );
         stdlib_complex128_reim( y, &re2, &im2 );
         printf( "cflipsign(%lf + %lfi, %lf) = %lf + %lfi\n", re1, im1, -1.0, re2, im2 );
+    };
 }
 ```
 
diff --git a/lib/node_modules/@stdlib/math/base/special/copysignf/README.md b/lib/node_modules/@stdlib/math/base/special/copysignf/README.md
index e623fd31548f..df90de263420 100644
--- a/lib/node_modules/@stdlib/math/base/special/copysignf/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/copysignf/README.md
@@ -124,7 +124,7 @@ Returns a [single-precision floating-point number][ieee754] with the magnitude o
 float v = stdlib_base_copysignf( -3.0f, 10.0f );
 // returns 3.0f
 
-float v = stdlib_base_copysignf( 3.0f, -1.0f );
+v = stdlib_base_copysignf( 3.0f, -1.0f );
 // returns -3.0f
 ```
 
diff --git a/lib/node_modules/@stdlib/math/base/special/cphase/README.md b/lib/node_modules/@stdlib/math/base/special/cphase/README.md
index c5f4989c9645..47b83a0e95ec 100644
--- a/lib/node_modules/@stdlib/math/base/special/cphase/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/cphase/README.md
@@ -158,17 +158,14 @@ int main( void ) {
     };
 
     stdlib_complex128_t v;
-    stdlib_complex128_t y;
-    double re1;
-    double im1;
-    double re2;
-    double im2;
+    double y;
+    double re;
+    double im;
     int i;
     for ( i = 0; i < 4; i++ ) {
         v = x[ i ];
         y = stdlib_base_cphase( v );
-        stdlib_complex128_reim( v, &re1, &im1 );
-        stdlib_complex128_reim( y, &re2, &im2 );
+        stdlib_complex128_reim( v, &re, &im );
         printf( "cphase(%lf + %lfi) = %lf\n", re, im, y );
     }
 }
diff --git a/lib/node_modules/@stdlib/math/base/special/flipsignf/README.md b/lib/node_modules/@stdlib/math/base/special/flipsignf/README.md
index f50febe785fd..b02f15b47f81 100644
--- a/lib/node_modules/@stdlib/math/base/special/flipsignf/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/flipsignf/README.md
@@ -127,7 +127,7 @@ Returns a [single-precision floating-point number][ieee754] with the magnitude o
 float v = stdlib_base_flipsignf( -3.0f, 10.0f );
 // returns -3.0f
 
-float v = stdlib_base_flipsignf( -3.0f, -1.0f );
+v = stdlib_base_flipsignf( -3.0f, -1.0f );
 // returns 3.0f
 ```
 
diff --git a/lib/node_modules/@stdlib/math/base/special/fmod/README.md b/lib/node_modules/@stdlib/math/base/special/fmod/README.md
index 5af3face8720..65dc6e4a5375 100644
--- a/lib/node_modules/@stdlib/math/base/special/fmod/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/fmod/README.md
@@ -132,7 +132,7 @@ Evaluates the modulus function.
 
 ```c
 double out = stdlib_base_fmod( 8.9, 3.0 );
-// returns 2.9
+// returns ~2.9
 
 out = stdlib_base_fmod( 4.0, 2.0 );
 // returns 0.0
diff --git a/lib/node_modules/@stdlib/math/base/special/negalucas/README.md b/lib/node_modules/@stdlib/math/base/special/negalucas/README.md
index 863cdb801dc2..2e6bb0fb637e 100644
--- a/lib/node_modules/@stdlib/math/base/special/negalucas/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/negalucas/README.md
@@ -192,7 +192,7 @@ Computes the nth [negaLucas number][lucas-number].
 
 ```c
 double out = stdlib_base_negalucas( 0.0 );
-// returns 0.0
+// returns 2.0
 
 out = stdlib_base_negalucas( -1.0 );
 // returns -1.0
diff --git a/lib/node_modules/@stdlib/math/base/special/negalucasf/README.md b/lib/node_modules/@stdlib/math/base/special/negalucasf/README.md
index c2e806115e4f..58f4d800e60c 100644
--- a/lib/node_modules/@stdlib/math/base/special/negalucasf/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/negalucasf/README.md
@@ -192,7 +192,7 @@ Computes the nth [negaLucas number][lucas-number] in single-precision floating-p
 
 ```c
 float out = stdlib_base_negalucasf( 0.0f );
-// returns 0.0f
+// returns 2.0f
 
 out = stdlib_base_negalucasf( -1.0f );
 // returns -1.0f
diff --git a/lib/node_modules/@stdlib/math/base/special/sec/README.md b/lib/node_modules/@stdlib/math/base/special/sec/README.md
index 7c87b1aed98d..8c1b685664a2 100644
--- a/lib/node_modules/@stdlib/math/base/special/sec/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/sec/README.md
@@ -119,8 +119,8 @@ Evaluates the [secant][trigonometric-functions] of `x` (in radians).
 double out = stdlib_base_sec( 0.0 );
 // returns 1.0
 
-out = stdlib_base_cos( 3.141592653589793 / 2.0 );
-// returns 6.123233995736766e-17
+out = stdlib_base_sec( 3.141592653589793 / 2.0 );
+// returns 16331239353195370.0
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/math/base/special/vercosf/README.md b/lib/node_modules/@stdlib/math/base/special/vercosf/README.md
index 0faba21cb70e..33f3c917f452 100644
--- a/lib/node_modules/@stdlib/math/base/special/vercosf/README.md
+++ b/lib/node_modules/@stdlib/math/base/special/vercosf/README.md
@@ -124,7 +124,7 @@ float out = stdlib_base_vercosf( 0.0f );
 // returns 2.0f
 
 out = stdlib_base_vercosf( 3.141592653589793f / 2.0f );
-// returns 1.0f
+// returns ~1.0f
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/ndarray/base/count-falsy/README.md b/lib/node_modules/@stdlib/ndarray/base/count-falsy/README.md
index 5d2cc147f6e2..6b2080661510 100644
--- a/lib/node_modules/@stdlib/ndarray/base/count-falsy/README.md
+++ b/lib/node_modules/@stdlib/ndarray/base/count-falsy/README.md
@@ -2198,7 +2198,7 @@ int main( void ) {
     struct ndarray *arrays[] = { x, y };
 
     // Perform operation:
-    int8_t status = stdlib_ndarray_count_falsy_b_i( arrays );
+    int8_t status = stdlib_ndarray_count_falsy_b_i( arrays, NULL );
     if ( status != 0 ) {
         fprintf( stderr, "Error during computation.\n" );
         exit( EXIT_FAILURE );
diff --git a/lib/node_modules/@stdlib/ndarray/base/count-truthy/README.md b/lib/node_modules/@stdlib/ndarray/base/count-truthy/README.md
index 549f5fc9f091..e825171c1f2b 100644
--- a/lib/node_modules/@stdlib/ndarray/base/count-truthy/README.md
+++ b/lib/node_modules/@stdlib/ndarray/base/count-truthy/README.md
@@ -2198,7 +2198,7 @@ int main( void ) {
     struct ndarray *arrays[] = { x, y };
 
     // Perform operation:
-    int8_t status = stdlib_ndarray_count_truthy_b_i( arrays );
+    int8_t status = stdlib_ndarray_count_truthy_b_i( arrays, NULL );
     if ( status != 0 ) {
         fprintf( stderr, "Error during computation.\n" );
         exit( EXIT_FAILURE );
diff --git a/lib/node_modules/@stdlib/ndarray/base/flatten-shape-from/README.md b/lib/node_modules/@stdlib/ndarray/base/flatten-shape-from/README.md
index 9a91d06d458b..7fe1cc131e00 100644
--- a/lib/node_modules/@stdlib/ndarray/base/flatten-shape-from/README.md
+++ b/lib/node_modules/@stdlib/ndarray/base/flatten-shape-from/README.md
@@ -153,7 +153,7 @@ const int64_t ndims = 3;
 const int64_t shape[] = { 2, 3, 10 };
 int64_t out[ 2 ];
 
-stdlib_ndarray_flatten_shape( ndims, shape, 1, out );
+stdlib_ndarray_flatten_shape_from( ndims, shape, 1, out );
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/number/float32/base/signbit/README.md b/lib/node_modules/@stdlib/number/float32/base/signbit/README.md
index 74224616c756..aa1de3440726 100644
--- a/lib/node_modules/@stdlib/number/float32/base/signbit/README.md
+++ b/lib/node_modules/@stdlib/number/float32/base/signbit/README.md
@@ -156,7 +156,7 @@ int main( void ) {
     int8_t out;
     int i;
     for ( i = 0; i < 9; i++ ) {
-        stdlib_base_float32_signbit( x[ i ], &out );
+        out = stdlib_base_float32_signbit( x[ i ] );
         printf( "%f => signbit: %" PRId8 "\n", x[ i ], out );
     }
 }
diff --git a/lib/node_modules/@stdlib/number/float64/base/signbit/README.md b/lib/node_modules/@stdlib/number/float64/base/signbit/README.md
index 916aae7ec6c9..8be728bd60ea 100644
--- a/lib/node_modules/@stdlib/number/float64/base/signbit/README.md
+++ b/lib/node_modules/@stdlib/number/float64/base/signbit/README.md
@@ -152,7 +152,7 @@ int main( void ) {
     int8_t out;
     int i;
     for ( i = 0; i < 9; i++ ) {
-        stdlib_base_float64_signbit( x[ i ], &out );
+        out = stdlib_base_float64_signbit( x[ i ] );
         printf( "%lf => signbit: %" PRId8 "\n", x[ i ], out );
     }
 }
diff --git a/lib/node_modules/@stdlib/stats/base/dists/arcsine/pdf/README.md b/lib/node_modules/@stdlib/stats/base/dists/arcsine/pdf/README.md
index 600b61c2d6b1..646ef97fc0a3 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/arcsine/pdf/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/arcsine/pdf/README.md
@@ -164,7 +164,7 @@ logEachMap( 'x: %0.4f, a: %0.4f, b: %0.4f, f(x;a,b): %0.4f', x, a, b, pdf );
 Evaluates the probability density function (PDF) for an arcsine distribution.
 
 ```c
-double out = stdlib_base_dists_arcsine_pdf( 0.5, 0.0, 4.0 );
+double out = stdlib_base_dists_arcsine_pdf( 2, 0.0, 4.0 );
 // returns ~0.159
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/arcsine/quantile/README.md b/lib/node_modules/@stdlib/stats/base/dists/arcsine/quantile/README.md
index 3d04a1bf8c9d..d10b493ed4d6 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/arcsine/quantile/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/arcsine/quantile/README.md
@@ -166,7 +166,7 @@ logEachMap( 'p: %0.4f, a: %0.4f, b: %0.4f, Q(p;a,b): %0.4f', p, a, b, quantile )
 Evaluates the quantile function for an arcsine distribution.
 
 ```c
-double out = quantile( 0.8, 0.0, 1.0 );
+double out = stdlib_base_dists_arcsine_quantile( 0.8, 0.0, 1.0 );
 // returns ~0.905
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/cauchy/pdf/README.md b/lib/node_modules/@stdlib/stats/base/dists/cauchy/pdf/README.md
index ec7fb2b5e34c..bbeda0ea4399 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/cauchy/pdf/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/cauchy/pdf/README.md
@@ -167,7 +167,7 @@ Evaluates the [probability density function][pdf] (PDF) for a [Cauchy][cauchy-di
 
 ```c
 double out = stdlib_base_dists_cauchy_pdf( 0.5, 0.0, 2.0 );
-// returns ~0.333
+// returns ~0.14979
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/stats/base/dists/chi/logpdf/README.md b/lib/node_modules/@stdlib/stats/base/dists/chi/logpdf/README.md
index b2c58a2e1b24..6afbe6235ab9 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/chi/logpdf/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/chi/logpdf/README.md
@@ -169,7 +169,7 @@ Evaluates the natural logarithm of the [probability density function][pdf] (PDF)
 
 ```c
 double out = stdlib_base_dists_chi_logpdf( 2.0, 2.0 );
-// returns ~-1.309
+// returns ~-1.307
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/stats/base/dists/frechet/cdf/README.md b/lib/node_modules/@stdlib/stats/base/dists/frechet/cdf/README.md
index 80b0e6915112..6215b936309f 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/frechet/cdf/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/frechet/cdf/README.md
@@ -192,7 +192,7 @@ for ( i = 0; i < 100; i++ ) {
 Evaluates the cumulative distribution function (CDF) for a Fréchet distribution with shape `alpha`, scale `s`, and location `m` at a value `x`.
 
 ```c
-double y = stdlib_base_frechet_cdf( 10.0, 2.0, 3.0, 2.0 );
+double y = stdlib_base_dists_frechet_cdf( 10.0, 2.0, 3.0, 2.0 );
 // returns ~0.869
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/gumbel/skewness/README.md b/lib/node_modules/@stdlib/stats/base/dists/gumbel/skewness/README.md
index 1e4afc5caec0..9811d0d49577 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/gumbel/skewness/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/gumbel/skewness/README.md
@@ -213,7 +213,7 @@ int main( void ) {
         beta = random_uniform( 0.0, 20.0 );
         y = stdlib_base_dists_gumbel_skewness( mu, beta );
         printf( "µ: %lf, β: %lf, Skew(X;µ,β): %lf\n", mu, beta, y );
-    s}
+    }
 }
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/levy/mean/README.md b/lib/node_modules/@stdlib/stats/base/dists/levy/mean/README.md
index 33eae0f143f5..94e3f51ea147 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/levy/mean/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/levy/mean/README.md
@@ -210,8 +210,8 @@ int main( void ) {
     int i;
 
     for ( i = 0; i < 25; i++ ) {
-        x = random_uniform( -5.0, 5.0 );
-        a = random_uniform( 0.0, 20.0 );
+        mu = random_uniform( -5.0, 5.0 );
+        c = random_uniform( 0.0, 20.0 );
         y = stdlib_base_dists_levy_mean( mu, c );
         printf( "µ: %lf, c: %lf, E(X;µ,c): %lf\n", mu, c, y );
     }
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/logcdf/README.md b/lib/node_modules/@stdlib/stats/base/dists/normal/logcdf/README.md
index 6299508c9150..7b1993779b45 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/logcdf/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/logcdf/README.md
@@ -155,7 +155,7 @@ logEachMap( 'x: %0.4f, µ: %0.4f, σ: %0.4f, ln(F(x;µ,σ)): %0.4f', x, mu, sigm
 Evaluates the natural logarithm of the [cumulative distribution function][cdf] (CDF) for a [normal][normal-distribution] distribution with parameters `mu` (mean) and `sigma` (standard deviation).
 
 ```c
-double out = stdlib_base_dists_normal_cdf( 2.0, 0.0, 1.0 );
+double out = stdlib_base_dists_normal_logcdf( 2.0, 0.0, 1.0 );
 // returns ~-0.023
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/logpdf/README.md b/lib/node_modules/@stdlib/stats/base/dists/normal/logpdf/README.md
index 9af71dfb742e..f5d25ce34bdf 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/logpdf/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/logpdf/README.md
@@ -169,7 +169,7 @@ logEachMap( 'x: %0.4f, µ: %0.4f, σ: %0.4f, ln(f(x;µ,σ)): %0.4f', x, mu, sigm
 Evaluates the natural logarithm of the probability density function (logPDF) for a normal distribution.
 
 ```c
-double out = stdlib_base_dists_normal_pdf( 2.0, 0.0, 1.0 );
+double out = stdlib_base_dists_normal_logpdf( 2.0, 0.0, 1.0 );
 // returns ~-2.919
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/poisson/logpmf/README.md b/lib/node_modules/@stdlib/stats/base/dists/poisson/logpmf/README.md
index feb0d40e49ba..60318cdda313 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/poisson/logpmf/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/poisson/logpmf/README.md
@@ -204,7 +204,7 @@ double stdlib_base_dists_poisson_logpmf( const double x, const double lambda );
 
 ```c
 #include "stdlib/stats/base/dists/poisson/logpmf.h"
-#include "stdlib/random/base/ceil.h"
+#include "stdlib/math/base/special/ceil.h"
 #include <stdlib.h>
 #include <stdio.h>
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/poisson/mode/README.md b/lib/node_modules/@stdlib/stats/base/dists/poisson/mode/README.md
index ac4482aad670..07f3a2892afb 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/poisson/mode/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/poisson/mode/README.md
@@ -144,7 +144,7 @@ logEachMap( 'λ: %0.4f, mode(X;λ): %0.4f', lambda, mode );
 Returns the mode of a Poisson distribution.
 
 ```c
-double out = sstdlib_base_dists_poisson_mode( 9.0 );
+double out = stdlib_base_dists_poisson_mode( 9.0 );
 // returns 9
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/poisson/skewness/README.md b/lib/node_modules/@stdlib/stats/base/dists/poisson/skewness/README.md
index 72c66fbccf81..eff87c2596be 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/poisson/skewness/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/poisson/skewness/README.md
@@ -144,7 +144,7 @@ logEachMap( 'λ: %0.4f, skew(X;λ): %0.4f', lambda, skewness );
 Returns the skewness of a Poisson distribution.
 
 ```c
-double out = sstdlib_base_dists_poisson_skewness( 9.0 );
+double out = stdlib_base_dists_poisson_skewness( 9.0 );
 // returns ~0.33
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/t/entropy/README.md b/lib/node_modules/@stdlib/stats/base/dists/t/entropy/README.md
index 940a5f4cda1c..914834f0b339 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/t/entropy/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/t/entropy/README.md
@@ -148,7 +148,7 @@ Returns the differential entropy of a Student's t distribution.
 
 ```c
 double out = stdlib_base_dists_t_entropy( 9.0 );
-// returns ~1.553
+// returns ~1.533
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/stats/strided/dmeanvarpn/README.md b/lib/node_modules/@stdlib/stats/strided/dmeanvarpn/README.md
index 0d09eb408095..676f1d46c23f 100644
--- a/lib/node_modules/@stdlib/stats/strided/dmeanvarpn/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/dmeanvarpn/README.md
@@ -258,7 +258,7 @@ Computes the [mean][arithmetic-mean] and [variance][variance] of a double-precis
 
 ```c
 const double x[] = { 1.0, -2.0, 2.0 };
-double out[] = { 0.0, 0.0 }
+double out[] = { 0.0, 0.0 };
 
 stdlib_strided_dmeanvarpn( 3, 1.0, x, 1, out, 1 );
 ```
@@ -282,7 +282,7 @@ Computes the [mean][arithmetic-mean] and [variance][variance] of a double-precis
 
 ```c
 const double x[] = { 1.0, -2.0, 2.0 };
-double out[] = { 0.0, 0.0 }
+double out[] = { 0.0, 0.0 };
 
 stdlib_strided_dmeanvarpn_ndarray( 3, 1.0, x, 1, 0, out, 1, 0 );
 ```
diff --git a/lib/node_modules/@stdlib/stats/strided/dmskmin/README.md b/lib/node_modules/@stdlib/stats/strided/dmskmin/README.md
index fbda72df9aff..66e7c4ec9f51 100644
--- a/lib/node_modules/@stdlib/stats/strided/dmskmin/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/dmskmin/README.md
@@ -232,7 +232,7 @@ Computes the minimum value of a double-precision floating-point strided array ac
 const double x[] = { 1.0, -2.0, 2.0 };
 const uint8_t mask[] = { 0, 1, 0 };
 
-double v = stdlib_strided_dmskmin( 3, x, 1, 0, mask, 1, 0 );
+double v = stdlib_strided_dmskmin_ndarray( 3, x, 1, 0, mask, 1, 0 );
 // returns 1.0
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/strided/dmskrange/README.md b/lib/node_modules/@stdlib/stats/strided/dmskrange/README.md
index a0717741a784..c4335b39d249 100644
--- a/lib/node_modules/@stdlib/stats/strided/dmskrange/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/dmskrange/README.md
@@ -233,7 +233,7 @@ Computes the [range][range] of a double-precision floating-point strided array a
 const double x[] = { 1.0, -2.0, 2.0 };
 const uint8_t mask[] = { 0, 1, 0 };
 
-double v = stdlib_strided_dmskrange( 3, x, 1, 0, mask, 1, 0 );
+double v = stdlib_strided_dmskrange_ndarray( 3, x, 1, 0, mask, 1, 0 );
 // returns 1.0
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/strided/dnanmaxabs/README.md b/lib/node_modules/@stdlib/stats/strided/dnanmaxabs/README.md
index f16b9bb81c3f..5d0f47127aad 100644
--- a/lib/node_modules/@stdlib/stats/strided/dnanmaxabs/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/dnanmaxabs/README.md
@@ -178,17 +178,17 @@ console.log( v );
 ### Usage
 
 ```c
-#include "stdlib/stats/strided/dnanmax.h"
+#include "stdlib/stats/strided/dnanmaxabs.h"
 ```
 
-#### stdlib_strided_dnanmax( N, \*X, strideX )
+#### stdlib_strided_dnanmaxabs( N, \*X, strideX )
 
 Computes the maximum absolute value of a double-precision floating-point strided array, ignoring `NaN` values.
 
 ```c
 const double x[] = { 1.0, -2.0, 0.0 / 0.0, -4.0 };
 
-double v = stdlib_strided_dnanmax( 4, x, 1 );
+double v = stdlib_strided_dnanmaxabs( 4, x, 1 );
 // returns 4.0
 ```
 
@@ -199,17 +199,17 @@ The function accepts the following arguments:
 -   **strideX**: `[in] CBLAS_INT` stride length for `X`.
 
 ```c
-double stdlib_strided_dnanmax( const CBLAS_INT N, const double *X, const CBLAS_INT strideX );
+double stdlib_strided_dnanmaxabs( const CBLAS_INT N, const double *X, const CBLAS_INT strideX );
 ```
 
-#### stdlib_strided_dnanmax_ndarray( N, \*X, strideX, offsetX )
+#### stdlib_strided_dnanmaxabs_ndarray( N, \*X, strideX, offsetX )
 
 Computes the maximum absolute value of a double-precision floating-point strided array, ignoring `NaN` values and using alternative indexing semantics.
 
 ```c
 const double x[] = { 1.0, -2.0, 0.0 / 0.0, -4.0 };
 
-double v = stdlib_strided_dnanmax_ndarray( 4, x, 1, 0 );
+double v = stdlib_strided_dnanmaxabs_ndarray( 4, x, 1, 0 );
 // returns 4.0
 ```
 
@@ -221,7 +221,7 @@ The function accepts the following arguments:
 -   **offsetX**: `[in] CBLAS_INT` starting index for `X`.
 
 ```c
-double stdlib_strided_dnanmax_ndarray( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX );
+double stdlib_strided_dnanmaxabs_ndarray( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX );
 ```
 
 </section>
diff --git a/lib/node_modules/@stdlib/stats/strided/dnanmeanpn/README.md b/lib/node_modules/@stdlib/stats/strided/dnanmeanpn/README.md
index 3688ad47ab3a..997197290e7d 100644
--- a/lib/node_modules/@stdlib/stats/strided/dnanmeanpn/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/dnanmeanpn/README.md
@@ -229,7 +229,7 @@ Computes the aithmetic mean of a double-precision floating-point strided array,
 ```c
 const double x[] = { 1.0, 2.0, 0.0/0.0, 3.0, 0.0/0.0, 4.0, 5.0, 6.0, 0.0/0.0, 7.0, 8.0, 0.0/0.0 };
 
-double v = stdlib_strided_dnanmeanpn( 6, x, 2, 0 );
+double v = stdlib_strided_dnanmeanpn_ndarray( 6, x, 2, 0 );
 // returns ~4.6667
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/strided/dnanmin/README.md b/lib/node_modules/@stdlib/stats/strided/dnanmin/README.md
index ac9ad0041ad7..9c52529599ba 100644
--- a/lib/node_modules/@stdlib/stats/strided/dnanmin/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/dnanmin/README.md
@@ -189,7 +189,7 @@ Calculate the minimum value of a double-precision floating-point strided array,
 const double x[] = { 1.0, -2.0, 0.0 / 0.0, -4.0 };
 
 double v = stdlib_strided_dnanmin( 4, x, 1 );
-// returns 1.0
+// returns -4.0
 ```
 
 The function accepts the following arguments:
@@ -210,7 +210,7 @@ Computes the minimum value of a double-precision floating-point strided array, i
 const double x[] = { 1.0, -2.0, 0.0 / 0.0, -4.0 };
 
 double v = stdlib_strided_dnanmin_ndarray( 4, x, 1, 0 );
-// returns 1.0
+// returns -4.0
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/stats/strided/drange/README.md b/lib/node_modules/@stdlib/stats/strided/drange/README.md
index 0cc69d537c4b..df9e57465e09 100644
--- a/lib/node_modules/@stdlib/stats/strided/drange/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/drange/README.md
@@ -182,7 +182,7 @@ Computes the [range][range] of a double-precision floating-point strided array `
 const double x[] = { 1.0, -2.0, 3.0, -4.0 };
 
 double v = stdlib_strided_drange( 4, x, 1 );
-// returns 4.0
+// returns 7.0
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/stats/strided/dstdev/README.md b/lib/node_modules/@stdlib/stats/strided/dstdev/README.md
index f6334c41457b..fdf0d321d052 100644
--- a/lib/node_modules/@stdlib/stats/strided/dstdev/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/dstdev/README.md
@@ -243,7 +243,7 @@ Computes the [standard deviation][standard-deviation] of a double-precision floa
 const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 };
 
 double v = stdlib_strided_dstdev( 4, 1.0, x, 2 );
-// returns 2.581989
+// returns ~2.581989
 ```
 
 The function accepts the following arguments:
@@ -265,7 +265,7 @@ Computes the [standard deviation][standard-deviation] of a double-precision floa
 const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 };
 
 double v = stdlib_strided_dstdev_ndarray( 4, 1.0, x, 2, 0 );
-// returns 2.581989
+// returns ~2.581989
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/stats/strided/dsvariance/README.md b/lib/node_modules/@stdlib/stats/strided/dsvariance/README.md
index 8e98a1900c9b..ed609ccfc1bd 100644
--- a/lib/node_modules/@stdlib/stats/strided/dsvariance/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/dsvariance/README.md
@@ -241,7 +241,7 @@ console.log( v );
 Computes the [variance][variance] of a single-precision floating-point strided array using extended accumulation and returning an extended precision result.
 
 ```c
-const float x[] = { 1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f, 8.0f }
+const float x[] = { 1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f, 8.0f };
 
 double v = stdlib_strided_dsvariance( 8, 1.0f, x, 1 );
 // returns 6.0
@@ -263,7 +263,7 @@ double stdlib_strided_dsvariance( const CBLAS_INT N, const float correction, con
 Computes the [variance][variance] of a single-precision floating-point strided array using extended accumulation and alternative indexing semantics and returning an extended precision result.
 
 ```c
-const float x[] = { 1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f, 8.0f }
+const float x[] = { 1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f, 8.0f };
 
 double v = stdlib_strided_dsvariance_ndarray( 4, 1.0f, x, 2, 0 );
 // returns ~6.666667
diff --git a/lib/node_modules/@stdlib/stats/strided/dvariance/README.md b/lib/node_modules/@stdlib/stats/strided/dvariance/README.md
index ddea8cdca0dc..eaaa545c0a77 100644
--- a/lib/node_modules/@stdlib/stats/strided/dvariance/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/dvariance/README.md
@@ -240,7 +240,7 @@ console.log( v );
 Computes the [variance][variance] of a double-precision floating-point strided array.
 
 ```c
-const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 }
+const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 };
 
 double v = stdlib_strided_dvariance( 8, 1.0, x, 1 );
 // returns 6.0
@@ -262,7 +262,7 @@ double stdlib_strided_dvariance( const CBLAS_INT N, const double correction, con
 Computes the [variance][variance] of a double-precision floating-point strided array using alternative indexing semantics.
 
 ```c
-const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 }
+const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 };
 
 double v = stdlib_strided_dvariance_ndarray( 4, 1.0, x, 2, 0 );
 // returns ~6.666667
diff --git a/lib/node_modules/@stdlib/stats/strided/dvariancepn/README.md b/lib/node_modules/@stdlib/stats/strided/dvariancepn/README.md
index 3f06234afb99..a826d6419ebc 100644
--- a/lib/node_modules/@stdlib/stats/strided/dvariancepn/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/dvariancepn/README.md
@@ -240,7 +240,7 @@ console.log( v );
 Computes the [variance][variance] of a double-precision floating-point strided array using a two-pass algorithm.
 
 ```c
-const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 }
+const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 };
 
 double v = stdlib_strided_dvariancepn( 8, 1.0, x, 1 );
 // returns 6.0
@@ -262,7 +262,7 @@ double stdlib_strided_dvariancepn( const CBLAS_INT N, const double correction, c
 Computes the [variance][variance] of a double-precision floating-point strided array using a two-pass algorithm and alternative indexing semantics.
 
 ```c
-const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 }
+const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 };
 
 double v = stdlib_strided_dvariancepn_ndarray( 4, 1.0, x, 2, 0 );
 // returns ~6.666667
diff --git a/lib/node_modules/@stdlib/stats/strided/dvarmpn/README.md b/lib/node_modules/@stdlib/stats/strided/dvarmpn/README.md
index 0810b2b28309..20e39b3f26e4 100644
--- a/lib/node_modules/@stdlib/stats/strided/dvarmpn/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/dvarmpn/README.md
@@ -267,7 +267,7 @@ Computes the [variance][variance] of a double-precision floating-point strided a
 const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 };
 
 double v = stdlib_strided_dvarmpn_ndarray( 4, 1.0, 1.25, x, 2, 0 );
-// returns 6.25
+// returns ~6.67
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/stats/strided/sdsnanmeanors/README.md b/lib/node_modules/@stdlib/stats/strided/sdsnanmeanors/README.md
index 948b3b977a6b..bcfc2ae99f9d 100644
--- a/lib/node_modules/@stdlib/stats/strided/sdsnanmeanors/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/sdsnanmeanors/README.md
@@ -202,7 +202,7 @@ Computes the [arithmetic mean][arithmetic-mean] of a single-precision floating-p
 ```c
 const float x[] = { 1.0f, -2.0f, 0.0f/0.0f, 2.0f };
 
-float v = stdlib_strided_sdsnanmeanors( x.length, x, 1 );
+float v = stdlib_strided_sdsnanmeanors( 4, x, 1 );
 // returns ~0.3333f
 ```
 
@@ -223,7 +223,7 @@ Computes the [arithmetic mean][arithmetic-mean] of a single-precision floating-p
 ```c
 const float x[] = { 1.0f, -2.0f, 0.0f/0.0f, 2.0f };
 
-float v = stdlib_strided_sdsnanmeanors_ndarray( x.length, x, 1, 0 );
+float v = stdlib_strided_sdsnanmeanors_ndarray( 4, x, 1, 0 );
 // returns ~0.3333f
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/strided/smidrange/README.md b/lib/node_modules/@stdlib/stats/strided/smidrange/README.md
index 66e1723c5be8..5d1990309829 100644
--- a/lib/node_modules/@stdlib/stats/strided/smidrange/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/smidrange/README.md
@@ -181,7 +181,7 @@ Computes the [mid-range][mid-range] of a single-precision floating-point strided
 const float x[] = { 1.0f, -2.0f, 3.0f, 4.0f };
 
 float v = stdlib_strided_smidrange( 4, x, 1 );
-// returns 2.5f
+// returns 1.0f
 ```
 
 The function accepts the following arguments:
@@ -202,7 +202,7 @@ Computes the [mid-range][mid-range] of a single-precision floating-point strided
 const float x[] = { 1.0f, -2.0f, 3.0f, 4.0f };
 
 float v = stdlib_strided_smidrange_ndarray( 4, x, 1, 0 );
-// returns 2.5f
+// returns 1.0f
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/stats/strided/smin/README.md b/lib/node_modules/@stdlib/stats/strided/smin/README.md
index f4714b8c2b3d..4cfeb4819580 100644
--- a/lib/node_modules/@stdlib/stats/strided/smin/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/smin/README.md
@@ -179,7 +179,7 @@ Computes the minimum value of a single-precision floating-point strided array.
 const float x[] = { 1.0f, -2.0f, 3.0f, -4.0f };
 
 float v = stdlib_strided_smin( 4, x, 1 );
-// returns 1.0f
+// returns -4.0f
 ```
 
 The function accepts the following arguments:
@@ -200,7 +200,7 @@ Computes the minimum value of a single-precision floating-point strided array us
 const float x[] = { 1.0f, -2.0f, 3.0f, -4.0f };
 
 float v = stdlib_strided_smin_ndarray( 4, x, 1, 0 );
-// returns 1.0f
+// returns -4.0f
 ```
 
 The function accepts the following arguments:
diff --git a/lib/node_modules/@stdlib/stats/strided/sminabs/README.md b/lib/node_modules/@stdlib/stats/strided/sminabs/README.md
index 39c864b17899..8a530b922232 100644
--- a/lib/node_modules/@stdlib/stats/strided/sminabs/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/sminabs/README.md
@@ -199,7 +199,7 @@ Computes the minimum absolute value of a single-precision floating-point strided
 ```c
 const float x[] = { 1.0f, -2.0f, 3.0f, -4.0f };
 
-float v = stdlib_strided_dmax_ndarray( 4, x, 1, 0 );
+float v = stdlib_strided_sminabs_ndarray( 4, x, 1, 0 );
 // returns 1.0f
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/strided/smskmax/README.md b/lib/node_modules/@stdlib/stats/strided/smskmax/README.md
index 89253cf323b2..e69664a7aaf4 100644
--- a/lib/node_modules/@stdlib/stats/strided/smskmax/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/smskmax/README.md
@@ -232,7 +232,7 @@ Computes the maximum value of a single-precision floating-point strided array ac
 const float x[] = { 1.0f, -2.0f, 2.0f };
 const uint8_t mask[] = { 0, 1, 0 };
 
-float v = stdlib_strided_smskmax( 3, x, 1, 0, mask, 1, 0 );
+float v = stdlib_strided_smskmax_ndarray( 3, x, 1, 0, mask, 1, 0 );
 // returns 2.0f
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/strided/smskmin/README.md b/lib/node_modules/@stdlib/stats/strided/smskmin/README.md
index 9887b4b6f193..6465dbef5fbd 100644
--- a/lib/node_modules/@stdlib/stats/strided/smskmin/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/smskmin/README.md
@@ -232,7 +232,7 @@ Computes the minimum value of a single-precision floating-point strided array ac
 const float x[] = { 1.0f, -2.0f, 2.0f };
 const uint8_t mask[] = { 0, 1, 0 };
 
-float v = stdlib_strided_smskmin( 3, x, 1, 0, mask, 1, 0 );
+float v = stdlib_strided_smskmin_ndarray( 3, x, 1, 0, mask, 1, 0 );
 // returns 1.0f
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/strided/smskrange/README.md b/lib/node_modules/@stdlib/stats/strided/smskrange/README.md
index 420c05ccbf18..2025a0b8dbd2 100644
--- a/lib/node_modules/@stdlib/stats/strided/smskrange/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/smskrange/README.md
@@ -234,7 +234,7 @@ Computes the [range][range] of a single-precision floating-point strided array a
 const float x[] = { 1.0f, -2.0f, 2.0f };
 const uint8_t mask[] = { 0, 1, 0 };
 
-float v = stdlib_strided_smskrange( 3, x, 1, 0, mask, 1, 0 );
+float v = stdlib_strided_smskrange_ndarray( 3, x, 1, 0, mask, 1, 0 );
 // returns 1.0f
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/strided/snanmeanpn/README.md b/lib/node_modules/@stdlib/stats/strided/snanmeanpn/README.md
index 1a92d6f66d35..dddb23449af9 100644
--- a/lib/node_modules/@stdlib/stats/strided/snanmeanpn/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/snanmeanpn/README.md
@@ -203,9 +203,9 @@ console.log( v );
 Computes the [arithmetic mean][arithmetic-mean] of a single-precision floating-point strided array, ignoring `NaN` values and using a two-pass error correction algorithm.
 
 ```c
-const double x[] = { 1.0f, -2.0f, 0.0f/0.0f, 2.0f };
+const float x[] = { 1.0f, -2.0f, 0.0f/0.0f, 2.0f };
 
-double v = stdlib_strided_snanmeanpn( 4, x, 1 );
+float v = stdlib_strided_snanmeanpn( 4, x, 1 );
 // returns ~0.33333
 ```
 
@@ -224,9 +224,9 @@ float stdlib_strided_snanmeanpn( const CBLAS_INT N, const float *X, const CBLAS_
 Computes the [arithmetic mean][arithmetic-mean] of a single-precision floating-point strided array, ignoring `NaN` values and using a two-pass error correction algorithm and alternative indexing semantics.
 
 ```c
-const double x[] = { 1.0f, -2.0f, 0.0f/0.0f, 2.0f };
+const float x[] = { 1.0f, -2.0f, 0.0f/0.0f, 2.0f };
 
-double v = stdlib_strided_snanmeanpn_ndarray( 4, x, 1, 0 );
+float v = stdlib_strided_snanmeanpn_ndarray( 4, x, 1, 0 );
 // returns ~0.33333
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/strided/svariance/README.md b/lib/node_modules/@stdlib/stats/strided/svariance/README.md
index b8cafd9a1509..a15d615ca6e5 100644
--- a/lib/node_modules/@stdlib/stats/strided/svariance/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/svariance/README.md
@@ -264,7 +264,7 @@ Computes the [variance][variance] of a single-precision floating-point strided a
 ```c
 const float x[] = { 1.0f, -2.0f, 2.0f };
 
-float v = stdlib_strided_svariance_ndarray( 4, 1.0f, x, 1, 0 );
+float v = stdlib_strided_svariance_ndarray( 3, 1.0f, x, 1, 0 );
 // returns ~4.3333f
 ```
 
diff --git a/lib/node_modules/@stdlib/stats/strided/svarianceyc/README.md b/lib/node_modules/@stdlib/stats/strided/svarianceyc/README.md
index d46115a72872..5aaa8340434e 100644
--- a/lib/node_modules/@stdlib/stats/strided/svarianceyc/README.md
+++ b/lib/node_modules/@stdlib/stats/strided/svarianceyc/README.md
@@ -264,7 +264,7 @@ Computes the [variance][variance] of a single-precision floating-point strided a
 ```c
 const float x[] = { 1.0f, -2.0f, 2.0f };
 
-float v = stdlib_strided_srange_ndarray( 3, 1, x, 1, 0 );
+float v = stdlib_strided_svarianceyc_ndarray( 3, 1, x, 1, 0 );
 // returns ~4.33333f
 ```
 
diff --git a/lib/node_modules/@stdlib/strided/base/stride2offset/README.md b/lib/node_modules/@stdlib/strided/base/stride2offset/README.md
index 53e18078d071..c6205c0d92bd 100644
--- a/lib/node_modules/@stdlib/strided/base/stride2offset/README.md
+++ b/lib/node_modules/@stdlib/strided/base/stride2offset/README.md
@@ -182,7 +182,7 @@ int main( void ) {
     int64_t stride = -2;
 
     // Compute the offset:
-    int64_t offset = stdlib_strided_stride2offset( N, stride )
+    int64_t offset = stdlib_strided_stride2offset( N, stride );
 
     // Print the results:
     printf( "stride: %"PRId64" => offset: %"PRId64"\n", stride, offset );