scores_2pnorm function

Calculating scores for the two-piece-normal distribution