using the same a and b, why solve(a,b) get a 1 by n matrix while solve(a)%*%b get a n by 1 matrix. I thought they should be equal to each other...