BondSquareWell

class BondSquareWell : public feasst::BondTwoBody

U(r) = 0 when |l-l0| < delta/2, otherwise infinity.