%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  magckt2.m - Analysis of Fig 3.40 magnetic circuit.
%              Calculates flux linkage & current arrays
%              for specified value of air gap length x.
%              Leakage flux is considered in calculations.
%              Used by soldes.m.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [LamI] = magckt1(x);
soldata,
if x==0; x=1e-12; else; end; del=0.01*x;
% Air gap & leakage permeances
Pg=4e-07*pi*(pi*d^2/4/x+0.52*pi*d);
Pl=4e-07*pi^2*(d+db)*lp/(db-d)/2;
% Generate left & right flux-mmf curves
npts=25; phimax=Bmax*pi*d^2/4; phi=linspace(0,phimax,npts);
for i=1:npts
   FL(i)=l1*Hcs(phi(i)/A1)+l2*Hcs(phi(i)/A2)+l3*Hcs(phi(i)/A3);
   phiL(i)=phi(i)+FL(i)*Pl;
   FR(i)=l6*Hcs(phi(i)/A6)+l7*Hcs(phi(i)/A7)+l8*Hcs(phi(i)/A8)+ ...
      l9*Hcs(phi(i)/A9);
   phiR(i)=phi(i)+FR(i)*Pl;
end
% Solve center loop looking up FL & FR at same values of flux
for i=1:npts
   mmfR=interp1(phiR, FR, phi(i));
   mmfL=interp1(phiL, FL, phi(i));
   F(i)=l4*Hcs(phi(i)/A4)+l5*Hcs(phi(i)/A5)+(lp/2-x)* ...
      Hcs(phi(i)/A10)+phi(i)/Pg+mmfR+mmfR;
   LamI(1,i)=N*phi(i); LamI(2,i)=F(i)/N;
   philR(i)=mmfR*Pl; philL(i)=mmfL*Pl;
end
% Activate to plot leakage & plunger fluxes
% Can handle from command window: magckt1(x)
% plot(F/N,phi,F/N,philL,F/N,philR); grid;
% ylabel('Flux, Wb'); xlabel('Coil current, A');