Standard methods for performing analytic perturbative calculations for the process of e⁺e^- → q? up to ?(α_s) are explained and results given. An emphasis is given to the organisation calculations using the Cutkosky cutting rules and the renormalisation of the massive quark propagator. Methods for numerical integration are presented including those used in VEGAS. The numerical methods used in the Beowulf program for calculating infra-red safe observables for jet events from electron-positron collisions are also explained. Cancellations of singularities required for numerical calculations are demonstrated using an example in ?³ theory both numerically and graphically. Renormalisation by subtraction of appropriate integrals is also covered. Adaptations of the Beowulf procedure required for the inclusion of massive fermions are developed and explained. An alternative method for including the quark self energy and its related cuts using scalar decomposition, numerically equivalent integrals and its spinor structure is introduced. The methods are used to calculate the ?(α_s) corrections to the process e⁺e^- → q? using VEGAS. Drawbacks of the smearing function required in the numerical integration due to the corrections dependence on the mass and centre of mass energy are discussed. Results of the ?(α_s) cross section using the numerical method verify the procedure. The method was then used to see the effects of mass on the thrust distribution and when using the Durham and JADE jet algorithms.