In generating non-normal data, H. Luo gives the coding to generate non-normal data using Fleishman Power method. 
www.diva-portal.org/smash/get/diva2:407995/FULLTEXT01.pdf. 
But this coding just can be used in specific skewness and kurtosis. 
I would like to simulate the non-normal data for more large values of skewness and kurtosis (for example skewness = 3.5 and kurtosis = -4.0). 

my question is how can i do this or are there any package in R that can generate the data?
Thanks...