How to generate tri-variate correlated normal random numbers in C++?

Does anyone know how to generate a trivariate correlated normal random numbers (X1, X2, X3) in c++? with correlations equal to rho_1, rho_2, and rho_3?