ols_hsp function

Hocking's Sp