Results involving the approximation of the difference between two integral means are utilised to obtain bounds on the integral mean divergence and the f-divergence due to Csiszár. The current work does not restrict the functions involved to be convex. If convexity is imposed then the integral mean divergence is the Hermite-Hadamard divergence introduced by Shioya and Da-te.