Résolution par éléments finis de PDE (matlab)

Considérons le problème aux limite suivant :
$$\begin{cases}
-u''(x)+6 u(x)=(-4x^2-6)\exp( x^2),\quad x\in \Omega=]-1,1[\\

u(-1)=u(1)=0
\end{cases}

$$ Pour un maillage de quatre éléments je veux écrire le scripte matlab qui permet de donner le tableau des connectivités et le tableau des coordonnées, et aussi écrire un programme matlab d’éléments qui permet de résoudre ce problème ?

Je suis un débutant en matlab.
Merci de votre cooperation.
Cordialement.