function f=liutrack3(t,X)%liu system's synchronization with new attractor
global p;
f=zeros(6,1);
x1=X(1);x2=X(2);x3=X(3);y1=X(4);y2=X(5);y3=X(6);
f(1)=10*(x2-x1);
f(2)=40*x1-x1*x3-(49+0.1*p)*x1+9*x2+x1*x3+0.1*p*y1-0.4*y3+0.4*y1*y2-0.15*y2*y3-0.1*y1*y3*y3+0.4*y1*y2*y2;
f(3)=-2.5*x3+4*x1*x1;
f(4)=y3-y1+y2*y3;
f(5)=2.5*y2-y1*y3;
f(6)=4*y1*y2-4*y3;