This paper describes an add-on procedure for the calculation of critical points in multicomponent mixtures according to the method by Heidemann and Khalil [1]. The benefit of this approach is that all equations of state (EOS) and mixing rules of an existing program for calculating phase equilibria, which can be a commercial product or self-made, are available for the calculation of critical points. The method by Heidemann and Khalil requires the first and second partial derivatives of the fugacities with respect to the male numbers (delta ln f(k)/delta n(i))(T,V,ni not equal j)(delta 2 ln f(k)/delta n(l)delta(nj))(T,V,ni not equal i,j). In our work the first partial derivatives are obtained numerically with a four-point differencing scheme and the second partial derivatives are obtained from the first partial derivatives with a numerical directional derivative as suggested by Michelsen [2]. This new approach allows the combination of the method by Heidemann and Khalil with any program for calculating phase equilibria and eliminates the need for the cumbersome determination of the analytical first partial derivatives of the fugacity with respect to the mole numbers. The implementation into the commercially available process simulator ASPEN PLUS Version 9.2 [3] is described and some calculation results are shown.