This paper provides a description of a package for solving elliptic partial differential equations in two spatial variables. The package is based on the global element method (see Delves and Phillips [1]) and is written in ALGOL 68. The paper presents an example of using GEM2 and sample input and output.