如何用matlab 拟合黑体fu辐出度曲线？

c1=93.61;c2=4.966;TR=3813;TB=7245;
G=inline('93.61./x.^5./(exp(4.966./x)-1)');
T=1500:100:6000; xr=T/TR; xb=T/TB;
for k=1:46; y (k)=quad(G,xb(k),xr(k)); end;
semilogy(T,y); hold on;
semilogy([3650 3650],[10^-3 1],'g'); semilogy(2780,0.083,'b*');
grid on; xlabel('T/K'); ylabel('Y/%')

黑体辐出度曲线？？