numpy / scipy中非线性函数的数值梯度

I'm trying to implement an numerical gradient calculation in numpy to be used as the callback function for the gradient in cyipopt. My understanding of the numpy gradient function is that it should return the gradient calculated at a point based on a I'm trying to implement an numerical gradient c