O teste de Dunnett é um teste de comparações múltiplas com um tratamento testemunha que controla simultaneamente a taxa de erro tipo I por experimento. A limitação para o seu uso é a dificuldade de obter as probabilidades da distribuição e os valores dos quantis da estatística do teste, pois as correlações possíveis entre os tratamentos têm larga amplitude. Neste trabalho são propostas funções R para obtenção das probabilidades do teste bilateral de Dunnett balanceado, usando quadraturas gaussianas para a resolução das integrais. Quando se comparam os resultados destas funções com aqueles obtidos pelas rotinas dos pacotes mvtnorm e mvtnormpcs do software R e com os apresentados na literatura observam-se resultados precisos.