clear
m1= 5;
m2= m1 + 1;
m = m1 - 1;
h=1/m1;
x=linspace(0,1,m2)';
kfun=@(x) 5*x;

%construct the symmetric system (2.72)
xmid=(x(1:end-1)+x(2:end))/2;
kv=kfun(xmid);

supdiag=kv(1:end-1);
maindiag=-(kv(1:end-1)+kv(2:end));
subdiag=kv(2:end);

A=spdiags([subdiag maindiag supdiag],[-1 0 1],m,m);
format short g
full(A)
