Simulation of multivariate diffusion bridges