function xd = my_fun(t,x)

xd1 = x(2);
xd2 = -9.81;

xd = [xd1; xd2]