The authors present an algorithm for determining the stiffness of the bone tissue for individual ranges of bone density. The paper begins with the preparation and appropriate mechanical processing of samples from the bovine femur and their imaging using computed tomography and then processing DICOM files in the MIMICS system. During the processing of DICOM files, particular emphasis was placed on defining basic planes along the sides of the samples, which improved the representation of sample geometry in the models. The MIMICS system transformed DICOM images into voxel models from which the whole bone FE model was built in the next step. A single voxel represents the averaged density of the real sample in a very small finite volume. In the numerical model, it is represented by the HEX8 element, which is a cube. All voxels were divided into groups that were assigned average equivalent densities. Then, the previously prepared samples were loaded to failure in a three-point bending test. The force waveforms as a function of the deflection of samples were obtained, based on which the global stiffness of the entire sample was determined. To determine the stiffness of each averaged voxel density value, the authors used advanced optimization analyses, during which numerical analyses were carried out simultaneously, independently mapping six experimental tests. Ultimately, the use of genetic algorithms made it possible to select a set of stiffness parameters for which the error of mapping the global stiffness for all samples was the smallest. The discrepancies obtained were less than 5%, which the authors considered satisfactory by the authors for such a heterogeneous medium and for samples collected from different parts of the bone. Finally, the determined data were validated for the sample that was not involved in the correlation of material parameters. The stiffness was 7% lower than in the experimental test.