... приведите в качестве контрпримера сколь-либо известное распределение, плотность которого определяется моментами порядка выше 4-го.
Примером плотности, для задания параметров которой требуется больше четырех моментов, является смесь двух одномерных нормальных распределений

.
Для задания параметров требуется пять моментов. Задача определение параметров через первые пять моментов приводит к задаче поиска подходящего решения алгебраического уравнения девятой степени.
Однако из самого дифференциального уравнения, определяющего семейство плотностей Пирсона, вытекает сильное качественное ограничение на плотности, которые могут входить в это семейство — это плотности с одной модой (поскольку производная плотности, при положительном значении плотности, может обращаться в ноль только при одном значении аргумента:

). Отсюда, кстати, вновь получаем, что приведенная выше смесь нормальных распределений не может, в общем случае, принадлежать семейству распределений Пирсона.