Fragment formation in heavy ion collisions at low to intermediate energy is described by a combined application of transport theory of the Boltzmann type and of a statistical program for the decay of the fragments at the late stage. The transport equations are solved by simulations using the test particle method as a finite element representation of the phase space distribution. The description of experimental data is reasonable overall, but the fragment velocity distributions point to the presence of other mechanisms and the role of fluctuations.