%exercise 7.9
%Runge-Kutta
v0=8;

%startvärde för u
u1=0;
u2=v0;
u=[u1 u2];

%startvärde för x
t=0;

U=u;
T=t;
h=0.5;
while u(2)>0
   f1=f79(u);
   f2=f79(u+h*f1/2);
   f3=f79(u+h*f2/2);
   f4=f79(u+h*f3);
   u=u+(h/6)*(f1+2*f2+2*f3+f4);
   U=[U; u];
   T=[T ;t];
   t=t+h;
end
U
t
plot(T, U(:,1), T, U(:,1),'o',T, U(:,2), '--',T, U(:,2), 'o'), text(2.5, 4, 'xprim(t)'), text(2, 14 ,'x(t)')

