Как посчитать частные производные функции по элементам матрицы?

dooble, вот пример функции для расчета производной в точке функции 2 перменных по формуле Тейлора:

import numpy as np

h = 1e-8
#----------   Функция расчета значений частных производных по формуле Тейлора   ---------- 
# f - функция 
# p - точка
def calc_gradient(f, p):
    (x, y) = p
    return np.array([
        (f(x + h, y) - f(x - h, y)) * 0.5 / h,
        (f(x, y + h) - f(x, y - h)) * 0.5 / h])