SC-Rutherford.txt 895 Bytes
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
function CrsSecRf(Z1, Z2: byte; Ein, M1, M2, Theta: double): double;
{****************************************************************************
    Calcula a seção de choque Rutherford
 ----------------------------------------------------------------------------
  Z1 -- número atômico do projétil
  Z2 -- número atômico do alvo
  Ein -- energia do projétil imediatamente antes da colisão em eV
  M1 -- massa do projétil em a.m.u.
  M2 -- massa do alvo em a.m.u.
 ****************************************************************************}
begin
  if M1 < M2 then
    result := sqr((Z1*Z2*e2*1e8)/Ein)/4 * power(sin(Theta),-4)
            * sqr(sqrt(1 - sqr((M1/M2)*sin(Theta))) + cos(Theta))
            / sqrt(1 - sqr((M1/M2)*sin(Theta))) {angstrom^2}
  else
    result := 0;
    
end;

// ----------------------------------------------------------------------------