You guys might find this article interesting.
It presents a linear algebra problem, then writes a program to solve it in Maple, Mathematica, and Maxima to see which one is better. What do you think?
Here is an example of manipulating an Array of pixels. I chose the x-rite ColorChecker as a model so there would be published results to check my work. A number of details about color spaces have become clear through this exercise. The color adaptation process was modeled by converting betweenXYZ and LMS. Different black points may be selected depending on how close to zero illuminance one would accept as a good model.
I look forward to extending this work to verify and improve the color calibration of my photography. Also some experimentation with demosaicing should be possible.