Motivation: To reconstruct metabolic pathways from biochemical and/or genome sequence data, stoichiometric and thermodynamic feasibility of the pathways has to be tested. This is achieved by characterizing the admissible region of flux distributions in steady state. This region is spanned by what can be called a convex basis. The concept of 'elementary flux modes' provides a mathematical tool to define all metabolic routes that are feasible in a given metabolic network. In addition, we define 'enzyme subsets' to be groups of enzymes that operate together in fixed flux proportions in all steady states of the system. Results: Algorithms for computing the convex basis and elementary modes developed earlier are briefly reviewed. A newly developed algorithm for detecting all enzyme subsets in a given network is presented. All of these algorithms have been implemented in a novel computer program named METATOOL, whose features are outlined here. The algorithms are illustrated by an example taken from sugar metabolism. Availability: METATOOL is available from ftp://bmsdarwin.brookes.ac.uk/pub/ software/ibmpc/metatool.