split up init an matrix methods
check
made T and P attrs
switched to residual formulation, but not seperate function
cleaned up a couple more nested loops