% Startvektor in Situation 1
N=500;
for i=1:2*N
    if i==1
        v0(i)=0.324;
    elseif i>N
        v0(i)=1;
    else v0(i)=0;
    end
end
%========================================================================
% Zeitspanne, in der berechnet werden soll
tspan=[0:100:15000];
% Aufrufen der Berechnung
tic
[t,v]=ode45('vdot1',tspan,v0)
%plot(t,v(:,1),'r-',t,v(:,N+1),'r.',t,v(:,50),'b-',t,v(:,N+50),'b.',t,v(:,100),'g-',t,v(:,N+100),'g.',t,v(:,N),'k-',t,v(:,N),'k.')
plot(t,v(:,1),t,v(:,3),t,v(:,5),t,v(:,7),t,v(:,10),t,v(:,12),t,v(:,15),t,v(:,20),t,v(:,25),t,v(:,30))
%plot(t,v(:,N+1),'.',t,v(:,N+3),'-',t,v(:,N+5),'.',t,v(:,N+7),'-',t,v(:,N+10),'.',t,v(:,N+12),'-',t,v(:,N+15),'.')
toc