Este documento trata sobre la solución numérica de un problema de valores en la frontera no local para una ecuación pseudoparabólica bidimensional que surge en muchos fenómenos físicos. Se investiga un método implícito de dirección alternante de tres capas para la solución de este problema. Este método generaliza el método ADI de Peaceman-Rachford para la ecuación parabólica 2D. La estabilidad del método propuesto se demuestra en una norma especial. Investigamos un problema de autovalores algebraicos con matrices no simétricas para demostrar esta estabilidad. Se presentan resultados numéricos.