空气涡轮机的matlab程序模拟

空气涡轮对阶跃压力调节器输入的响应。这是从一本外文文献中看到的，书中给出的源代码如下：
if (nargin<1)
kP=1
kT=2
tauP=10
tauT=40
c=eye(2)
b=[kP/tauP;0]
a=[-1/tauP,0;kT/tauT,-1/tauT]
xDot=struct('a',a,'b',b,'c',c,'u',0)
if (nargout==0)
disp('RHSAirTurbine struct:')
end
return
end
xDot=d.ax+d.bd.u
第二部分：
tEnd=1000
d=RHSAirturbine
d.u=100
dT=0.02
n=ceil(tEnd/dT)
x=[0;0]

xP=zeros(2,n)
t=(0:n-1)*dT
for k=1:n
xP(:,k)=x
x=RungeKutta(@RHSAirTurbine,t(k),x,dT,d)
end

[t,tL]=TimeLabel(t)
yL={'p (N/m^2)' '\omega (rad/s)'}
tTL='Air Turbine Simulation'
PlotSet(t,xP,'x label',tL,'y label',yL,'figure title',tTL)
程序实现后，得到图像如下

完全按照这个来，却不能实现。请问是哪里出问题了？