A numerical method for solving differential equations in which the original equations are approximated by a set of algebraic equations that are solved on a network of discrete grid points.