Function H

bem/analytic_integrals.cpp:67–70

  }

  double H (double x, double y, double z)
  {
    return 1/(4*M_PI) * (k1(x,y,z) - y * k0(x,y,z));
  }

  double dHdx (double x, double y, double z)
  {