A=2;
B=1;
fun=@(t,x) [x(2);A-B/x(1)^3];
ts=[0,10];
x0=[1 0]; %x 和 dx/dt的初值
[t x]=ode45(fun,ts,x0);
plot(t,x(:,1));
xlabel('t');
ylabel('x');
