如何计算一个没有辛的函数的数值混合偏导数？

Question

如何计算函数的混合偏导数：？我知道有一种类似的方法可以做到这一点，但我更喜欢使用标准和已建立的数值技术进行数值计算，以计算任何阶的偏导数（Abramowitz，Stegun Ch 25）scipy.misc.derivative可以对单个变量执行，但不能对混合变量执行

import numpy as np
def f(x, y):
    return (np.sin(x) + np.tanh(y))*x*y