%%
%Uppgift 4 IV)
%Beräknar deltaTMax på samma sätt som i Uppgift 4 I)
%Får värdet deltaTMax = 1.1025e-04
%Upprepar beräkningen i uppgift 4 III) med deltaT = 0.1*deltaTmax

clc; clear all; close all;
k1 = 5300; k2 = 136000;
deltaTmax = 1.1025e-04;
v = 65/3.6;
L = 1;
t0 = 0;
t1 = 2;
T = (L/v)*30;
vVec = [0; 0; 0; 0]; %Försäkra oss om korrekt storlek på vVec
h = 0.1*deltaTmax;
n = T/h;
for i = 1:n
    t = (i-1)*h;
    vVec(:,i+1) = vVec(:,i) + h*quartercar(t,vVec(:,i),k1,k2); %Euler
end

tid = 0:h:T;
plot(tid,vVec(1,:))
hold on;
plot(tid,vVec(2,:));
xlabel('Tid [s]')
ylabel('Förskjutning [m]');
legend('Z1', 'Z2');