We develop some simple simulation algorithms for CIR and Wishart processes. We investigate rigorously the square of a matrix valued Ornstein–Uhlenbeck process, the main idea being to split the generator and to reduce the problem to the simulation of the square of a matrix valued...