Scribd
Upload
5 views

Chapter3_Prob31

The document presents a mathematical problem involving the force between a charged particle and a charged disk, providing the equation to calculate the force. It outlines three methods for finding the distance z where the force equals 0.3 N, using user-defined functions and MATLAB's built-in function. The solutions obtained from all methods yield a distance of approximately 0.1692 m.

Uploaded by

josa.palaganas.up
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
5 views

Chapter3_Prob31

The document presents a mathematical problem involving the force between a charged particle and a charged disk, providing the equation to calculate the force. It outlines three methods for finding the distance z where the force equals 0.3 N, using user-defined functions and MATLAB's built-in function. The solutions obtained from all methods yield a distance of approximately 0.1692 m.

Uploaded by

josa.palaganas.up
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

1

3.31 The force F acting between a particle with a charge q and a round disk with a radius R and a charge
Q is given by the equation:
F = ----------  1 – --------------------- 
Qqz z
2ε 0  z + R2 
2

where ε 0 = 0.885 × 10 –12 C2/(Nm2) is the permittivity constant and z is the distance to the particle. Deter-
mine the distance z if F = 0.3 N, Q = 9.4 × 10 –6 C, and q = 2.4 × 10 –5 C, and R = 0.1 m.
(a) Use the user-defined function BisectionRoot that was developed in Problem 3.16 with a starting
interval of [ 0.1, 0.2 ] .
(b) Use the user-defined function SteffensenRoot from Problem 3.24.
(c) Use MATLAB’s built-in function fzero.
Solution
The solution is the zero of the function: 0.15

0.1

f ( z ) = ----------  1 – ---------------------  – 0.3


Qqz z 0.05

2ε 0  z + R2 
2 0

f(z)
-0.05

The following script file solves the problem. When the script is -0.1

executed the figure on the right is displayed. The figure shows -0.15

that the solution is between z = 0.1 m and z = 0.2 m. -0.2


0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
z (m)

% HW 3.31 3ed
clear, clc
e0=8.85E-12; R=0.1;
Adisc=pi*R^2;
Q=9.4E-6; q=2.4E-5;R=0.1; F=0.3;
Fun=@ (x) Q*q*x/(2*e0)*(1-x/sqrt(x^2+R^2))-F;
fplot(Fun,[0.1 0.5])
xlabel('z (m)'); ylabel('f(z)')
disp('Part (a)')
z_a = BisectionRoot(Fun,0.1,0.2)
disp('Part (b)')

Excerpts from this work may be reproduced by instructors for distribution on a not-for-profit basis
for testing or instructional purposes only to students enrolled in courses for which the textbook
has been adopted. Any other reproduction or translation of this work beyond that permitted by
Sections 107 or 108 of the 1976 United States Copyright Act without the permission of the
copyright owner is unlawful.
2

z_b = SteffensenRoot(Fun,0.2)
disp('Part (c)')
z_c=fzero(Fun,0.2)

The following is displayed in the Command Window when the script is executed:

Part (a)
z_a =
0.1692
Part (b)
z_b =
0.1692
Part (c)
z_c =
0.1692

Excerpts from this work may be reproduced by instructors for distribution on a not-for-profit basis
for testing or instructional purposes only to students enrolled in courses for which the textbook
has been adopted. Any other reproduction or translation of this work beyond that permitted by
Sections 107 or 108 of the 1976 United States Copyright Act without the permission of the
copyright owner is unlawful.

You might also like