dr11_dt = gt - 2*(gp + gt)*r11 - E1*i*r13 + E1c*i*r31 + G*r33;
dr12_dt = (-gp - 2*gt - d1*i + da*i)*r12 - i*(E2*d(t) + E3*dc(t))*r13 + E1c*i*r32;
dr13_dt = -(E1c*i*r11) - i*(E3c*d(t) + E2c*dc(t))*r12 + (-G - gp - 2*gt - d1*i)*r13 + E1c*i*r33;
dr22_dt = gt - 2*gt*r22 - i*(E2*d(t) + E3*dc(t))*r23 + i*(E3c*d(t) + E2c*dc(t))*r32 + G*r33;
dr23_dt = -(E1c*i*r21) - i*(E3c*d(t) + E2c*dc(t))*r22 + (-G - 2*gt - da*i)*r23 + i*(E3c*d(t) + E2c*dc(t))*r33;
dr33_dt = 2*gp*r11 + E1*i*r13 + i*(E2*d(t) + E3*dc(t))*r23 - E1c*i*r31 - i*(E3c*d(t) + E2c*dc(t))*r32 - 2*(G + gt)*r33;