In this work, an algorithm based on method of lines coupled with radial basis functions namely meshless method of lines (MMOL) is presented for the numerical solution of Kawahara, modified Kawahara and KdV Kawahara equations. The motion of a single solitary wave, interaction of two and three solitons and the phenomena of wave generation is discussed. The results are compared with the exact solution and with the results in the relevant literature to show the efficiency of the method.