Ta=20; L=10; hconv=0.01;
T1=40; TL=200;
dx=1; coeff1=(2+hconv*dx^2); coeff2=hconv*dx^2*Ta;

%A= [coeff1 -1 0 0 0 0 0 0 0; -1 coeff1 -1 0 0 0 0 0 0; 0 -1 coeff1 -1 0 0 0 0 0; 0 0 -1 coeff1 -1 0 0 0 0; 0 0 0 -1 coeff1 -1 0 0 0; 0 0 0 0 -1 coeff1 -1 0 0; 0 0 0 0 0 -1 coeff1 -1 0; 0 0 0 0 0 0 -1 coeff1 -1; 0 0 0 0 0 0 0 -1 coeff1]

B = [coeff2 + T1; coeff2; coeff2; coeff2; coeff2; coeff2; coeff2; coeff2; coeff2 + TL]
x=[1:dx:9];
n=length(x);
for ii=1:n
	for jj=1:n
		if (ii==jj)
			A(ii,jj)=coeff1;
		elseif (jj==ii+1)
			A(ii,jj) = -1;
		elseif (jj==ii-1)
			A(ii,jj) = -1;
		else
			A(ii,jj) = 0;
		end
	end
end