%
%  bbdemo.m
%  Demonstration of bbspec
%                   bbody
%                   bbd
%
%
disp(' Demonstration of bbspec');
[t,lambda]=meshgrid([100,300,500,1000,3000,5000,10000],[0.01:0.01:30]);
sre=bbspec(t,lambda);
figure;loglog(lambda,sre);axis([0,30,1e-1,1e10]);
legend(num2str(t(1,:)'),2);
xlabel('\lambda, Wavelength, \mu m');
ylabel('m_\lambda, Spect. Radiant Exitance, W/m^2/\mu m');
title('Demonstration of bbspec');1

disp('Calling... [m,frac,sre,lpeak]=bbody(3000,0.4,0.7);');
[m,frac,sre,lpeak]=bbody(3000,0.4,0.7)

disp('Calling... bbd(3000,0.4,0.7);');
bbd(3000,0.4,0.7);