各位Matlab的高手，本人有一道疑难程序无法运行得到图型，希望各位高手帮帮忙，再次先谢谢了，程序如下：
clear
global a w wp p;
a=0.25;w=1;wp=0.6667;
p=input('0.69');
tspan=[0 600];
y0=[0 0.2 0]';
%[t,y]=ode23('xxdb',tspan,y0);
[t,y]=ode23('fcdo',tspan,y0);
subplot(2,1,1)
plot(t,y(:,1),'b');
hold on
y0=[0.01 0.2 0]';
%[t,y]=ode23('xxdb',tspan,y0);
[t,y]=ode('fcdo',tspan,y0);
plot(t,y(:,1),'k')
xlabel('时间');ylabel('角位移');
function yp=fcdo(t,y)
global a w wp p;
yp=[y(2)-2*a*y(2)-w*sin(y(1))+p*cos(y(3)) wp]';
function yp=xxdb(t,y)
global a w wp p;
yp=[y(2)-2*a*y(2)-w*(y(1))+p*cos(y(3)) wp]';