Função de probabilidade do desvio no átomo de Thomson

u=NIntegrate[ Exp[-x^2/(2 4 10^-4)]Sin[x],{x,0,Pi}]

f[q_]:=NIntegrate[ Exp[-x^2/(2 4 10^-4)]Sin[x],{x,0,q}]/u

Plot[f[q],{q,0,0.08},AxesLabel->{"θ(rd)","∫P(θ) Sin[θ]dθ"},BaseStyle->{FontSize->14}]

Plot[f[q Degree],{q,0,5},AxesLabel->{"θ(graus)","∫P(θ) Sin[θ]dθ"},BaseStyle->{FontSize->14}]

f[3.5Degree]

átomo de Rutherford

v=NIntegrate[Sin[x]/Sin[x/2]^4,{x,3/10^3,Pi}]

r[q_]:=NIntegrate[Sin[x]/Sin[x/2]^4,{x,3/10^3,q}]/v

Plot[r[q ],{q,3/10^3,0.1},PlotRange->All,AxesLabel->{"θ(rd)","∫P(θ) Sin[θ]dθ"},BaseStyle->{FontSize->14}]

Plot[{r[q  Degree],f[q Degree]},{q,0.18,5 },PlotRange->All,AxesLabel->{"θ(graus)","∫P(θ) Sin[θ]dθ"},BaseStyle->{FontSize->14}, PlotLegends->{"Rutherford","Thomson"}]

Plot[{r[q  Degree],f[q Degree]},{q,30,180},PlotRange->All,AxesLabel->{"θ(graus)","∫P(θ) Sin[θ]dθ"},BaseStyle->{FontSize->14}, PlotLegends->{"Rutherford","Thomson"}]

Plot[{r[q  Degree],f[q Degree]},{q,3,30},PlotRange->All,AxesLabel->{"θ(graus)","∫P(θ) Sin[θ]dθ"},BaseStyle->{FontSize->14}, PlotLegends->{"Rutherford","Thomson"}]