sample_sd function

Compute the sample standard deviation