Compute Black body illuminants

Using the black body radiation concept, we can compute XYZ of all black-body illuminant. Prediction of this illuminants seems to be very close to CIE recommended illuminants such as CIE illuminant A

Let's show a quick example. We know illuminant A 2856K has 10 degree observer based X =111.144, Y=100,  Z=35.200

Now, lets compute using our black body illuminant formula

get_XYZ_of_blackbody(2856).*100

ans =

  111.1385  100.0000   35.2127

Pretty close! :)

Have fun!