% This file calculates parameters for an LBL style coil.
% Zachary Wolf
% 11/9/00

% Coil parameters
r1 = .035005;
r2 = .024997;
r3 = .024569;
r4 = .019447;
Nm = 140;
Nb = 350;

% Main coil sensitivity factors
n = 1:1:18;
Sm = 1 - (-r3/r1).^n;

% Main + bucking coil sensitivity factors
Sb = 1 - (-r3/r1).^n - (Nb/Nm)*(r2/r1).^n + (Nb/Nm)*(-r4/r1).^n;

% Display the results
fprintf('\n\n')
fprintf('  n        Sm          Sb    \n')
fprintf('-----  ----------  ----------\n')
for nn = n
   fprintf('%5i  %10.7f  %10.7f\n', nn, Sm(nn), Sb(nn))
end

% Coil constant (main coil, r0 = r1)
K = 6 / (Nm * r1.^3 * Sm(3));

% Display the result
fprintf('\nMain coil constant, K = %f 1/m^3\n', K)