Numerical solution of reaction–diffusion systems of lambda-omega type by trigonometrically fitted methods