We investigate the existence and uniqueness of positive solutions to an integral equation involving convex or concave nonlinearities. A numerical algorithm based on Picard iterations is provided to obtain an approximation of the unique solution. The main tools used in this work are based on partial-ordering methods and fixed-point theory. Our results are supported by examples.