symsolve function

solve(A, B) where A is symmetric