In[]:=
T=5800;λmin=50*10^(-9);λmax=3000*10^(-9);L=(2*h*c^2)/(λ^5)*1/(Exp[(h*c)/(k*λ*T)]-1);h=6.626*10^(-34);c=2.998*10^8;k=1.381*10^(-23);Plot[L,{λ,λmin,λmax},AxesLabel{"λ in m","L in J/(m^3*s)"},PlotRangeAll]gesamt=Integrate[L,{λ,0,Infinity}];Framed[Labeled[gesamt,"Die gesamte Strahldichte pro Raumwinkeleinheit ist:",Top]]
Out[]=
Out[]=