Magnetic field of a finite solenoid. Figure 4.9(a) shows a solenoid (cylindrical coi...

Magnetic field of a finite solenoid. Figure 4.9(a) shows a solenoid (cylindrical coil) consisting of N turns of an insulated thin wire wound uniformly and densely in one layer on a cylindrical nonmagnetic support with a circular cross section of radius

The length of the solenoid is l, the current through the wire is I, and the medium is air. The

current over an elemental length dz along the solenoid, shown in Fig.4.9(b), dI = NI dz/l, can be viewed as the current of an equivalent circular current loop [Fig.4.4(b)]. Using Eq.(4.8) and superposition (integration), the total magnetic flux density vector at an arbitrary point P at the solenoid axis [Fig.4.9(b)] comes out to be

where the position of the point is defined by angles θ1 and θ2. Based on this expression, write a function BzFiniteSolenoid() in MATLAB to compute B along the z-axis. (BzFiniteSolenoid.m on IR)