plot the 3d magnetic field generated by 1a of current passing through a 6-turn coil in matlab

main.m
% Constants
I = 1; % Current in Amperes
N = 6; % Number of turns

% Create a grid of points where you want to visualize the magnetic field
[x, y, z] = meshgrid(-1:0.1:1, -1:0.1:1, -1:0.1:1);

% Magnetic field calculation
Bx = zeros(size(x));
By = zeros(size(y));
Bz = zeros(size(z));

for i = 1:N
    r = sqrt(x.^2 + y.^2 + z.^2); % Calculate the distance from each point to the coil
    th = acos(z./r); % Calculate the angle theta
    phi = atan2(y, x); % Calculate the angle phi
    
    dBx = (-mu0 * I * a^2 * sin(th) / (2 * pi * r.^2)) .* cos(phi);
    dBy = (-mu0 * I * a^2 * sin(th) / (2 * pi * r.^2)) .* sin(phi);
    dBz = (mu0 * I * a^2 * sin(th) / (4 * pi * r.^3)) .* cos(th);
    
    Bx = Bx + dBx;
    By = By + dBy;
    Bz = Bz + dBz;
end

% Plot the magnetic field using quiver3
figure;
quiver3(x, y, z, Bx, By, Bz);
xlabel('x');
ylabel('y');
zlabel('z');
title('3D Magnetic Field Generated by 1A Current in a 6-turn Coil');
942 chars
34 lines