Shape-dependent strain energy of an inhomogeneous coherent inclusion with general tetragonal misfit strains