diff --git a/src/stdlib_experimental_stats.md b/src/stdlib_experimental_stats.md
index b09f9a0ec..00d06b801 100644
--- a/src/stdlib_experimental_stats.md
+++ b/src/stdlib_experimental_stats.md
@@ -33,7 +33,7 @@ If *dim* is absent, a scalar with the mean of all elements in *array* is returne
 
 ```fortran
 program test
-    use stdlib_experimental_stat, only: mean
+    use stdlib_experimental_stats, only: mean
     implicit none
     real :: x(1:6) = [ 1., 2., 3., 4., 5., 6. ]
     print *, mean(x)                            !returns 21.