Mit der Inversen ist es ganz einfach.
Du suchst ja a,b,c,d mit
M * (a,b,c,d)^T = e1 . Multipliziere die Gleichung mit M^(-1)
M^(-1) * M * (a,b,c,d)^T = M^(-1) * e1
==> (a,b,c,d)^T = M^(-1) * e1
Wenn du das entsprechend mit den anderen Basisvektoren
auch machst , dann siehst du, dass für die Matrix N , in deren
Spalten die gesuchten Koordinaten stehen gilt
N = M^(-1) * E (Einheitsmatrix) , also
N = M^(-1).
Die gesuchten Koordinaten sind also in den Spalten
deiner - richtig bestimmten - Matrix M^(-1).