Enter the energy  in Joules per particle (kilocalories per mol times 4186 J / Cal divided by Avogadro's Number).

k=1.381*10^(-23);

energy=6.0*10^(-19);

Plot[E^(-energy/(k*T)),{T,1,100000},
	PlotRange->{0,1}];