Coord dx(Coord x1, Coord x2) { return(x2 - x1); }
Coord s(Coord x1, Coord x2) { return(A * Gamma(x1, x2)); }