some modification

Author: Jim-215-Fisher

doc/specs/stdlib_stats_distribution_normal.md

@@ -64,7 +64,7 @@ program demo_normal_rvs
 ! 0.563655198
 
     print *, norm(1.0, 2.0)
-    !normal random variate miu=1.0, sigma=2.0
+    !normal random variate mu=1.0, sigma=2.0
 
 ! -0.633261681
 
@@ -86,8 +86,8 @@ program demo_normal_rvs
     loc = (-1.0, 2.0)
     scale = (2.0, 1.0)
     print *, norm(loc, scale)
-    !single complex normal random variate with real part of miu=-1, sigma=2;
-	!imagainary part of miu=2.0 and sigma=1.0
+    !single complex normal random variate with real part of mu=-1, sigma=2;
+	!imagainary part of mu=2.0 and sigma=1.0
 
 ! (1.22566295,2.12518454)
 
@@ -149,7 +149,7 @@ program demo_normal_pdf
 ! 0.241970733
 
     print *, norm_pdf(2.0,-1.0, 2.0)
-    !a probability density at 2.0 with miu=-1.0 sigma=2.0
+    !a probability density at 2.0 with mu=-1.0 sigma=2.0
 
 !6.47588000E-02
 
@@ -170,7 +170,7 @@ program demo_normal_pdf
     scale = (1.0, 2.)
     print *, norm_pdf((1.5,1.0), loc, scale)
     ! a complex normal probability density function at (1.5,1.0) with real part
-	! of miu=1.0, sigma=1.0 and imaginary part of miu=-0.5, sigma=2.0
+	! of mu=1.0, sigma=1.0 and imaginary part of mu=-0.5, sigma=2.0
 
 ! 5.30100204E-02
 
@@ -232,7 +232,7 @@ program demo_norm_cdf
 ! 0.841344714
 
     print *, norm_cdf(2.0, -1.0, 2.0)
-    ! a cumulative at 2.0 with miu=-1 sigma=2
+    ! a cumulative at 2.0 with mu=-1 sigma=2
 
 ! 0.933192849
 
@@ -253,7 +253,7 @@ program demo_norm_cdf
     scale = (0.5,1.0)
     print *, norm_cdf((0.5,-0.5),loc,scale)
     !complex normal cumulative distribution at (0.5,-0.5) with real part of
-	!miu=1.0, sigma=0.5 and imaginary part of miu=0.0, sigma=1.0
+	!mu=1.0, sigma=0.5 and imaginary part of mu=0.0, sigma=1.0
 
 !4.89511043E-02
 

src/stdlib_stats_distribution_normal.fypp

@@ -10,9 +10,9 @@ module stdlib_stats_distribution_normal
     private
 
     real(dp), parameter  :: HALF = 0.5_dp, ONE = 1.0_dp, TWO = 2.0_dp
-    integer, save :: kn(0:127)
-    real(dp), save :: wn(0:127), fn(0:127)
-    logical, save  :: zig_norm_initialized = .false.
+    integer :: kn(0:127)
+    real(dp) :: wn(0:127), fn(0:127)
+    logical  :: zig_norm_initialized = .false.
 
     public :: rvs_normal
     public :: pdf_normal
@@ -21,7 +21,7 @@ module stdlib_stats_distribution_normal
 
 
     interface rvs_normal
-    !! Version experimental
+    !! version: experimental
     !!
     !! Normal Distribution Random Variates
     !! ([Specification](../page/specs/stdlib_stats_distribution_normal.html#
@@ -41,7 +41,7 @@ module stdlib_stats_distribution_normal
 
 
     interface pdf_normal
-    !! Version experimental
+    !! version: experimental
     !!
     !! Normal Distribution Probability Density Function
     !! ([Specification](../page/specs/stdlib_stats_distribution_normal.html#
@@ -82,8 +82,6 @@ contains
     ! This is an electronic journal which can be downloaded from:
     ! http://www.jstatsoft.org/v05/i08
     !
-    ! N.B. It is assumed that all integers are 32-bit.
-    !
     ! Latest version - 1 January 2001
     !
         real(dp), parameter :: M1 = 2147483648.0_dp, vn = 0.00991256303526217_dp