Problemas con valores en la frontera con Python

Hacer un problema de ecuación con valores en la frontera de Neumann en Python.
-U’’(x) + ku(x)=f(x), con a<x<b
Si u’(a)=u’(b)=0

Hacer un esquema de diferencias finitas , sobre una partición uniforme de nodos, aproximando u’’ por la diferencia centrada y las condiciones de contorno, introducir los nodos fantasmas fuera del intervalo, la diferencia hacia adelante en b y hacía atrás en a.
Y en ambos casos representarlo de forma matricial