matlab解一级倒立摆的微分方程组

[行走在时间里]

一级倒立摆的微分方程组为：
x'=x'
x''=-0.088x'+0.629y+0.88u
y'=y'
y''=-0.2356x'+27.825y+2.35u

我是这样写的，但是发现运行的结果不对
function yuan
Y0 = [0 0 0 0];

% 解方程
[t,Y]= ode45(@ode,[0 10],Y0);
t1=Y(:,1);
t2=Y(:,3);

% 绘图
subplot 211
plot(t,t1);
subplot 212
plot(t,t2);

% 微分方程定义
function dY = ode(u, Y)
y1=Y(2);
y2=Y(3);
dY = [
Y(2);
-0.088316*y1+0.629317*y2+0.88*u;
Y(4);
-0.235655*y1+27.8285*y2+2.35*u;
];
初学matlab, 在网上看了一个例子开始做，发现网上的例子是二元二阶的，而这里却多了一个参数u 不知该如何处理？