Solution of phase space diffusion equations using interacting trajectory ensembles