# c) * c) ** Skriv et program som beregner intensiteten fra Plancks str?lingslov
# som funksjon av en innlest temperatur, og gir som resultat b?lgelengden med
# st?rst intensitet. Gj?r dette ved ? lage en l?kke der b?lgelengden endres
# logaritmisk fra 10-16 til 106 m med 100 tall per dekade, og lagre den til enhver
# tid st?rste intensitetens b?lgelengde.
# Sammenlign med resultatet fra programmet i b).

#Vi trenger mattepakken:
import math               # Importerer hele pakken

# Vi definerer konstanter og tilordner dem tabellverdier:
c = 2.99792E8
h = 6.626E-34
pi = 3.14159
k = 1.381E-23
e = 2.7182818

# Vi definerer variabler og tilordner dem default-verdier:
f = 1.0
T = 1.0
Lambda = 1.0
Lambdamax = 1.0
Intensitet = 1.0
Intensitetmax = 1.0
mangler_verdi = True

# Be om temperatur:
mangler_verdi = True
while mangler_verdi:
    try:
    	T = float(raw_input('Oppgi temperaturen i kelvin: '))
    	mangler_verdi = False
    except:
    	print "Du maa oppgi et tall!"

# Startverdier for l?kke:
Lambda = 1E-7
Lambdamax = Lambda
Intensitetmax = 0.0

# Gj?r l?kka:
while Lambda <= 1E6:
  f = c / Lambda
  Intensitet = 2*pi*h*f/(Lambda**2)/(e**(h*f/k/T)-1)
  if Intensitet > Intensitetmax: 
      Intensitetmax = Intensitet
      Lambdamax = Lambda
  Lambda = Lambda*1.01 # ?ker Lambda med med 1%
  # print Lambda, Intensitet

print
print 'Planck lambda max:  ',Lambdamax,'   Intensitet: ',Intensitetmax
print

# Wiens lov, fra b):
a = 2.9E-3
print 'Wien lambda max:  ',a / T