%普朗克黑體單色輻射斯特潘常數(shù)和維恩常數(shù)求法
clear                                  %清除變量
k=1.38054e-23;                         %玻爾茲曼常數(shù)
h=6.626e-34;                           %普朗克常數(shù)
c=2.997925e8;                          %光速
syms x                                 %定義符號(hào)變量
y=x^3/(exp(x)-1);                      %被積函數(shù)
i=int(y,0,inf)                         %求積分
sigma=eval(2*pi*k^4/h^3/c^2*i)         %求斯特潘常數(shù)
y=x^5/(exp(x)-1);                      %普朗克約化公式
d=diff(y)                              %求符號(hào)導(dǎo)數(shù)
s=solve(d)                             %求符號(hào)零解
e=eval(s)                              %求零解的數(shù)值
b=h*c/e/k                              %求維恩常數(shù)