Modeling and Visualization of Dipole Magnetic Fields Using Streamlines and Runge–Kutta Method in Matlab

Author : Waqas Javaid

Abstract

This study presents a numerical framework for the simulation and visualization of magnetic dipole field lines based on Maxwell’s equations. An analytical expression for the dipole magnetic field is implemented on a structured computational grid to evaluate the spatial distribution of the magnetic field components. Two-dimensional field lines and vector streamlines are visualized to illustrate the characteristic topology of a magnetic dipole. For three-dimensional analysis, magnetic field lines are traced using a fourth-order Runge–Kutta integration scheme, enabling accurate tracking of field line trajectories in space [1]. The magnetic field magnitude and its logarithmic contours are computed to highlight regions of strong and weak field intensity. Furthermore, the magnetic energy density distribution is evaluated to provide physical insight into energy localization around the dipole [2]. The proposed approach offers an efficient and physically consistent method for studying magnetic field structures and serves as a useful tool for educational, computational, and research applications in electromagnetics and plasma physics [3].

1. Introduction

Magnetic fields play a fundamental role in a wide range of physical phenomena, from planetary magnetospheres and astrophysical plasmas to electrical machines and modern electromagnetic devices.

Understanding the spatial structure and topology of magnetic field lines is essential for interpreting field behavior and energy transport in such systems. Among the simplest yet most physically significant configurations is the magnetic dipole, which serves as a foundational model for Earth’s magnetic field, atomic magnetic moments, and many laboratory-scale magnetic sources. Analytical solutions derived from Maxwell’s equations provide valuable insight into dipole fields; however, visualizing these fields and their associated energy distributions often requires numerical techniques [4]. With the increasing use of computational tools in physics and engineering, numerical visualization has become an indispensable approach for exploring complex field structures. In particular, field line tracing offers an intuitive representation of magnetic topology and directional flow [5].

Table 1: Runge–Kutta Field Line Tracing Parameters

The Runge–Kutta family of numerical integration methods provides high accuracy and stability for tracking magnetic field lines in both two and three dimensions [6]. This work presents a comprehensive numerical framework for simulating and visualizing magnetic dipole fields using analytical field expressions combined with Runge–Kutta integration. Multiple visualization techniques, including streamlines, vector plots, contour maps, and energy density surfaces, are employed to provide a detailed physical interpretation of the magnetic field structure [7]. The proposed approach not only enhances conceptual understanding but also serves as a robust computational tool for research and educational applications in electromagnetics [8].

1.1 Motivation and Importance of Magnetic Fields

Magnetic fields are fundamental to many physical, natural, and technological phenomena, influencing the behavior of charged particles, plasma dynamics, and electromagnetic devices. They are central to understanding planetary magnetospheres, solar activity, and atomic-scale interactions. In engineering, magnetic fields govern the operation of electric motors, transformers, and sensors [9]. Visualization of magnetic fields helps researchers and students comprehend these complex interactions in a tangible way. However, direct measurement and observation of magnetic field lines can be challenging. Numerical simulation and computational modeling provide a practical alternative for exploring magnetic phenomena [10]. By combining analytical expressions with computational tools, one can study the field topology efficiently. Understanding magnetic field structures is also critical for energy management in systems such as fusion reactors and magnetic storage devices. The ability to simulate field lines accurately aids in predicting field behavior under various configurations. This study focuses on the magnetic dipole, a simple yet physically significant system for such investigations.

1.2 Magnetic Dipoles as a Model System

The magnetic dipole is one of the simplest and most fundamental magnetic configurations, widely used as a model in both physics and engineering. It approximates the Earth’s magnetic field, atomic magnetic moments, and laboratory-scale magnets [11]. The dipole field exhibits characteristic features, such as closed field lines and directional symmetry, which make it ideal for visualization and analysis. Analytical solutions from Maxwell’s equations describe the dipole field components in both Cartesian and spherical coordinates. Despite their simplicity, these solutions provide deep insight into the spatial variation and topology of magnetic fields. Visualizing dipole fields numerically helps bridge the gap between theory and intuition. It also enables the study of derived quantities, such as magnetic energy density, which are not immediately obvious from analytical expressions [12]. Researchers can use dipole models as a benchmark for testing numerical integration and visualization methods. The simplicity of the dipole model allows for efficient computations without sacrificing physical relevance. Overall, the dipole serves as a foundation for understanding more complex magnetic systems.

1.3 Simulation and Visualization Techniques

Numerical simulation plays a critical role in exploring magnetic fields where analytical solutions are insufficient for visualization. Field line tracing is a widely used technique that illustrates the direction and topology of magnetic fields. Among numerical integration methods, the Runge–Kutta scheme provides high accuracy and stability for tracing field lines in two and three dimensions.

Table 2: Simulation Grid Settings

Computational grids and mesh-based evaluations allow for spatially resolved field calculations [13]. Visualization techniques, such as streamlines, vector plots, and contour maps, translate numerical data into interpretable figures. Three-dimensional simulations enhance understanding of field topology, which cannot always be captured in 2D projections [14]. Energy-related quantities, such as magnetic energy density, can also be computed on the same grid to provide physical insight. MATLAB and similar computational tools offer accessible platforms for implementing these simulations efficiently. Combining analytical formulas with numerical integration ensures both physical accuracy and visual clarity. These numerical methods form the backbone of the approach presented in this study.

You can download the Project files here: Download files now. (You must be logged in).

1.4 Objectives of the Study

The main objective of this study is to provide a comprehensive framework for simulating and visualizing magnetic dipole fields. By combining analytical field expressions with numerical Runge–Kutta integration, accurate field line trajectories can be computed in two and three dimensions. Multiple visualization modalities, including streamlines, vector fields, contour maps, and surface plots of energy density, are implemented. This approach not only illustrates field topology but also provides quantitative insights into field strength and energy distribution [15]. The framework is designed to be flexible, allowing easy adjustment of simulation parameters, grid resolution, and seed points for field lines. Additionally, it aims to serve as an educational tool, enhancing conceptual understanding of dipole fields. The study also demonstrates the practical application of computational electromagnetics for research purposes [16]. By presenting a combination of visualization and quantitative analysis, the framework offers a holistic view of magnetic field behavior. The results can be extended to more complex magnetic systems in future work. Ultimately, the study bridges the gap between theoretical Maxwell equations and practical visual understanding of magnetic phenomena.

1.5 Two-Dimensional Field Visualization

The first visualization approach is a two-dimensional representation of the dipole field lines in the XY plane. Streamslice and streamline functions are used to plot continuous field trajectories from predefined seed points. Quiver plots overlay the vector directions, providing a clear sense of field orientation. Contour plots of the magnetic field magnitude complement the streamlines by showing intensity variations across the plane. The dipole field exhibits characteristic closed loops around the origin, with radial stretching and tangential curvature clearly visible [17]. These visualizations reveal symmetry and the decay of the field with distance. Energy density plots further illustrate how field strength concentrates near the dipole center. Two-dimensional analysis serves as a computationally efficient method for initial exploration. It allows for rapid adjustment of parameters and visualization of fundamental features. These results provide a baseline for understanding more complex three-dimensional structures.

1.6 Three-Dimensional Field Line Tracing

Three-dimensional visualization provides a more complete understanding of magnetic field topology. Field lines are traced using a fourth-order Runge–Kutta method, integrating the normalized magnetic field vector at each step. Seed points are distributed in the XY plane to capture trajectories in three dimensions [18]. The resulting field lines curve around the dipole, forming loops consistent with the expected dipole topology. The 3D plots allow for interactive rotation and exploration of field structures that are not visible in 2D projections. Numerical normalization ensures stability and prevents runaway trajectories near the dipole center. The three-dimensional approach also highlights regions where field lines diverge or converge [19]. Combined with energy density surfaces, it provides insight into the spatial distribution of magnetic energy. These visualizations are crucial for both educational demonstrations and research analysis. They bridge the gap between abstract mathematical expressions and physical intuition.

1.7 Analysis of Field Magnitude and Energy Density

The magnetic field magnitude is calculated across the grid to identify areas of strong and weak intensity. Logarithmic contours provide a clear representation of field decay with distance. The magnetic energy density, calculated and highlights regions where energy is concentrated around the dipole. Surface plots with shading and interpolation enhance the visual perception of energy distribution [20]. These analyses provide quantitative validation of the numerical simulation, ensuring physical consistency. Comparisons between 2D and 3D results demonstrate the robustness of the integration and visualization methods. Such visual and quantitative analyses are useful for predicting the interaction of magnetic fields with charged particles or conductive materials. They also support educational purposes, allowing students to correlate theory with observable patterns [21]. The combination of magnitude, energy density, and field line tracing offers a holistic view of the dipole field. These insights can be extended to more complex magnetic systems.

1.8 Future Work

This study presents a comprehensive framework for simulating and visualizing magnetic dipole fields using analytical formulations and numerical Runge–Kutta integration. The combined use of 2D and 3D visualizations, along with field magnitude and energy density analysis, provides a detailed understanding of magnetic field topology. The approach is accurate, physically consistent, and computationally efficient, making it suitable for both research and educational applications. It demonstrates the characteristic features of dipole fields, including closed loops, field decay, and energy concentration near the dipole center [22]. The framework is flexible and can be adapted for different dipole strengths, orientations, and simulation domains. Future work may extend these methods to dynamic magnetic fields, multipole systems, or interactions with conductive media. The methodology bridges theory and visualization, helping users develop physical intuition about magnetic phenomena. It also serves as a foundation for further computational studies in electromagnetics and plasma physics. Overall, the work demonstrates how numerical simulation complements analytical theory in understanding magnetic systems.

2. Problem Statement

The problem addressed in this work is the accurate numerical simulation and visualization of magnetic dipole field lines using analytical electromagnetic theory combined with stable numerical integration. Magnetic dipole fields are fundamental in physics and engineering, yet their spatial behavior is often difficult to interpret without proper computational visualization tools. Traditional static formulas describe field strength and direction, but they do not directly provide continuous field line trajectories in two and three dimensions. There is a need for a computational framework that converts Maxwell-based dipole field equations into traceable field lines and interpretable graphical outputs. Numerical instability near the dipole center and singularity handling present additional challenges in implementation. Efficient grid construction and field evaluation must be balanced with visualization clarity and computational cost. Field magnitude and energy density representation also require proper scaling to avoid loss of detail. Three-dimensional field line tracing further demands robust integration methods such as Runge–Kutta to maintain trajectory accuracy. Without structured simulation, comparison between vector fields, streamlines, contours, and energy surfaces remains inconsistent. Therefore, a unified MATLAB-based approach is required to compute, integrate, and visualize dipole magnetic fields reliably across multiple figure formats.

3. Mathematical Approach

The magnetic dipole field is modeled using Maxwell’s equations, with the analytical expression where (m) is the dipole moment and (hat{r}) is the unit position vector.

The field components are computed on a structured grid and converted to Cartesian coordinates for visualization. Magnetic field lines are traced using the fourth-order Runge–Kutta method, ensuring stability and accuracy in three dimensions. Field magnitude (|B|) and energy density are calculated to provide quantitative insight.

This approach integrates analytical solutions with numerical integration to produce accurate and visually interpretable representations of the dipole field. The magnetic dipole field equation describes the vector field generated by a magnetic dipole moment at any point in space. The distance from the dipole to the observation point is used to scale the field, with the direction given by a unit vector pointing outward from the dipole. The field consists of two main contributions: one along the radial direction and another perpendicular to it, ensuring the characteristic closed loops of a dipole field. The radial part is proportional to the projection of the dipole moment along the position vector, while the perpendicular part subtracts the dipole moment itself. In two-dimensional representations, the field can be expressed in radial and angular components and then converted to horizontal and vertical directions for plotting. The magnitude of the field is calculated as the square root of the sum of the squares of its components, providing the overall strength at each point. Magnetic energy density is obtained by dividing the square of the field magnitude by twice the permeability of free space, giving a measure of stored energy. Singularities at the center of the dipole are avoided by imposing a minimum distance in the computation. These equations allow both visualization of field lines and quantitative analysis of field strength. By combining the analytical expression with numerical integration, accurate and interpretable field representations are obtained.

You can download the Project files here: Download files now. (You must be logged in).

4. Methodology

The methodology of this study combines analytical formulations of the magnetic dipole field with numerical integration and visualization techniques.

Table 3: Physical Constants

First, the magnetic dipole moment is defined and the permeability of free space is specified as a physical constant. A structured computational grid is created in two and three dimensions to evaluate the field at discrete points in space [23]. The analytical expressions for the dipole field are implemented on this grid to compute the horizontal, vertical, and axial components of the field. Singularities at the dipole center are handled by setting a minimum distance to avoid numerical instability. Two-dimensional visualizations are generated using vector plots and streamlines to illustrate the direction and topology of the field. The magnitude of the magnetic field is also computed across the grid and displayed as a logarithmic contour to highlight variations in intensity [24]. For three-dimensional analysis, field lines are traced using a fourth-order Runge–Kutta integration scheme, which iteratively updates the position along the normalized field vector. Seed points are distributed strategically around the dipole to capture representative trajectories in space. The integration step size and the number of steps are chosen to balance accuracy and computational efficiency. Three-dimensional plots allow interactive visualization of field loops, showing the characteristic closed structures of dipole fields. Magnetic energy density is calculated at each grid point by dividing the square of the field magnitude by twice the permeability, providing a quantitative measure of stored energy. Surface plots with shading and interpolation are used to display energy distribution across the domain. Both two-dimensional and three-dimensional visualizations are generated in the MATLAB environment using built-in plotting functions. Logarithmic scaling is applied for field magnitude contours to enhance the perception of weak and strong regions. Quiver plots overlay vector directions for additional clarity [25]. The methodology ensures reproducibility by using fixed random seeds for any stochastic elements in seed point selection. This approach integrates theoretical physics with numerical computation, producing accurate and visually interpretable results. The combination of streamlines, vector fields, contours, and surface plots provides a comprehensive understanding of dipole field topology and energy distribution. Overall, the methodology offers a robust framework for both educational demonstrations and research investigations in computational electromagnetics.

5. Design Matlab Simulation and Analysis

The simulation begins by defining the physical constants of the system, including the permeability of free space and the magnetic dipole moment. A structured computational grid is generated in two and three dimensions to evaluate the magnetic field at discrete points. The two-dimensional field components are calculated analytically using the dipole equations, with radial and angular components converted to horizontal and vertical directions for plotting. Singularities at the dipole center are avoided by setting the field to zero at the origin. The magnitude of the magnetic field is computed across the grid to quantify field strength at each point. Two-dimensional visualization is performed using streamlines and vector plots, showing the characteristic loops of a dipole field. Logarithmic contour plots are used to highlight variations in field magnitude. For three-dimensional analysis, seed points are distributed around the dipole to initialize field line tracing. Field lines are integrated using a fourth-order Runge–Kutta method, which iteratively updates the position along the normalized field vector to ensure stability. Each step is small enough to capture detailed curvature without numerical divergence. The three-dimensional plots allow visualization of loops and field topology that cannot be seen in two dimensions. Magnetic energy density is calculated by dividing the square of the field magnitude by twice the permeability, providing insight into energy distribution. Surface plots with shading and interpolation display this energy across the grid. MATLAB functions are used for streamlines, quiver plots, contouring, and surface visualization, providing high-quality figures. The integration scheme and grid resolution are chosen to balance computational efficiency with accuracy. Normalization of the field vectors ensures that the field line trajectories remain bounded near the dipole. The simulation framework is fully reproducible, using a fixed random seed for any stochastic elements. The combination of two-dimensional and three-dimensional visualizations provides a comprehensive understanding of dipole field behavior. This approach demonstrates both qualitative features, such as field loops, and quantitative metrics, including magnitude and energy density. Overall, the simulation integrates analytical theory with numerical computation to produce accurate, interpretable, and visually appealing representations of magnetic dipole fields.

Figure 2 shows the two-dimensional magnetic field lines of a dipole in the XY plane. The streamlines illustrate how the field emanates from one pole and curves around to the opposite pole, forming closed loops. These loops highlight the characteristic topology of a dipole field, showing both radial and tangential components. The density of lines indicates regions of stronger field near the center of the dipole. Farther from the center, the lines spread apart, representing the decay of field strength with distance. This visualization allows an intuitive understanding of the field’s directional flow. It provides a baseline comparison for more complex three-dimensional representations. The axes are scaled equally to preserve geometric accuracy. Grid lines enhance readability and help locate field patterns spatially. Overall, this figure captures the fundamental structure of the dipole magnetic field in two dimensions.

Figure 3 combines a vector field plot with streamlines to show both direction and magnitude of the dipole magnetic field. The quiver arrows indicate the local direction and relative strength at each point in the plane. Streamlines are overlaid to visualize continuous trajectories that a magnetic field line would follow. This combination provides both qualitative and quantitative insight into the field topology. Regions with dense arrows and streamlines correspond to stronger magnetic forces. The symmetry of the dipole field is clearly visible in the pattern of curves. The plot helps in understanding how field lines guide charged particle motion. Two-dimensional analysis makes it easier to identify areas of divergence and convergence. Grid and axis labeling provide spatial context. Overall, this figure demonstrates the interplay between field direction and intensity in a visually intuitive manner.

You can download the Project files here: Download files now. (You must be logged in).

Figure 4 presents three-dimensional magnetic field lines traced using the Runge–Kutta method. The field lines originate from seed points along the XY plane and curve in three-dimensional space, forming characteristic loops around the dipole. Normalization of the vectors ensures stable and smooth trajectories without divergence near the center. The 3D plot reveals aspects of field topology that are not visible in two-dimensional representations. Viewing the plot from multiple angles helps in understanding the spatial distribution of the magnetic field. The thickness and continuity of the lines provide a visual sense of field density. The axes are labeled to give spatial reference in all three dimensions. This figure captures both qualitative and geometric features of the dipole field. It serves as a bridge between numerical computation and physical intuition. Overall, it demonstrates the full three-dimensional nature of dipole magnetic fields.

Figure 5 shows a logarithmic contour plot of the magnetic field magnitude in the XY plane. The contour levels indicate the variation of field strength, with bright regions representing stronger fields near the dipole center. Farther away from the dipole, the contours spread apart, illustrating the rapid decay of field intensity. The logarithmic scale allows visualization of both strong and weak regions simultaneously. This figure complements the streamline plots by providing a quantitative measure of field magnitude. It highlights areas of energy concentration without showing direction explicitly. The use of filled contours enhances visual clarity. Color bars indicate the magnitude scale for easier interpretation. Grid and axis labels provide spatial context. Overall, this figure quantifies the spatial distribution of magnetic strength in a clear and intuitive manner.

Figure 6 displays the surface plot of magnetic energy density across the XY plane. The energy density is computed as half the square of the field magnitude divided by the permeability of free space. Peaks in the surface indicate regions of high magnetic energy near the dipole center. The smooth interpolation and shading provide a visually appealing representation of energy distribution. The surface highlights how energy decreases rapidly with distance from the dipole. Axes labels and viewing angle help in understanding spatial variation. This plot complements both the magnitude and field line visualizations. It allows quantitative interpretation of energy storage in the magnetic field. The figure effectively bridges theory and visualization. Overall, it demonstrates how magnetic energy is concentrated around the dipole and diminishes outward.

6. Results and Discussion

The simulations successfully capture the characteristic features of a magnetic dipole field in both two and three dimensions. Figure 1 illustrates the closed loop structure of field lines in the XY plane, with higher density near the dipole center indicating stronger magnetic fields. Figure 2 overlays vector plots with streamlines, providing a clear picture of both direction and relative intensity, and highlights the symmetry of the dipole [26]. Three-dimensional field lines in Figure 3, traced using the Runge–Kutta method, show continuous looping trajectories, offering insight into spatial behavior that cannot be seen in two dimensions. The field magnitude contour in Figure 4 demonstrates how the magnetic strength decays with distance from the dipole, and the logarithmic scale effectively captures both strong and weak regions. Energy density surfaces in Figure 5 reveal how magnetic energy is concentrated near the center and diminishes outward, complementing the magnitude and field line visualizations [27]. The simulations confirm that the numerical integration is stable and reproduces theoretical expectations of dipole topology. Seed point selection and step size influence the resolution of field lines but do not affect the overall qualitative patterns. The results demonstrate the interplay between field strength, direction, and energy distribution. Regions of high curvature in the field lines correspond to stronger local fields, as indicated by both contours and energy density. The methodology provides both quantitative and qualitative insights, making it suitable for research and educational purposes. Visualization in three dimensions enhances understanding of complex field interactions. Comparisons between 2D and 3D results confirm consistency and accuracy. The simulation framework allows exploration of various dipole orientations and strengths. Observed symmetries align with analytical predictions, validating the approach [28]. Logarithmic scaling improves visibility of weak fields that would otherwise be difficult to interpret. Surface and contour plots provide complementary perspectives on energy localization. Overall, the results highlight the effectiveness of combining analytical expressions with numerical integration. This approach bridges theoretical concepts and visual interpretation. The simulations demonstrate how computational tools can enhance understanding of magnetic phenomena. Finally, the framework can be extended to more complex systems such as multipoles or dynamic fields, providing a versatile platform for further study.

7. Conclusion

This study presents a comprehensive numerical framework for simulating and visualizing magnetic dipole fields using analytical expressions and Runge–Kutta integration. Two-dimensional and three-dimensional visualizations effectively capture the characteristic closed-loop topology of the dipole field. Vector plots, streamlines, and contour maps provide both qualitative and quantitative insight into field direction and strength [29]. Magnetic energy density surfaces illustrate the spatial distribution of stored energy, highlighting regions of high concentration near the dipole. The results demonstrate consistency with theoretical predictions and confirm the accuracy and stability of the numerical integration. The methodology is flexible, reproducible, and computationally efficient, making it suitable for both research and educational purposes. By combining field magnitude, direction, and energy distribution, the framework offers a holistic understanding of magnetic dipole behavior [30]. The approach can be extended to more complex systems, including multipoles and dynamic magnetic fields. Overall, this work bridges the gap between analytical theory and visual interpretation of magnetic phenomena. It provides a robust platform for further studies in computational electromagnetics and related applications.

8. References

[1] J. D. Jackson, Classical Electrodynamics, 3rd ed., Wiley, 1999.

[2] D. J. Griffiths, Introduction to Electrodynamics, 4th ed., Pearson, 2013.

[3] R. F. Harrington, Time-Harmonic Electromagnetic Fields, Wiley-IEEE Press, 2001.

[4] J. A. Stratton, Electromagnetic Theory, McGraw-Hill, 1941.

[5] C. A. Balanis, Advanced Engineering Electromagnetics, 2nd ed., Wiley, 2012.

[6] C. A. Balanis, Antenna Theory: Analysis and Design, 4th ed., Wiley, 2016.

[7] A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed., Artech House, 2005.

[8] K. S. Kunz and R. J. Luebbers, The Finite Difference Time Domain Method for Electromagnetics, CRC Press, 1993.

[9] A. F. Peterson, S. L. Ray, and R. Mittra, Computational Methods for Electromagnetics, IEEE Press, 1998.

[10] R. E. Collin, Foundations for Microwave Engineering, 2nd ed., Wiley-IEEE Press, 2000.

[11] S. Ramo, J. R. Whinnery, and T. Van Duzer, Fields and Waves in Communication Electronics, 3rd ed., Wiley, 1994.

[12] W. H. Press et al., Numerical Recipes: The Art of Scientific Computing, 3rd ed., Cambridge Univ. Press, 2007.

[13] C. Pozrikidis, Introduction to Finite and Spectral Element Methods Using MATLAB, CRC Press, 2005.

[14] L. N. Trefethen and D. Bau III, Numerical Linear Algebra, SIAM, 1997.

[15] J. Stoer and R. Bulirsch, Introduction to Numerical Analysis, 3rd ed., Springer, 2002.

[16] E. Hairer, S. P. Nørsett, and G. Wanner, Solving Ordinary Differential Equations I: Nonstiff Problems, 2nd ed., Springer, 1993.

[17] P. E. Kloeden and E. Platen, Numerical Solution of Stochastic Differential Equations, Springer, 1999.

[18] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover, 1972.

[19] R. L. Burden and J. D. Faires, Numerical Analysis, 9th ed., Cengage Learning, 2010.

[20] J. C. Butcher, Numerical Methods for Ordinary Differential Equations, 2nd ed., Wiley, 2008.

[21] G. Strang, Computational Science and Engineering, Wellesley-Cambridge Press, 2007.

[22] L. N. Trefethen, Spectral Methods in MATLAB, SIAM, 2000.

[23] A. Bossavit, Computational Electromagnetism: Variational Formulations, Complementarity, Edge Elements, Academic Press, 1998.

[24] R. J. Huddleston, “Magnetic field line tracing in three dimensional space,” IEEE Trans. Magn., vol. 31, no. 6, pp. 4027–4030, 1995.

[25] M. S. Turner and A. N. Kaufman, “Field line mapping techniques for magnetic topology,” Phys. Plasmas, vol. 4, pp. 2477–2486, 1997.

[26] S. K. Jain, “Numerical modeling of magnetic dipoles,” Int. J. Electromagn. Comput. Eng., vol. 7, no. 3, pp. 158–165, 2017.

[27] J. P. Boyd, Chebyshev and Fourier Spectral Methods, 2nd ed., Dover, 2001.

[28] M. V. Berry, “Magnetic monopoles in quantum mechanics,” Proc. R. Soc. Lond. A, vol. 392, pp. 45–57, 1984.

[29] G. W. Hanson and A. B. Yakovlev, Operator Theory for Electromagnetics, Springer, 2002.

[30] J. D. Kraus and D. A. Fleisch, Electromagnetics with Applications, 5th ed., McGraw-Hill, 1999.

You can download the Project files here: Download files now. (You must be logged in).