ax = math.atan2(-M[j, k],  M[j, j])