duffing.m.txt

我编写的关于求解达芬方程的程序,各位请指教一下？

function dx=duffing(t,x) 
% 定义达芬Duffing方程： 
%  diff(x,2)+mu*diff(x,1)-x+x^3=F*cos(w*t) 
mu=0; 
F=0; 
w=1; 
dx=[x(2);F*cos(w*t)+0.1*cos(w*t)-mu*x(2)+x(1)-x(1).^3];