function OUT=new(t,X)

a=0.04;b=2;c=0.16;

x(1)=X(1);x(2)=X(2);x(3)=X(3);x(4)=X(4);

Q=[ X(5),  X(9),  X(13), X(17);
    X(6),  X(10), X(14), X(18);
    X(7),  X(11), X(15), X(19);
    X(8),  X(12), X(16), X(20);];
    
    
K=(abs(x(2)+x(3)+1)-abs(x(2)+x(3)-1))-(abs(x(2)+x(3)+3)-abs(x(2)+x(3)-3))-0.25*(abs(x(2)+x(3)+5)-abs(x(2)+x(3)-5));    

dx(1)=-a*x(1)+b*x(2)+K*b;
dx(2)=-b*x(1)+c*x(2)+c*K;
dx(3)=-c*x(3)-x(4)-c*K;
dx(4)=x(3)-c*x(4)+K;





DX1=[dx(1);dx(2);dx(3);dx(4)];	

K=(abs(x(2)+x(3)+1)-abs(x(2)+x(3)-1))-(abs(x(2)+x(3)+3)-abs(x(2)+x(3)-3))-0.25*(abs(x(2)+x(3)+5)-abs(x(2)+x(3)-5));     

J=[  -a b+b*K b*K 0;
      -b c+c*K c*K 0;
       0 -c*K -c-c*K -1;
       0 K 1 -c];
  
F=J*Q;

OUT=[DX1; F(:)];