clc, clear all;
v0=8;
v=v0;
t=0;
x0=0;
x=x0;
X=x;
T=t;
dt=0.1;
while v>0
    U=[x;v];
    k1=dt*dUdt(t,U);
    k2=dt*dUdt(t+dt/2,U+k1/2);
    k3=dt*dUdt(t+dt/2,U+k2/2);
    k4=dt*dUdt(t+dt,U+k3);
    U=U+(k1+2*k2+2*k3+k4)/6;
    t=t+dt;
    x=U(1);
    v=U(2);
    X=[X x];
    T=[T t];
end
plot(T, X)
X(end)