Unfortunately, I only have Express, so I can't help with the presentation of the symbolic result.
Perhaps, try putting the variables in a vector?
[x_1]
solve, |y_1|
[z_1]
You can assign the result to a variable and transpose it.
However, is there any particular reason you want to use the symbolic processor to solve the equation, when a numeric solution seems to give reasonable results?
The symbolic solve will always provide the solutions row-wise and there is no option to set which will change this behaviour. You will get further rows for non-linear systems and each row represents a solution set.
I also think that in case you are just looking for the numeric values its better to use a numeric method and not the symbolics. Stuart has already shown some ways. Other options are the usage of the "lsolve" function
or using a solve block with "find" (you have to provide arbitrary guess values) which displays its results column-wise as you like it.
Using the symbolics makes sense if you are looking for symbolic solutions (in case your system contains some variables)
and, yes, in this case displaying the results column-wise sure makes sense 😉
Another reason for using the symbolics could be that you want 'exact' results, displayed in fractions (and maybe pi, e, etc) and not a decimal approximation. Unfortunately the symbolics automatically switches into float-mode as soon as just one decimal point is present. So you must also provide your information with fractions, i.e. 3/2, even though 1.5 is also an exact value.
You may also assign the result of the symbolic "solve" a row vector (1 x 3 matrix) with three variables which you then can evaluate symbolically or numerically and use for further calculations.
EDIT: If you really want to use the symbolic engine and insist on a column vector right from the start, there is a way to do it: You have to use a solve block with "find" which is evaluated symbolically:
