How can I derive omega (scale) in R (for use in sn::psn)?

Question

I\'m not good at reading math notation. I have an understanding of stats, but unfortunately it stops short of skewed distributions (really only covering normal distributions an