Etude de la répartition de la température

%*****************************************
% Etude de la répartition de température *
% dans un plaque *
%*****************************************
tic;
flops(0);
clear all; clc;clf;
eps=1e-4;k1=300;
% Données initiales
L=input ('Introduire la valeur de L '); %la longueur d'une plaque
w=input ('Introduire la valeur de w'); %la largeur d'une plaque
T0=20
Tc=80
x=0:0.2:2;
y=0:0.3:3;
for i=1:11;
for j=1:11
s=0;
for n=1:50
s=s+((1-(-1)^n)/(n*pi)).*(sinh(n*pi.*y(i)/2)./sinh(n*pi*3/2)).*sin(n*pi.*x(j)/2)
n=n+1
end
T(i,j)=2*s*(Tc-T0)+T0
end
end;
% Tracé des courbes de température en 3D
hold off;
if 11==11
figure(1);
T
i=1:11;
j=1:11;
grid on;
mesh(x,y,T);
Title('Evolution de la température dans un plaque');
xlabel('i');ylabel('j');zlabel('T (en °C)');
end
% Tracé des courbes isothermes en 2D
figure(2);
T
i=1:11;j=1:11;
grid on;
contour(i,j,T(i,j),11);
title('Lignes isothermes dans un plaque');
xlabel('i');ylabel('j');
t_mis=toc
Nb_opt=flops