function y=f(t,y)

m=70;
k=20;
a1=1;   a2=1;
g=10;
y0=[-30;0];

x=y(1);       %  x代表蹦極者的位移（在以橋梁面為基準面、以指向地面方向為正方向的坐標系中）
x_dot=y(2);   %  x_dot 代表 x 的速度
if x>0
    y=[0,1;-k/m,-a1/m-(a2/m)*abs(x_dot)]*y+[0;g];
else
    y=[0,1;0,-a1/m-(a2/m)*abs(x_dot)]*y+[0;g];
end