积分后，可以用vpa（）函数求出其数值解。
注：h为普朗克常数，c为光速，k为玻尔兹曼常数没有考虑其单位的统一，请你自行修改
h = 6.626196*10^-34; %h为普朗克常数
c=299792458; %c为光速
k = 1.3806505*10^-23; %k为玻尔兹曼常数
A=h*c/k;
int(x^(-5)/(exp(A/(400*x))-1),x,550,625)
vpa(ans)
