We prove existence of a positive solution for a system of non-variational bi-harmonic equations. Furthermore, we give some a priori estimates of solutions and a non-existence result. In addition we compute numerical solutions to illustrate the theoretical results.