In this work we approximate the solution of initialboun-dary value problem using a Galerkin-finite element method for the spatial discretization and Implicit Runge-Kutta method for the spatial discretization and implicit Runge-Kutta methods for the time stepping. To deal with the nonlinear term f(x, t, u), we introduce the well-known extrapolation sheme which was used widely to prove the convergence in $L_2$-norm. We present computational results showing that the optimal order of convergence arising under $L_2$-norm will be preserved in $L_{infty}$-norm.