This work aims to present the development of a three-dimensional numerical code for simulating fracturing and fragmentation processes in geomaterials. The development of such code is important due to the interest of several engineering fields in understanding such processes. The code uses an extrinsic cohesive zone model to characterize the fracturing process, which differentiates it from most similar codes. The use of such model eliminates the artificial compliance introduced by the initial slope of intrinsic cohesive models and helps to reduce the computational cost. At first, the details of the code development are presented. Then, the implementations made on the code are verified upon the simulation of two types of laboratory tests commonly employed in rock mechanics (Brazilian tensile strength test and uniaxial compressive strength test), in which the parameters used and the results for comparison were obtained from the literature. Next, the results exhibiting the fracture patterns for each type of test and its stress-displacement, or stress-strain, curves are displayed. The results obtained with the developed code show a good agreement with the values obtained from the literature, with variations of less than 10% in the strengths values. At last, the final considerations are presented.