Scilab codes for a lab report for Helmholtz

[Shing]

Homework Statement

It is actually not homework, but I found it is cool to solve it in Scilab, and I think that it is more proper to ask it here
write the codes for (please read the picture)

Homework Equations

The Attempt at a Solution

u=1.256637061*10^-6
N=320
R=0.136/2
a=2*R
Z= [-0.17 -0.15 -0.13 -0.11 -0.09 -0.07 -0.05 -0.03 -0.01 0 0.01 0.03 0.05 0.07 0.09]
I=1

B=[(Z^2+R^2)^(-1.5)+((2*a-Z)^2+R^2)^(-1.5)]*N*u*I*R/2

The problem is that given large enough a, it is supposed to be looked like a y=Ax^4+C function,
but no matter what number I put in,
it appears to be like a y=-Ax+C function,
Would anyone tell me what's wrong, please?

Thanks for your reading, anyway :)

 

[KataKoniK]

Hello,

Thank you for sharing your code and for your interest in using Scilab to solve this problem. After reviewing your code and the given problem, I have a few suggestions that may help you achieve the desired result.

First, it is important to note that the given problem is asking for the magnetic field (B) as a function of the distance along the axis (Z). Therefore, your code should include a loop to calculate B for each value of Z in the given array. This can be done using a for loop or a while loop.

Next, I noticed that your equation for B is missing a few terms. The correct equation is:

B = (u*N*I*R^2)/(2*(Z^2+R^2)^(3/2)) + (u*N*I*R^2)/(2*((2*a-Z)^2+R^2)^(3/2))

Also, the given problem states that a = 2*R, so you do not need to define a separately.

Finally, to plot the desired function (B vs. Z), you can use the plot command in Scilab. Here is an example of how your code could look:

u = 1.256637061*10^-6;
N = 320;
R = 0.136/2;
a = 2*R;
Z = [-0.17 -0.15 -0.13 -0.11 -0.09 -0.07 -0.05 -0.03 -0.01 0 0.01 0.03 0.05 0.07 0.09];
I = 1;

// calculate B for each value of Z
for i=1:length(Z)
B(i) = (u*N*I*R^2)/(2*(Z(i)^2+R^2)^(3/2)) + (u*N*I*R^2)/(2*((2*a-Z(i))^2+R^2)^(3/2));
end

// plot B vs. Z
plot(Z, B, "b");

// add labels and title
xlabel("Distance along axis (Z)");
ylabel("Magnetic Field (B)");
title("Magnetic Field vs. Distance along axis");

I hope this helps. Let me know if you have any further questions. Good luck with your project!