import SimPy.MarkovClasses as Cls

probMatrix = [
    [0.721,	0.202,	0.067,	0.010],
    [0.000,	0.581,	0.407,	0.012],
    [0.000,	0.000,	0.750,	0.250],
    [0.000,	0.000,	0.000,	1.000]
]

# discrete to continuous
transRateMatrix = Cls.discrete_to_continuous(
    trans_prob_matrix=probMatrix, delta_t=1)
print(transRateMatrix, '\n')

# continuous to discrete
newProbMatrix, prob2events = \
    Cls.continuous_to_discrete(trans_rate_matrix=transRateMatrix, delta_t=1)
print(newProbMatrix)
print('Upper bound for the probability of 2 transitions withing delta_t:', prob2events)


