I want to solve the following system of equations:
$$XY-YX=-aY+bZ$$
$$ZX-XZ=cY+aZ$$
$$YZ-ZY=dX$$
Where $X,Y,Z$ are all real $3\times3$ matrices which are all different and not a scalar multiple of eachother. $a,b,c,d$ are given and take the values of either $0,-1,1$. Moreover
I desire $X,Y,Z$ to be as simple as possible, most likely having less then three entries, but this is not a necessity.
The problem arises from trying to find a basis of the standard representation of a three dimensional Lie algebra with the above commutator relations.
Currently the best way for me to solve this is to write a program which iterates through possible matrices with coefficient $0,1,-1$ and while this usually gives a solution it takes a very long time.
See here for an example of sage doing it.