changed to numpy.minimum in E_bond_pot function