loreq47(t,u)

function udot = loreq47(t,u)
global U V t0 k1 k2 k3 k4 k21 k22 k23 k24 a m1 m2 m3 m4 mu d0 d2 m
udot =[-i.*((t0 +  a * (  u(5)- u(4))).*u(2)  +( t0 + a * ( u(4)- u(6))).*u(3)  +(U./2.*(conj(u(2)).*u(2)+conj(u(3)).*u(3))+mu).*u(1));
       -i.*((t0 +  a * (  u(6)- u(5))).*u(3)  +( t0 + a * ( u(5)- u(4))).*u(1)  +(U./2.*(conj(u(3)).*u(3)+conj(u(1)).*u(1))+mu).*u(2));
       -i.*((t0 +  a * (  u(4)- u(6))).*u(1)  +( t0 + a * ( u(6)- u(5))).*u(2)  +(U./2.*(conj(u(1)).*u(1)+conj(u(2)).*u(2))-mu).*u(3));   
 i.*( -k1*(u(5)-2.*u(4) +u(6)) + m.*u(4))-a.*(-u(2)+u(3)-u(1)+u(1));
 i.*( -k1*(u(6)-2.*u(5) +u(4)) + m.*u(5))-a.*(-u(3)+u(3)-u(1)+u(2));
 i.*( -k1*(u(4)-2.*u(6) +u(5)) + m.*u(6))-a.*(-u(3)-u(3)-u(1)+u(2));
];