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)"},PlotRangeAll]
Out[]=