clear
%% Constantes physiques
Ta=20; % oC % température à gauche
Tb=-20; % oC % température à droite
k=100; % W/m.oC % conductivité thermique
f_source=1e7; % W/m^3 % terme source constant
L=0.04; % m % longueur de la tige

%% Paramètres du modèle
nn=5; % nombre de noeuds
dx=???? %m % discrétisation en espace

%% Assemblage de la matrice A
A=[????]

%% Assemblage du vecteur b
% terme source
b= [ ???? ]

%% Résolution du système matriciel A{T}=b

T= ????

%% Impression des résultats à l'écran 
FigHandle = figure(1);
set(FigHandle, 'Position', [200, 200, 800, 600]);
plot(dx*[0:1:4],T,'b-')
hold on
plot(dx*[0:1:4],T,'ro','MarkerSize',18)
xlabel('x (m)','fontsize',16)
ylabel(['Température (' char(176) 'C)'],'fontsize',16)
title('Profil de température dans la tige','fontsize',16)
set(gca,'fontsize',14)