R, 7272 61 bytes
function(A,V=t(A/sum(diag(A%*%t(A))))){for(i in 1:1e4)V=2*V-V%*%A%*%V;V}
Porting Sisyphus' answer is not futile at all...and thanks to Sisyphus for -11 bytes.
Observes that \$Tr(AA^T)=\sum\limits_{i,j}a_{ij}^2\$.
R, 94 bytes
function(M)outer(k<-1:dim(M),k,Vectorize(function(j,i)det(M[-i,-j,drop=F])*(-1)^(i+j)))/det(M)
Thanks to Robin Ryder for fixing a bug and making this actually work.
Calculates \$A^{-1}\$ using the adjugate/determinant method.