% Svartkroppsstrålning enligt Planck
  clear, clf
  alfa=3.7415e-16; beta=0.014388;
  ds=1e-8;  s=ds*(1:250)';  
  for T=3000:1000:6000
    p=beta./(s*T);
    w=alfa./(s.^5.*(exp(p)-1)); plot(s,w), hold on
  end
  disp('Klicka in textplaceringarna!')
  for i=3:6, gtext(['T=' int2str(i*1000)]), end
  
  
  
