from CieXYZ import *
cie=CIEXYZ("../ciexyz/ciexyz31.txt")
lumen=[cie.ener2lum(357,700),cie.ener2lum(6.83,546.1),cie.ener2lum(4.95,435.8)]
            
 
r=lumen[0]
for k in range(3):
    lumen[k]=lumen[k]/r
            
from pylab import *
[r,g,b]=cie.rgbColorim()
clf()
xlabel('lambda (nm)')
plot(cie.Lambda,r,c='r',label='rc')
plot(cie.Lambda,g,c='g',label='gc')
plot(cie.Lambda,b,c='b',label='bc')
title(u'Fonctions colorim\xe9triques')
legend()
grid(True)