Abstract: | In this paper, we propose a new high accuracy numerical method of O(k2 + k2h2 + h4) based on off-step discretization for the solution of 3-space dimensional non-linear wave equation of the form utt = A(x,y,z,t)uxx + B(x,y,z,t)uyy + C(x,y,z,t)uzz + g(x,y,z,t,u,ux,uy,uz,ut), 0 < x,y,z < 1,t > 0 subject to given appropriate initial and Dirichlet boundary conditions, where k > 0 and h > 0 are mesh sizes in time and space directions respectively. We use only seven evaluations of the function g as compared to nine evaluations of the same function discussed in and . We describe the derivation procedure in details of the algorithm. The proposed numerical algorithm is directly applicable to wave equation in polar coordinates and we do not require any fictitious points to discretize the differential equation. The proposed method when applied to a telegraphic equation is also shown to be unconditionally stable. Comparative numerical results are provided to justify the usefulness of the proposed method. |