MATLAB ode45 杜芬方程的相图 哪里出错了？

omega是个啥？常数还是变量？

clear all
clc

gamma=0.1;
kapa=1;
xi=1;
F=0.1;
omega=5;

f=@(t,x)([x(2);-gamma*x(2)+kapa*x(1)-xi*x(1)^3+F*cos(x(3));omega]);
[t,X]=ode45(f,0:.1:4,[0 .1 0]);
plot(X(:,1),