We consider a quasilinear degenerate diffusion–reaction system that describes biofilm formation. The model exhibits two non-linear diffusion effects: a power law degeneracy as one of the dependent variables vanishes and a super diffusion singularity as it approaches unity. Discretisation of the PDE in space by a standard finite volume scheme leads to a singular system of ordinary differential equations. We show that regularisation of this system allows the application of error controlled adaptive integration techniques to solve the underlying PDE. This overcomes the major limitation of existing methods for this type of problem which work with fixed time-steps. We apply the resulting numerical method to study the effect of signal diffusion in the aqueous phase on quorum sensing induction in a biofilm.