Added functions for expectation values based on density matrix
Rewrote aiTEBD file to work with class systems, for iMPS and for iTEBD algorithm