theta = linspace(0.1,pi);

dcs = sin(theta/2).^(-4); % differential cross section;



plot(180/pi*theta,dcs);

ylim([0 300]);

xlabel('Scattering angle, \phi');

ylabel('Differential cross section ~ sin(\theta/2)^-^4');

title('Expected counting rate curve (Rutherford differential cross section)');

grid on;