def rsolve(a,b): ''' wraps solve, applies to right hand solution solves system x A = B ''' return scipy.linalg.solve(b.T,a.T).T