Removed duplicate plot of the initial cube
No need to multiply by an identity matix here
Changed the simulation time back to 0~1s
No need to have this be a matrix; is constant.
Clarified comments and fixed initial deviatoric stress

Initial deviatoric stress is zero, not an identity.
Reuse already established variable instead of recreating it
Rewritten for readability
Added relevant units to plots
Fixed E_e:IE_0:E_e calculation
Fixing theta_dot calculation; still doesn't work?
Fixing plastic strain; still broken?
Fixing densification; should be working
Optional debugging line; remove later