get numerical partial derivatives of matrix
[px py]=egradient(z[,dx[,dy]])
z :
z-matrix
dx :
delta x
dy :
delta y
px :
px=dz/dx
py :
py=dz/dy