The first thing to notice is that you have four equations in five variables -- this means you will have an infinite number of solutions or zero solutions. In the case of an infinite number, you can solve for the set of solutions in terms of one parameter.
The easiest way to do this by hand is to reduce the augmented matrix $\left( A | b\right)$ to reduced row-echelon form. This means proceed with Gaussian elimination until it reaches the form $ \left( \begin{array}{cccccc} 1 && 0 && 0 && 0 && \star && \star\\ 0 && 1 && 0 && 0 && \star && \star \\ 0 && 0 && 1 && 0 && \star && \star\\ 0 && 0 && 0 && 1 && \star && \star \end{array} \right) $
where $\star$ is some unknown value.
For your matrix, this results in $ \left( \begin{array}{cccccc} 1 && 0 && 0 && 0 && 16.5 && -52.5\\ 0 && 1 && 0 && 0 && -5 && 15\\ 0 && 0 && 1 && 0 && 2 && -2\\ 0 && 0 && 0 && 1 && 1 && -1 \end{array} \right) $
Thus, we see that $ \begin{align*} x_1 + 16.5x_5 &= -52.5 \\ x_2 - 5 x_5 &= 15 \\ x_3 + 2 x_5 &= -2 \\ x_4 + x_5 &= -1 \end{align*} $
which I assume is the "general" solution. For a particular solution you can just pick any $x_5$ and solve for the remaining variables.