Implementation of a triangular probabilistic distribution for optimal parametrization of the SEIR model recovery rates with delay