yxt2dtQ1(t)Q2(t)=13a1a2RD(U2+a1b2,U2+a2b1,U2)+(xy/U),