Ta=20; hconv=0.01; L=10;
T1=40; TL=200;
h=1;
x=[0:h:L];

diag=2+h^2*hconv;
rhs=-h^2*hconv*Ta;

A=zeros(9,9);
n=length(x)-2;

for ii=1:n
	A(ii,ii)=-diag;

	if ii==1
		B(ii)=rhs-T1;
	elseif ii==n
		B(ii)=rhs-TL;
	else
		B(ii)=rhs;	
	end
end
for ii=1:n-1
	A(ii,ii+1)=1;
end
for ii=2:n
	A(ii,ii-1)=1;
end