fun_liu.m

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

function dy=fun_liu(t,y)
dy=zeros(3,1);
dy(1)=10*(y(2)-y(1));
dy(2)=40*y(1)-y(1)*y(3);
dy(3)=-2.5*y(3)+4*y(1).^2;