duffing混沌振子形式如下：

k,a,c,f为自定义系数，将初值设为,k=0.5,a=c=1

此时可通过更改f的值从0到1来改变duffing混沌系统状态，从固定点状态，小周期状态，混沌状态到大周期状态。例如f=0.6时处于混沌状态，如下：

上图横坐标x，纵坐标x的一阶导数

上图横坐标为时间t，纵坐标x

clc;clear;
[x,y]=ode113(@duffingFunc,[0:0.05:200],[0;1]);
figure(1);
plot(x,y(:,1));
grid on;
figure(2);
plot(y(:,1),y(:,2));
grid on;
function res=duffingFunc(x,y)
    k=0.5;
    a=1;
    c=1;
    f=0.75;
    res=zeros(2,1);
    res(1)=y(2);
    res(2)=a*y(1)-c*y(1)^3-k*y(2)+f*cos(x);
end

求解利用matlab中的ode113函数，duff ingFunc函数步骤为