#*************************************************************************
# programme de tracé de la loi de Planck
# Dominique Lefebvre  février 2013
# www.tangenteX.com
#*************************************************************************

# initialisation des variables
assume(T>0):

# définition des constantes
h := 6.63e-34:
kb := 1.38e-23: 
c := 3.00e8:

# définition de la loi de Planck
f := (lambda,T)->8*Pi*(h*c/(lambda^5))*1/(exp(h*c/(lambda*kb*T)) - 1);

# calcul des courbes
for k from 0 to 6 do
  p||k := plot(f(lambda,1000*(1 + k*0.5)), lambda=1e-7..3e-6, 'f'=0..2e5):
od:

#tracé des courbes
with(plots):
display([p||(0..6)]);