hconv=0.01;T1=40; T11=200; Ta=20;
dx=1;
diag=-2-hconv*dx^2;
rhs=-hconv*dx^2*Ta;
A=zeros(9,9);

for i=1:9
	A(i,i)=diag;
end
for i=1:9-1
	A(i,i+1)=1;
end
for i=2:9
	A(i,i-1)=1;
end
for i=1:9
	if i==1
		B(i)=rhs-T1;
	elseif i==9
		B(i)=rhs-T11;
	else
		B(i)=rhs;
	end
end