pub fn wsum(n: usize) -> f64
Sum of linear weights 1..n. Also the size of an upper or lower triangle of a square array (including the diagonal)
