Compare Spectral Radiance Approximations to Planck's Radiation Law 

Compare the Wien distribution law to Rayleigh–Jeans's law and Planck's radiation law.

In[1]:=

X wiendistributionlaw[t_, \[Nu]_] := UnitConvert[ Quantity[2, "PlanckConstant"/"SpeedOfLight"^2/"Steradians"]*\[Nu]^3* E^(-Quantity[1, "PlanckConstant"/"BoltzmannConstant"]*\[Nu]/ t), "Watts"/"Meters"^2/"Hertz"/"Steradians"]

In[2]:=

X wien = Table[{10^x, QuantityMagnitude@ wiendistributionlaw[Quantity[7000, "Kelvins"], Quantity[10^x, "Hertz"]]}, {x, 13, 15, 0.1}]; rayleighjean = Table[{10^x, QuantityMagnitude@ UnitConvert[ Quantity[2, "BoltzmannConstant"/"SpeedOfLight"^2/"Steradians"]* Quantity[10^x, "Hertz"]^2*Quantity[7000, "Kelvins"], "Watts"/("Hertz"*"Meters"^2*"Steradians")]}, {x, 13, 15, 0.1}]; planck = Table[{10^x, QuantityMagnitude@ PlanckRadiationLaw[Quantity[7000, "Kelvins"], Quantity[10^x, "Hertz"]]}, {x, 13, 15, 0.1}];

In[3]:=

X ListLogLogPlot[{wien, rayleighjean, planck}, Joined -> True, Frame -> True, Axes -> False, PlotLegends -> Placed[{"Wien Distribution Law", "Rayleigh-Jeans's Law", "Planck's Radiation Law"}, Below]]

Out[3]=