I have a pdf that doesn't yield trivial derivatives, so I cannot differentiate it and find the root to determine where its max exactly occurs. However, I have a general formula to express all its moments (mean, var, skew, kurtosis, etc).
The mean doesn't coincide with the mode because the function is not symmetrical. Is there a way to use the mean and the skewness (and/or higher moments) to characterize the mode?