如何用R算部标准的skew normal distribution的四阶矩呀？skew normal distribution的密度函数形式本身就是不定积分，这时候应该怎么用R处理呢？