The parabolic systems for migration of groundwater contaminants are studied through the semi-analytic method (SAM). For systems, a modification of the classic trigonometric basis function (TBF) is proposed that successfully obtains the solution on the available boundary data and improves accuracy. The method is based on a theory that removes instability and keeps the same accuracy. However, especially when using SAM schemes with acceptable stability properties, one is still faced with the considerable task of linearizing the systems, consequently, the linearization of the systems which approxi...