Renormalisation group recursion relations are obtained to two-loop order for an n-component Landau-Ginzburg-Wilson model containing both hypercubic anisotropy and quenched random impurities. The fixed points are enumerated and their stability properties analysed, using the epsilon expansion about four spatial dimensions to second order. The existence of a new fixed point corresponding to an anisotropic impure system is reported. It is found that both isotropic and anisotropic pure systems are unstable to the addition of random impurities for n<n_R identical to 4-4 epsilon +O( epsilon ²), while both pure and random isotropic systems are unstable to anisotropy for n>n_C identical to 4-2 epsilon +O( epsilon ²). Singularities in the epsilon expansions for the exponents of the random fixed points are shown to be associated with the onset of focal behaviour (the acquisition of complex eigenvalues). These, together with the slow convergence of the expansions leads one to doubt the reliability of extrapolations to three dimensions.