got muj working
added epsilon for zero concentrations
sized the outputs for chmatrix and pi derivatives in n2pi
complex now optional, for better performance
renamed n2pi to n2ls for n to linear system
computed an intermediate derivative
did s0 and cp0
Merge remote-tracking branch 'origin/master'

Conflicts:
	CEAFS.py
	test_CEA.py