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

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