def simple_bet(pressure, n_monolayer, c_const):
    """A simple BET equation returning loading at a pressure."""
    return (n_monolayer * c_const * pressure / (1 - pressure) / (1 - pressure + c_const * pressure))