liu_fun.m

根据《物理学报》2006年第10期上的4篇论文

function dy = Liu_fun(t,y)
dy = zeros(6,1);
a=1.0;b=40;c=2.5;h=4;k=3;
dy(1) = a*(y(2) - y(1));
dy(2) = b*y(1) - y(1)*y(3);
dy(3) = -c*y(3) + h*y(1)*y(1);
dy(4) = a*(y(5) - y(4)) + k*(y(1)-y(4));
dy(5) = b*y(4) - y(4)*y(6) + k*(y(2)-y(5));
dy(6) = -c*y(6) + h*y(4)*y(4) + k*(y(3)-y(6));