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