R.DST(x,n,shape,give_log)
x: observation
n: the number of degrees of freedom of the distribution
shape: the shape parameter of the distribution
give_log: if true, log of the result will be returned instead
This function returns the probability density function of the skew-t distribution.