matlab画出一个载有电流的圆形回路磁偶极子的磁力线

以下是matlab代码实现，画出一个载有电流的圆形回路磁偶极子的磁力线：

% 圆形回路半径为r，电流为I
r = 1;
I = 1;

% 定义磁场计算函数
B = @(x,y) (mu0*I/(2*pi*r^3)) * [(y*sin(theta))/(x^2+y^2-2*r*x*cos(theta)+r^2)^1.5, ...
                                  (x*cos(theta)-r)/(x^2+y^2-2*r*x*cos(theta)+r^2)^1.5, ...
                                  0];

% 生成网格
x = linspace(-2*r,2*r,20);
y = linspace(-2*r,2*r,20);
[X,Y] = meshgrid(x,y);

% 计算磁场
theta = 0;
mu0 = 4*pi*1e-7;
Bx = zeros(size(X));
By = zeros(size(Y));
Bz = zeros(size(X));
for i = 1:numel(X)
    Bxy = B(X(i),Y(i));
    Bx(i) = Bxy(1);
    By(i) = Bxy(2);
    Bz(i) = Bxy(3);
end

% 画图
quiver(X,Y,Bx,By);
axis equal

解释一下代码，首先定义圆形回路半径为r，电流为I。然后定义磁场计算函数B，其中theta为磁场角度，mu0为真空磁导率。接着生成网格，计算出每个点的磁场。最后，用quiver函数画出磁力线。