In this paper, we report three level implicit method of high accuracyschemes for the numerical solution of fourth order nonhomogeneousparabolic partial dierential equation, that governs the behavior of avibrating beam. Parametric septic spline is used in space and nitedierence discretization in time. The linear stability of the presentedmethod is investigated. The computed results for three examples arecompared wherever possible with those already available in literaturewhich shows the superiority of the proposed method.