#include <cmath>
#include <stdexcept>
#include "coupling.h"
#include "constants_math.h"
Implementation of coupling