In this study, the beta function that is encountered in computational mathematics and physics is analyzed. The correct evaluation of this function also affects the accuracy of other mathematical functions in quantum mechanical calculations. Especially in recent years, there is an interest in studies related to the beta function for zero and negative p and q integers. This study, considering the neutrix limits of the beta function, presents new relations for the numerical computation of the beta function, especially for negative integers p and q. In addition, taking into account the definition of the beta function for positive p and q integer values, an algorithm is created to calculate the function for all integer values. Finally, numerical results obtained with the help of our new recurrence relations and algorithm are presented.