This paper introduces a robust algorithm to determine the interfacial tension (gamma) from pendant drop profiles using the Galerkin finite element method (gamma-PD-FEM) to solve the axisymmetric form of the Young-Laplace (YL) equation. In this algorithm, the theoretical profiles are generated by solving the spherical coordinate form of the YL equation. gamma-PD-FEM also solves for the parameter estimates by minimizing the difference between the theoretical and experimental surface functions, f (theta). This technique is compared to the widely used method of converting the YL equation to the three are length-based (ALB) first-order ODEs developed by Bashforth and Adams (BA) in 1883, or as denoted in this paper, the gamma-PD-BA method. The drop apex is the initial condition for the gamma-PD-BA algorithm and the integration is terminated at a specified location along the drop profile. In contrast to techniques based on the BA approach, computation of the theoretical drop profile in gamma-PD-FEM is obtained from a second-order ordinary differential equation and requires boundary conditions at the drop apex and at the contact line of the drop to the nozzle. By incorporating both boundary conditions into the problem formulation, the algorithm can also determine if the drop shape is at static equilibrium. Results to be presented include an outline of the computer algorithm, and comparison of gamma values obtained from the gamma-PD-FEM and the traditional gamma-PD-BA method using simulated and experimental drop profile data sets. (c) 2005 Elsevier Inc. All rights reserved.