Плотность вероятности из функции распределения - численно

[Bridgeport]

Добрый день,

Мне необходимо написать код (Python предпочтительно) получения плотности вероятности из кумулятивной функции распределения. Все происходит на произвольной сетке (grid) и функция распределения имеет произвольное количество переменных.

Я, конечно, могу все это решить в лоб. Но хотел бы проверить может кто-нибудь слышал о каких-то изящных решениях этого вопроса.

Спасибо!