Método matricial para la solución de EDs