% The Stefan-Boltzmann constant 
function val = sig( )     % val=W/m^2/K^4
  val = 5.6704e-08;
end
