This thesis reports onreview and research work carried out on the numerical analysis of elastomers. The two numerical techniques investigated for this purpose are the finite and boundary element methods. The finite element method is studied so that existing theory is used to develop a finite element code both to review the finite element method as applied to the stress analysis of elastomers and to provide a comparison of results and numerical approach with the boundary element method. The research work supported on in this thesis covers the application of the boundary element method to the stress analysis of elastomers. To this end a simplified regularization approach is discussed for the removal of strong and hypersingularities generated in the system on non-linear boundary integral equations. The necessary programming details for the implementation of the boundary element method are discussed based on the code developed for this research. Both the finite and boundary element codes developed for this research use the Mooney-Rivlin material model as the strain energy based constitutive stress strain function. For validation purposes four test cases are investigated. These are the uni-axial patch test, pressurized thick wall cylinder, centrifugal loading of a rotating disk and the J-Integral evaluation for a centrally cracked plate. For the patch test and pressurized cylinder, both plane stress and strain have been investigated. For the centrifugal loading and centrally cracked plate test cases only plane stress has been investigated. For each test case the equivalent results for an equivalent FEM program mesh have been presented. The test results included in this thesis prove that the FE and BE derivations detailed in this work are correct. Specifically the simplified domain integral singular and hyper-singular regularization approach was shown to lead to accurate results for the test cases detailed. Various algorithm findings specific to the BEM implementation of the theory are also discussed.