A program has been developed in order to fit analytical power series expansions (Dunham, Simon-Parr-Finlan, Ogilvie and their exponential variants) and Pade approximants to discrete ab initio potential energy surfaces of non-linear triatomic molecules. The program employs standard least-squares fitting techniques using the singular decomposition method in order to dampen the higher-order coefficients (if deemed necessary) without significantly degrading the fit. The program makes full use of the symmetry of a triatomic molecule and so addresses the D3h, C2v and C(S) cases.