(*******************************************************
	Aquest programa dibuixa l´espectre rotacional teòric
per a una molècula diatómica.
	Els nivells d´energia rotacional es calculen mitjançant
la següent expressió, 
  on s´ha rebutjat el terme de distorsió
centrífuga:
 					E(J) = h c B J (J+1),
	
	Cal recordar que les transicions rotacionals per a molècules 
diatòmiques estan restringides per la següent regla de 
      sel.lecció:

						J1 - J0 = +1
						
	Tenint en compte que la intensitat d´una banda espectral
corresponent a una determinada transició des d´un nivell 
    d´energia
E dependeix de la població del nivell, 
  podem fer ús de la distribució
de Boltzman per tal d´establir les intensitats.

Data: 26/10/00 

********************************************************)

Off[General::spell, General::spell1];
ClearAll["Global`*"];
ClearAll[n,J,B,k,T];
(*Paràmetres*)

k=1.38066*10^(-23);

T=298; 
h=6.62618 10^(-34);
c=2.997925 10^8;

B=0.8*10^(-22)/(h c);
(*Determinació de la població relativa dels nivells espectrals*)

n[J_]:= (2*J+1)*Exp[-h c B*J*(J+1)/(k*T)];
f0[J_]:= 2 B (J+1);
loren[J_,f_]:=aa n[J]/(aa+(f-f0[J])^2);



Plot[n[J],{J,0,15},AxesLabel->{"J","NJ/N0"}];
(*Dibuix de les 15 primeres transicions  l´espectre amb lorentzianes *)
aa=5000;
f[x_]=0;
For[i=1,i<=15,i++,
aux=loren[i,x];
f[x_]=f[x]+aux];
Plot[f[x],{x,0,15000},PlotRange->All,AxesLabel->{"ν","I"}];