Local matrix calculation worked on.  Added Planck integration functions