Thermal conductivity problem...

Having a question regarding using band matrixes. We basically get no information in the document so it's hard to understand what to do.

g0 and gl is the temperature in the edges of the rod, c is the thermal diffusivity and f(x) is the heat source.

Let c=2, L=1 g0=40, gl=60 f(x) = 200exp(−(x-L/2)^2) Use n=30

c=2; L=1; g0=40; gl=60; h=L/n+1;
f=@(x)200.*exp(-(x-(L./2)).^2);
n=30; ett=ones(n,1);
A=spdiags([-ett 2*ett -ett],[-1 0 1],n,n);
b=[(h^2/c)*f(x)+g0;(h^2/c)*f(x)]

Solve the problem and sketch the solution.

Above is how far I've gotten.