29 #include "component.h"
71 if (frequency * h >= 0.03 *
C0) {
73 "freq/C0*h < 0.03 (is %g)\n", frequency /
C0 * h);
78 nr_double_t res = R *
sqrt (1 + frequency * fs);
80 nr_double_t ind =
MU0 * (h *
log ((h + a) / r) + 1.5 * (r - a));
81 return Z =
rect (res, frequency * ind);
91 nr_double_t v = h /
M_PI / (
sqr (r) -
sqr (r - t));
100 nr_double_t g = 1.0 / r;
130 nr_double_t y =
real (1.0 / Z);
132 nr_double_t f =
kelvin (T) /
T0 * 4.0 * y;