function xdot=dlorenz(t,x);

global r b s c N


u=x(1);v=x(2);w=x(3);
u1=x(4);v1=x(5);w1=x(6);

xdot=zeros(N,1);

xdot(1)=s*(v-u)-c(1)*(u-u1);
xdot(2)=r*u-v-u*w-c(2)*(v-v1);
xdot(3)=u*v-b*w-c(3)*(w-w1);



xdot(4)=s*(v1-u1)+c(1)*(u-u1);
xdot(5)=r*u1-v1-u1*w1+c(2)*(v-v1);
xdot(6)=u1*v1-b*w1+c(3)*(w-w1);