For gases I always found PR very good. L-K in fugacity is ok, just makes life a bit more complicated by adding mixing rules etc. If you look at our book "Thermophysical Properties of Fluids. An Introduction to their Prediction", at Imperial College Press, we have there fortran routines that will calculate both for hydrocarbon mixtures.
right now developed a computerized calculation module for comparing results of different models of equations of state (VdW, RK, RKS, LK, PR, etc.) applied to gas mixtures to estimate thermodynamic properties ... But I could not locate the application of the Lee-Kesler equation to determine the fugacity coefficient of each component in the mixture ...
is the solution in the book "Thermophysical Properties of Fluids. An Introduction to Their Prediction"?
Yep, in FORTRAN. If you give me an Email I can send you the examples computer code. but you need the book for the example headings. Yes I do compare RKS, PR, BWR and LK.
Using the L-K may be important in specific cases, such as associating species in the vapor phase. As mentioned by many others if this is not the case using PR, or or RKS (Redlich-Kwong-Soave) may be if not better just as good.