急求duffing方程lyapunov指数程序

shmily137

function [dx]=duffing(t,x,F)
% global F;
w0=1.0;
gama=0.168;
dx=[x(2);F*sin(w0*t)-gama*x(2)+1/2*x(1)*(1-x(1)^2)];

84shmily987

不好编啊！