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;