i calculated the bulk modulus by changing a value, and at that time i wondered that could the b and c value could be changed for the same purpose, then I realized a serious problem, if we only changed one parameter (a), it may trapped into local optimum which the bulk value may not be accuracy.