Numerical methods for multiserver discrete time queues with batch markovian arrivals