A mathematical procedure was developed to utilize the complementary energy method, by minimization, in order to obtain an approximate analytical solution to the 31) stress distributions in bonded interfaces of dissimilar materials. The stress solutions obtained predict the Stress jumps at the interfaces, which cannot be captured by the current FEA methods. As a novel method, the penalty function is used to enforce the displacement boundary conditions at the interfaces. Furthermore, the mathematical procedure developed enables the integration of different interfacial topographies into the solution procedure. In order to incorporate the effects of surface topography, the interface is expressed as a general surface in Cartesian coordinates, i.e. F(x, y, z) = 0. In this paper, the scarf interface problem, i.e. y = x/2 surface is considered for verification of the method by comparison with finite element analysis (FEA) results. Comparison of the results reveals our new mathematical procedure to be a promising and efficient method for optimizing interface topographies.