% File: rutherfordsx.m
% Author: Scott Sewell
% Date: 01 Aug 2005
% 
% Models Rutherford Scattering Cross Section 
% Melissinos Equation 2.17 (1st edition page 235)


e1 = 2;  % z of alpha particle
e2 = 21; % z of titanium
E=5.4e6;   % energy of alpha particle from Am-241
theta=[5*pi/180:5*pi/180:pi/2]
prefactor = (e1*e2/2)^2 * (1/(4*E^2));
angulardependence = 1./(sin(theta/2).^4);
xs= prefactor* angulardependence*10^4;   %10^4 factor is to convert cross section from m^2 to cm^2
subplot(1,2,1);
semilogy(theta*(180/pi),angulardependence,'+');
subplot(1,2,2);
semilogy(theta*(180/pi),xs,'+');