solve(eqn1,...,eqnN) returns solutions to the system of equations with respect to variables found by findsym.

solve(eqn1,eqn2,...,eqnN, var1, var2,..., varN) returns solutions with respect to given variables.

If no closed-form solution is found and the number of equations is equal to number of variables, numeric solve is attempted.

Three different types of output are possible.  For one equation and one output, the resulting solution is returned with multiplicities for a polynomial equation in a vector.  For several variables and an equal number of outputs, the results are sorted and assigned to the outputs.  For several equations and a single output, a MTM[struct] containing the solutions is returned.