Stabilised the models by forcing the production of k and epsilon to be positive