用matlab怎么写这个程序呀

这是典型的常微分方程，可以使用ode45求出数值解。

m = 1;
g = 9.8;
dfun = @(t,x) [x(2); -g-0.1*abs(x(2))*sign(x(2))/m];
[T,X] = ode45(dfun,[0 100],[0 0]);
figure
plot(T,X(:,2))
xlabel('时间/s')
ylabel('速度m/s')

结果：

速度最终趋于稳定，此时重力与空气阻力相等。