AngleSquareWell

class AngleSquareWell : public feasst::BondThreeBody

U(theta) = 0 when |theta-theta0| < delta/2, otherwise infinity. where theta is in degrees