The dataset is about the codes of finite element and finite difference methods, which can be used to study pattern formation instabilities in the generalized reaction-diffusion systems. These codes can simulate a wide spectrum of patterns in geomaterials and geosystems.