function positions = createPositionsLin(TimeSteps, TotalShift)  
step       = ((TimeSteps)/TotalShift)^-1;
positions2 = round([0:step:TotalShift; 0:step:TotalShift]');
positions  = positions2(1:TimeSteps, 1:2);

