The given equations could be written in matrix form as:

#((x,y,z,"|",c),(5,-1,-7,"|",-13),(-3,0,8,"|",22),(0,9,1,"|",77))#

#color(white)("XXX")#actually the top row shouldn't be there, but I though it might make the translation more clear.

Cramer's Rule says that if you take the derivatives:

#color(white)("XXX")D = # the matrix composed of the variable coefficients (the left side of the above)

#color(white)("XXX")D_x = # the matrix composed of the variable coefficients with the #x# column replaced with the #c# column.

#color(white)("XXX")D_y = # the matrix composed of the variable coefficients with the #y# column replaced with the #c# column.

#color(white)("XXX")D_z = # the matrix composed of the variable coefficients with the #z# column replaced with the #c# column.

then

#color(white)("XXX")x=|D_x|/|D|color(white)("XXX")y=|D_y|/|D|color(white)("XXX")z=|D_z|/|D|#

Here's what the solution looks like done on a spreadsheet: