This study presents a collocation scheme using modified Lucas wavelets to approximate solutions for important Benjamina Bona Mohany equations applied in diverse physical applications. The proposed scheme transforms the problem into an algebraic system, and comparative analysis with existing methods demonstrates its effectiveness in providing accurate solutions for Benjamina Bona Mohany partial differential equations.