Our main goal is to study the behavior of the error in finite element approximations of partial differential equations. The error typically has two components -- local error and global (pollution) error. We also will discuss the very interesting phenomena of superconvergence, and in particular, how to determine superconvergence points, and what advantages can be derived from them.