Use a non diagonal metric tensor when integrating the equations of motion of the matter fields. Use general relativity with the matter fields to solve for the metric tensor. Do the above two steps simultaneously.