Half measure for factor function

Fixed free_symbols attributes of BosonOp and MultiBosonOp so as to work for factor function. (Just a temporary expedient)