Se presentó un enfoque de diferencias finitas con mallas no uniformes para simular respuestas magnetotelúricas en estructuras 2D. Presentamos la fórmula de cálculo de este esquema a partir del problema de valor límite del campo eléctrico y campo magnético, y comparamos las soluciones de diferencias finitas con los resultados numéricos de elementos finitos y soluciones analíticas de un modelo 1D. En primer lugar, se probó un modelo de semiespacio homogéneo y el enfoque de diferencias finitas puede proporcionar una precisión muy buena para la modelización magnetotelúrica 2D. Luego las comparamos con las soluciones analíticas para el modelo geoeléctrico de dos capas; los errores relativos de la resistividad aparente y la fase de impedancia aumentaron cuando se incrementó la frecuencia. Para concluir, comparamos nuestros resultados de simulación de diferencias finitas con el modelo COMMEMI 2D-0 con las soluciones de elementos finitos. Ambos resultados están en estrecho acuerdo entre sí. Estas comparaciones pueden confirmar la validez y confiabilidad de nuestro algoritmo de diferencias finitas. Además, un proyecto futuro ampliará las estructuras 2D a 3D, donde las mallas no uniformes deberían funcionar especialmente bien.