rawmoment: small change in spec

Author: jvdp1

src/stdlib_experimental_stats.md

@@ -105,21 +105,15 @@ program demo_moment
     use stdlib_experimental_stats, only: mean, moment
     implicit none
     real :: x(1:6) = [ 1., 2., 3., 4., 5., 6. ]
-    print *, moment(x, 2, center = mean(x))
-    !returns 2.9167
-
-    print *, moment( reshape(x, [ 2, 3 ] ), 2,&
-                     center = mean(reshape(x, [ 2, 3 ])))
-    !returns 2.9167
-
-    print *, moment( reshape(x, [ 2, 3 ] ), 2, 1,&
-                     center = mean(reshape(x, [ 2, 3 ] ), 1))
-    !returns [0.25, 0.25, 0.25]
-
-    print *, moment( reshape(x, [ 2, 3 ] ), 2, 1,&
-                     center = mean(reshape(x, [2, 3]), 1, reshape(x, [ 2, 3 ] ) > 3.),&
-                     mask = (reshape(x, [ 2, 3 ] ) > 3.))
-    !returns [NaN, 0., 0.25]
+    real :: y(1:2, 1:3) = reshape([ 1., 2., 3., 4., 5., 6. ], [ 2, 3])
+    print *, moment(x, 2, center = mean(x))                 !returns 2.9167
+    print *, moment( y, 2,&
+                     center = mean(y))                      !returns 2.9167
+    print *, moment( y, 2, 1,&
+                     center = mean(y, 1))                   !returns [0.25, 0.25, 0.25]
+    print *, moment( y, 2, 1,&
+                     center = mean(y, 1, y > 3.),&
+                     mask = (y > 3.))                       !returns [NaN, 0., 0.25]
 end program demo_moment
 ```