new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/corrcoef-1.m
@@ -0,0 +1,4 @@
+x = rand (10);
+cc1 = corrcoef (x);
+cc2 = corrcoef (x, x);
+size (cc1) == [10, 10] && size (cc2) == [10, 10] && abs (cc1 - cc2) < eps

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/corrcoef-2.m
@@ -0,0 +1,1 @@
+corrcoef ()

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/corrcoef-3.m
@@ -0,0 +1,1 @@
+corrcoef (1, 2, 3)

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/cov-1.m
@@ -0,0 +1,4 @@
+x = rand (10);
+cx1 = cov (x);
+cx2 = cov (x, x);
+size (cx1) == [10, 10] && size (cx2) == [10, 10] && cx1 == cx2

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/cov-2.m
@@ -0,0 +1,1 @@
+cov ()

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/cov-3.m
@@ -0,0 +1,1 @@
+cov (1, 2, 3)

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/kurtosis-1.m
@@ -0,0 +1,3 @@
+x = [-1; 0; 0; 0; 1];
+y = [x, 2*x];
+all (abs (kurtosis (y) - [-1.4, -1.4]) < sqrt (eps))

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/kurtosis-2.m
@@ -0,0 +1,1 @@
+kurtosis ()

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/kurtosis-3.m
@@ -0,0 +1,1 @@
+kurtosis (1, 2)

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/mahalanobis-2.m
@@ -0,0 +1,1 @@
+mahalanobis ()

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/mahalanobis-3.m
@@ -0,0 +1,1 @@
+mahalanobis (1, 2, 3)

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/mean-1.m
@@ -0,0 +1,4 @@
+x = -10:10;
+y = x';
+z = [y, y+10];
+mean (x) == 0 && mean (y) == 0 && mean (z) == [0, 10]

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/mean-2.m
@@ -0,0 +1,1 @@
+mean ()

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/mean-3.m
@@ -0,0 +1,1 @@
+mean (1, 2)

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/median-1.m
@@ -0,0 +1,9 @@
+x = [1, 2, 3, 4, 5, 6];
+x2 = x';
+y = [1, 2, 3, 4, 5, 6, 7];
+y2 = y';
+
+(median (x) == median (x2) && median (x) == 3.5
+ && median (y) == median (y2) && median (y) == 4
+ && median ([x2, 2*x2]) == [3.5, 7]
+ && median ([y2, 3*y2]) == [4, 12])

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/median-2.m
@@ -0,0 +1,1 @@
+median ()

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/median-3.m
@@ -0,0 +1,1 @@
+median (1, 2)

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/skewness-2.m
@@ -0,0 +1,1 @@
+skewness ()

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/skewness-3.m
@@ -0,0 +1,1 @@
+skewness (1, 2)

--- a/test/octave.test/stats/stats.exp
+++ b/test/octave.test/stats/stats.exp
@@ -1,10 +1,100 @@
 ## The following functions need to have tests written for them:
 
-# mean
-# median
-# std
-# cov
-# corrcoef
-# kurtosis
-# mahalanobis
-# skewness
+# mahalanobis (test for success)
+# skewness (test for success)
+
+set test mean-1
+set prog_output "ans = 1"
+do_test mean-1.m
+
+set test mean-2
+set prog_output "usage:.*"
+do_test mean-2.m
+
+set test mean-3
+set prog_output "usage:.*"
+do_test mean-3.m
+
+set test median-1
+set prog_output "ans = 1"
+do_test median-1.m
+
+set test median-2
+set prog_output "usage:.*"
+do_test median-2.m
+
+set test median-3
+set prog_output "usage:.*"
+do_test median-3.m
+
+set test std-1
+set prog_output "ans = 1"
+do_test std-1.m
+
+set test std-2
+set prog_output "usage:.*"
+do_test std-2.m
+
+set test std-3
+set prog_output "usage:.*"
+do_test std-3.m
+
+set test cov-1
+set prog_output "ans = 1"
+do_test cov-1.m
+
+set test cov-2
+set prog_output "usage:.*"
+do_test cov-2.m
+
+set test cov-3
+set prog_output "usage:.*"
+do_test cov-3.m
+
+set test corrcoef-1
+set prog_output "ans = 1"
+do_test corrcoef-1.m
+
+set test corrcoef-2
+set prog_output "usage:.*"
+do_test corrcoef-2.m
+
+set test corrcoef-3
+set prog_output "usage:.*"
+do_test corrcoef-3.m
+
+set test kurtosis-1
+set prog_output "ans = 1"
+do_test kurtosis-1.m
+
+set test kurtosis-2
+set prog_output "usage:.*"
+do_test kurtosis-2.m
+
+set test kurtosis-3
+set prog_output "usage:.*"
+do_test kurtosis-3.m
+
+# set test mahalanobis-1
+# set prog_output "ans = 1"
+# do_test mahalanobis-1.m
+
+set test mahalanobis-2
+set prog_output "usage:.*"
+do_test mahalanobis-2.m
+
+set test mahalanobis-3
+set prog_output "usage:.*"
+do_test mahalanobis-3.m
+
+# set test skewness-1
+# set prog_output "ans = 1"
+# do_test skewness-1.m
+
+set test skewness-2
+set prog_output "usage:.*"
+do_test skewness-2.m
+
+set test skewness-3
+set prog_output "usage:.*"
+do_test skewness-3.m

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/std-1.m
@@ -0,0 +1,3 @@
+x = ones (10, 2);
+y = [1, 3];
+std (x) == [0, 0] && abs (std (y) - sqrt (2)) < sqrt (eps)

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/std-2.m
@@ -0,0 +1,1 @@
+std ()

new file mode 100644
--- /dev/null
+++ b/test/octave.test/stats/std-3.m
@@ -0,0 +1,1 @@
+std (1, 2)