Random orthogonal matrix simulation targeting Kollo skewness and kurtosis