A new method to solve a non-steady-state multispecies biofilm model