Magnetic flux density vector due to a moving charge. Expand the code developed in...

Magnetic flux density vector due to a moving charge.

Expand the code developed in the previous MATLAB exercise to also calculate and plot the magnetic force on the charge Q₁ (due to Q₂), as well as the magnetic flux density vector B₁ due to the moving charge Q₁ at the location of the charge Q₂, based on Eq.(4.4), and vice versa (vector B₂ due to Q₂). The plot should contain the charges and vectors v₁, v₂, B₁, B₂, F_m12, and F_m21. (ME4 3.m on IR)

Reference: Equation(4.4)