comment about finding inverse with solve() rather than ginv()