dev_skewnorm function

Skew Normal Deviances