Benjamin Morgan have VASP-Utilities code that can calculate B-M equation in bash/ruby. Please check their Github address: https://github.com/bjmorgan/VASP-Utilities
has different functions and among these the calculation of bulk moduli as well as its derivatives is allowed, according to different equation of states. The only thing needed is Energy vs Volume datafile so in principle it can be used with any periodic DFT code. For an example of its application for the calculation of the Bulk modulus of some metal you can see the attached reference.
Roberto Schimmenti
Gibbs2 user.
Article Investigation of Polyol Adsorption on Ru, Pd, and Re Using v...