Octave 3.0.0 (i686-pc-msdosmsvc) The statistics function seems to be incorrectly computing the median: > x = 0:1; > x_stat = statistics( x ); > x_med = x_stat( 3 ); > x_med == median( x ) ans = 0 -Miguel