Se presenta un método libre de mallas con derivadas difusas y estabilización por penalización. Un análisis de error para la aproximación de la solución de una ecuación elíptica general en múltiples dimensiones, con condiciones de frontera tipo Neumann es desarrollado. Resultados numéricos y teóricos muestran que el error de aproximación y la velocidad de convergencia son mejores que en el método de elementos difusos.