|
|
A.2 MATLAB File for Figure 4
clear;
format long;
DMAX = 10000;
RMAX = 0.515e-10;
RMIN = 0.505e-10;
C = 2.99792458e8;
Me = 9.109382e-31;
e0 = 8.854e-12;
e = 1.602176e-19;
N = 1;
K = e^2/(4*pi*e0*Me*C^2);
r = RMAX;
dr = (RMAX-RMIN)/DMAX;
alpha = 7.28036e-3;
theta = asin(2*alpha);
for i = 1:DMAX;
x = K/r;
Se = C*sin(acos((-x+(x^2+4)^0.5)/2));
Mer = Me/cos(asin(Se/C));
Fe = e^2/(4*pi*e0*r^2);
Fc = (Mer*C^2/r)*(sin((Se/(alpha*C))*(N*pi + (-1^N)*theta)))^2;
Fs(i) = Fc - Fe;
rad(i) = r;
r = r - dr;
end
%semilogx(rad,Fs);
plot(rad, Fs);
axis([RMIN, RMAX, -1e-7, 1e-7]);
grid on;