%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  magckt1.m - Analysis of Fig 3.34 magnetic circuit.
%              Calculates flux linkage & current arrays
%              for specified value of air gap length x.
%              Used by coenFe.m.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [LamI] = magckt1(x);
% Assumes that input dimensions are MKS units.
w=0.05; d=0.05;      % Core width & depth
lm1=0.20; lm2=0.40;  % Mean length paths
N=100;               % Coil turns
npts=50; B=linspace(0,1.6,npts); LamI=zeros(2,npts);
Pg=4e-7*pi*(w*d/x+0.52*(w+d)+0.308*x);
for i=1:npts
   LamI(1,i)=N*B(i)*w*d;
   LamI(2,i)=(B(i)*w*d/Pg+(lm1-x)*Hcs(B(i))+lm2*Hcs(B(i)))/N;
end