Classical Calculation for Mutual Inductance of Two ... - Kurt Nalty

deltar.x = source_coord.x - sensor_coord.x;
deltar.y = source_coord.y - sensor_coord.y;
deltar.z = source_coord.z - sensor_coord.z;
mag_deltar = sqrt(dot(deltar,deltar));
d_mu = dot(source_tangent, sensor_tangent)/mag_deltar;
mu += d_mu*ds;
}
return (mu);
}
double mu(void) // get mutual inductance
{
double theta; // angle along path, cylindrical coordinates
double dtheta; // differential angle along path
double r=0.5; // radius of loop = 1 m diameter
double z_source= 0.5; // position of source loop
double z_sensor=-0.5; // position of sensor loop
double s; // path length along coil
double ds; // differential path length
double mag_deltar; // difference between source and sensor
double mu; // mu_prime
double d_mu; // incremental mu_prime
vector source_coord, source_tangent,deltar;
vector sensor_coord, sensor_tangent;
int nturns = 1; // number of turns
long i,j;
s = 0.0;
mu = 0.0;
for (i=0; i