1 | initial version |
I don't think there is a built-in solution that does exactly what you want, but it's easy enough to extract a vector solution:
sage: s = solve([x1==4/3, x2==-1/6, x3==-1/6], (x1,x2,x3))
sage: [ x.rhs() for x in s[0] ]
[4/3, -1/6, -1/6]