%Solution (Matlab)
TA=100 ;
TB=200 ;
L=0.02 ;
q=1000000 ;
k=0.5 ;
A=[375,-125,0,0,0;-125,250,-125,0,0;0,-125,250,-125,0;0,0,-125,250,-125;0,0,0,-125,375]
B=[29000;4000;4000;4000;54000] ;
T=A\B
x=[0 0.002 0.006 0.01 0.014 0.018 0.02];
%T1=[100 150 218 254 258 230 200];
T1=T';
T1=[TA T' TB] ;
x1=0:0.0001:0.02;
y = ((TB-TA)/L + (L-x1)*(q/(2*k))).*x1+TA;
plot(x,T1,'o',x1,y)
xlabel ('x');
ylabel ('T (x)');
title ('Plaque infinie avec source de chaleur interne');
grid on;