(* ROZKŁAD PLANCKA W FUNKCJI DŁ. FALI *)
h=6.62607004*10^(-34);
c=299792458;
kb=1.38064852*10^(-23);
PlanckLambda[λ_,T_]:=(2*Pi*h*c^2 )/( λ^5 (Exp[( h*c)/(λ*kb*T)]-1));

Plot[{PlanckLambda[λ,5700],PlanckLambda[λ,5000],PlanckLambda[λ,4500]},{λ,0,3*10^(-6)},PlotRange->Full,AxesLabel->{"λ  [m]","(Subscript[e, λ]) [W/(m^3)]"},PlotLegends->{"T=5700 K","T=5000 K","T=4100 K"}]
(* dł. fali λmax odpowiadająca maksimum (Subscript[e, λ]) *)
λmax=Quantity[x/.Last[FindMaximum[PlanckLambda[x,5700],{x,1*10^(-6)} ]],"Meters"]

(* ROZKŁAD PLANCKA W FUNKCJI CZĘSTOŚCI *)
PlanckNu[ν_,T_]:=(2*Pi*h*ν^3 )/( c^2 (Exp[( h*ν)/(kb*T)]-1));

Plot[{PlanckNu[ν,5700],PlanckNu[ν,5000],PlanckNu[ν,4500]},{ν,0,2*10^(15)},PlotRange->Full,AxesLabel->{"ν  [Hz]","(Subscript[e, ν]) [W/(m^3)]"},PlotLegends->{"T=5700 K","T=5000 K","T=4100 K"}]
(* częstość νmax odpowiadająca maksimum (Subscript[e, ν]) oraz odpowiadająca jej dł. fali *)
νmax=Quantity[x/.Last[FindMaximum[PlanckNu[x,5700]*10^12,{x,2*10^14} ]],"1/Seconds"]
v:=Quantity[c,"Meters/Seconds"]
v/νmax