Code corresponding to paper "Robust time-discretisation and linearisation schemes for singular and degenerate evolution systems modelling biofilm growth" implemented in python using FEniCSx. It covers the biofilm model and the porous medium equation.