Morten's code added
New matrix elements