s=sqrt((sum(samples)-mean(samples))/(num_samples-1))

##</code><code>##

  s=sqrt((num_samples*sum(squares(samples))-square(sum(samples)))/(num_samples*(num_samples-1)))