Simulation of fourth order Bi-harmonic equation