Given three noncolinear points $P_1=(x_1,y_1)$, $P_2=(x_2,y_2)$ and $P_3=(x_3,y_3)$, an equation of the circle through them is $$\det\begin{pmatrix}x^2+y^2 & x & y & 1 \\ x_1^2+y_1^2 & x_1 & y_1 & 1 \\ x_2^2+y_2^2 & x_2 & y_2 & 1 \\ x_3^2+y_3^2 & x_3 & y_3 & 1 \end{pmatrix} = 0.$$ If this isn’t simple enough, another way is to find the circle’s center $(h,k)$ by intersecting the perpendicular bisectors of two pairs of the points. Once you have that, it’s a simple matter of computing the (square of) the radius and plugging into the generic equation $(x-h)^2+(y-k)^2=r^2$.
This intersection of lines can be computed directly, without solving any equations: Let $\mathbf n_{12}=P_2-P_1$ and $\mathbf n_{23} = P_3-P_2$. Then, compute the cross product $$\left[\mathbf n_{12};-\frac12(P_1+P_2)\cdot\mathbf n_{12}\right] \times \left[\mathbf n_{23};-\frac12(P_2+P_3)\cdot\mathbf n_{23}\right].$$ The two terms in this cross product are the homogeneous vectors that represent the perpendicular bisectors of $P_1P_2$ and $P_2P_3$, respectively, and are closely related to the point-normal form of their equations. The result is the homogeneous coordinates of the circle’s center; divide through by the last component to convert to inhomogeneous Cartesian coordinates. (If the last component is $0$, the two lines are parallel—the three points are colinear.)
Using your updated example, we have $\mathbf n_{12}=(-8,-6)$ and $\mathbf n_{23}=(4,8)$, and the center of the circle is therefore $$[-8,-6,-(0,3)\cdot(-8,-6)]\times[4,8,-(-2,4)\cdot(4,8)] = [-8,-6,18]\times[4,8,-24] = [0,-120,-40],$$ which is the point $(0,3)$. The radius is easily found by computing the distance to the point $(0,8)$, so an equation of the circle is $x^2+(y-3)^2=25$.