%
% Spektrales Emissionsvermögen El(T) 
%   Kurve El(T) für Temperatur T = 6000 K  berechnen und 
%   in die Excel-Datei mit dem Namen xlsFile schreiben.
%   In Spalte A die lambda-Werte von 50 nm bis 3000 nm, in 50 nm-Schritten
%   in Spalte B die zugehörigen E-Werte für T = 6000 K, 
%
function xlsPlanck( xlsFile )

   c = 2.998e8;    % m/s  
   h = 6.626e-34;  %  Js   
   k = 1.38e-23;   % J/K 
   hc_k = h*c/k;  

   lambda = ( 50:50:3000 ) * 10^(-9);  % m
   lambda = lambda' ;  % transponierter Vektor 
   
   T = 6000;
   E = 2*pi*h*c^2 ./ ( lambda.^5 .* ( exp( hc_k./(T*lambda)) - 1 ) );  
   
   xlswrite( xlsFile, lambda, 1, 'A' ); 
   xlswrite( xlsFile, E, 1, 'B' );
   
   plot( lambda, E ); 
 