Finite difference method on boundary value problem

I want to solve boundary value problem in n points, $\displaystyle x_{n}\in(0,1)$:

$\displaystyle -u'(x)(ku'(x))=12x^{2}$

where:

$\displaystyle u(x)\in[0,1]$

$\displaystyle u(0)=0$

$\displaystyle ku'(1)+hu'(1)=hu_{z}$

$\displaystyle x\in R$

with finite difference method instead of finite element method.

I substituted finite difference:

$\displaystyle u'(x)=\frac{u(x_{i-1})-u(x_{i+1})}{2h}$

$\displaystyle h=\frac{1}{n}$

and obtained:

$\displaystyle -kn^{2}[u(x_{i-1})-u(x_{i+1})]^{2}=3x_{i}^{2}$

$\displaystyle i=1,...,n$

And I don't know how to turn that into the matrix equation. Is there a way to do it?