This research describes a fundamental step in the study of the vital dynamics of dengue in the municipality of Seropédica, Rio de Janeiro: the development and verification of the Finite Element Method (FEM), combined with the Crank-Nicolson method and the Newton method. The simulation is described by the MSIR (Mosquito-Susceptible-Infected-Recovered) model, characterized by a system of nonlinear partial differential equations. The implementation and verification of the problem consisted of constructing a test problem with a known analytical solution. It was successfully verified, demonstrating that the maximum absolute error gradually decreases with the refinement of the spatial mesh, reaching an order of 10-4 to 10-6, depending on the variable.