[t,z]=ode45(@zdot,[0,2], [0,0,0,0]);
plot(t,z(:,1), 'r')
hold on
plot(t,z(:,2), 'g')
xlabel('tid [s]')
ylabel('höjd [m]')
legend('x1(t)', 'x2(t))')