from datetime import timedelta
import pandas as pd
import numpy as np

A = np.array([0,0,0,1,0,0,1,0,1,0,0,1,1,1,0,1,0,1,0,1,0,1,0,0,1,0,1,0,1,0])
B = np.array([0,0,1,1,0,0,1,0,1,1,0,1,1,0,1,1,0,1,1,0,1,1,1,1,0,1,1,1,1,0])

range = 100 minuts.
λ_A = 4 per 100 mints.
λ_B = 5 per 100 mints.

d = {'Rest_A':A, 'Rest_B':B}
Time_Interval = pd.timedelta_range('11:00:00', periods=30, freq='10T')
df = pd.DataFrame(d, Time_Interval)

period_one = poisson.pmf(0, λ_A/30)  = 0.875173
period_two = poisson.pmf(0, λ_A*2/30) = 0.765928
period_three = poisson.pmf(0, λ_A*3/30) = 0.670320

period_four = P[X4 = 1| X3 = 0] = P[X4 = 1] = poisson.pmf(1, λ_A*4/30) = 0.312878

period_five = P[X5 = 0 | X4 = 1] = P[(X5 = 0) n (X4 = 1)] / P[X4 = 1] =
              P[X5 = 0] / P[X4 = 1] = poisson.pmf(1, λ_A*4/30) / poisson.pmf(0, λ_A*5/30)
              = 1.640950 ??

period_six = P[X6 = 1 | X5 = 0] =
             P[(X6 = 1) n (X5 = 0)] / P[X5 = 0] = P[X5 = 0] / P[X5 = 0] = 1 ???

P[X6 = 1 | X5 = 0, X4 = 1, X3 = 0, X2 = 0, X1 = 0] ??