La ecuación de Poisson surge con frecuencia en muchos campos de la ciencia y la ingeniería. Dado que las soluciones exactas rara vez son posibles, los enfoques numéricos son de gran interés. A pesar de esto, puede ser difícil encontrar una discusión concisa de un enfoque sistemático para construir un solucionador numérico de Poisson flexible y general. En este artículo introductorio, se presenta una discusión exhaustiva sobre cómo construir un solucionador de matriz de diferencias finitas que puede resolver la ecuación de Poisson para geometrías y condiciones de contorno arbitrarias. Las condiciones de contorno se implementan de manera sistemática, lo que permite la modificación sencilla del solucionador para diferentes problemas. También se discute un enfoque de definición de geometría basado en imágenes. El código Python de la receta numérica está disponible públicamente. Se presentan ejemplos numéricos que muestran cómo configurar el solucionador para diferentes problemas.