implemented asymmetry = epsi / pol in matrix code