Return to Answer

function times = f()
  lambda = 1.5
  list times 
  t = 0
  do 
    dt = random_exponential(lambda)
    t = t + dt
    st = scale(t)
    if st<3
      times.add(st)
    end if
  while st < 3
end

function st = scale(t)
  % scale time according to inteinsities
  scale2 = 2.1/1.5
  scale3 = 3.4/1.5
  if t < 1
    st = t
  else if t < 1 + scale2
    st = 1 + (t-1) / scale2 
  else
    st = 1 + 1 / scale2  + (t - (1 + scale2)) / scale3
  end if
end

function times = f()
  lambda = 1.5
  list times 
  t = 0
  do 
    dt = random_exponential(lambda)
    t = t + dt
    st = compress(t)
    if st<3
      times.add(st)
    end if
  while st < 3
end

% compress slow time to faster times after 1st hour
function st = compress(t)
  scale2 = 2.1/1.5
  scale3 = 3.4/1.5
  t2 = 1 + scale2;
  st = min(t,1) + (min(t,t2)-min(t,1))/scale2 + (max(t,t2)-t2)/scale3
end

function times = f()
  lambda = 1.5
  list times 
  t = 0
  do 
    dt = random_exponential(lambda)
    t1 = t + dt
    if t1 < 1
      t = t1
    else if t1 < 2
      dt = dt * lambda / 2.1
      if t < 1
        dt = (1 - t) * (1 - lambda / 2.1)
      end if
      t = t + dt
    else
      dt = dt * lambda / 3.4
      if t < 2
        dt = (2 - t) * (1 - lambda / 3.4)
      end if
      t = t + dt
    end if
    if t<3
      times.add(t)
    end if
  while t < 3
end

function times = f()
  lambda = 1.5
  list times 
  t = 0
  do 
    dt = random_exponential(lambda)
    t = t + dt
    st = scale(t)
    if st<3
      times.add(st)
    end if
  while st < 3
end

function st = scale(t)
  % scale time according to inteinsities
  scale2 = 2.1/1.5
  scale3 = 3.4/1.5
  if t < 1
    st = t
  else if t < 1 + scale2
    st = 1 + (t-1) / scale2 
  else
    st = 1 + 1 / scale2  + (t - (1 + scale2)) / scale3
  end if
end

function times = f()
  lambda = 1.5
  list times 
  t = 0
  do 
    dt = random_exponential(lambda)
    t1 = t + dt
    if t1 < 1
      t = t1
    else if t1 < 2
      t = t + dt * lambda / 2.1
    else
      t = t + dt * lambda / 3.4
    end if
    if t<3
      times.add(t)
    end if
  while t < 3
end

function times = f()
  lambda = 1.5
  list times 
  t = 0
  do 
    dt = random_exponential(lambda)
    t1 = t + dt
    if t1 < 1
      t = t1
    else if t1 < 2
      dt = dt * lambda / 2.1
      if t < 1
        dt = (1 - t) * (1 - lambda / 2.1)
      end if
      t = t + dt
    else
      dt = dt * lambda / 3.4
      if t < 2
        dt = (2 - t) * (1 - lambda / 3.4)
      end if
      t = t + dt
    end if
    if t<3
      times.add(t)
    end if
  while t < 3
end