diff --git a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/cdf/README.md b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/cdf/README.md
index 45a25b2c13f9..f798cb011a17 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/cdf/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/cdf/README.md
@@ -153,6 +153,112 @@ for ( i = 0; i < 10; i++ ) {
+
+
+* * *
+
+
+
+## C APIs
+
+
+
+
+
+
+
+
+
+
+
+### Usage
+
+```c
+#include "stdlib/stats/base/dists/pareto-type1/cdf.h"
+```
+
+#### stdlib_base_dists_pareto_type1_cdf( x, alpha, beta )
+
+Evaluates the cumulative distribution function (CDF) for a Pareto (Type I) distribution with shape parameter `alpha` and scale parameter `beta` at a value `x`.
+
+```c
+double out = stdlib_base_dists_pareto_type1_cdf( 2.0, 1.0, 1.0 );
+// returns 0.5
+```
+
+The function accepts the following arguments:
+
+- **x**: `[in] double` input value.
+- **alpha**: `[in] double` shape parameter.
+- **beta**: `[in] double` scale parameter.
+
+```c
+double stdlib_base_dists_pareto_type1_cdf( const double x, const double alpha, const double beta );
+```
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+### Examples
+
+```c
+#include "stdlib/stats/base/dists/pareto-type1/cdf.h"
+#include
+#include
+
+static double random_uniform( const double min, const double max ) {
+ double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
+ return min + ( v*(max-min) );
+}
+
+int main( void ) {
+ double alpha;
+ double beta;
+ double x;
+ double y;
+ int i;
+
+ for ( i = 0; i < 25; i++ ) {
+ x = random_uniform( 0.0, 8.0 );
+ alpha = random_uniform( 0.0, 5.0 );
+ beta = random_uniform( 0.0, 5.0 );
+ y = stdlib_base_dists_pareto_type1_cdf( x, alpha, beta );
+ printf( "x: %lf, α: %lf, β: %lf, F(x;α,β): %lf\n", x, alpha, beta, y );
+ }
+}
+```
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+