code: def adet(z, x):
    
    adjugate = numpy.linalg.det(x) * numpy.linalg.pinv(x)
    d[x] = d[z] * numpy.transpose(adjugate)
docs: d|A|/dA = adj(A).T
  
    See  Jacobi's formula: https://en.wikipedia.org/wiki/Jacobi%27s_formula