Proponemos un método de cuadratura diferencial (DQM) basado en funciones de base de B-spline hiperbólicas cúbicas para computar ecuaciones de onda en 3D. Este método convierte el problema en un sistema de EDOs. Utilizamos un esquema óptimo de Runge-Kutta de cinco etapas y orden cuatro (SSPRK-(5,4)) para resolver el sistema de EDOs obtenido. También se investiga el análisis de estabilidad de la matriz. La precisión y eficiencia del método propuesto se demuestran a través de tres ejemplos numéricos. Se ha encontrado que el método propuesto proporciona resultados más precisos que los métodos existentes. El propósito principal de este trabajo es presentar una técnica precisa, económica, fácil de implementar y estable para resolver ecuaciones diferenciales parciales hiperbólicas.