T=0:0.1:10;syms h,c,k,E,b,R,G,I(注意：上面的几个变量需要你自己输入给定的值,我这个是按照变量处理的)F=1/(2*pi)*(h*c/kT.*exp(E/(b*R*T.))*exp(G/(R*T.))+pi/8*(I/(k*T.))^0.5)^(-1);plot(T,F)
