Programmation de l'équation de la chaleur

Bonjour,
je veux programmer l'équation de la chaleur discrétisée en coordonnées cylindriques mais j'ai des problémes aprés avoir déclaré mes données. je ne peut plus continuer.
Je mets ci joint cette équation et la programmation débutée.
T_i^(n-1)=-p.(1-1/(2.i)).T_(i-1)^n+(1+2.p).T_i^n-p.(1+1/(2.i)).T_(i+1)^n
B.T_(n+1)=T_n+A.S_(n+1)

clear;clf;clc;hold off;
k=80.2;ro=7870;Cp=447e3;TL=0;TR=200;
alpha=k/ro/Cp;dx=(50-1)/49;
%*******************************
% Construction de la matrice A *
%*******************************
A(1,1:2)=[-2 1];A(1,3:50)=0;
A(50,1:48)=0;A(50,49:50)=[1 -2];
for i=2:49
for j=1:50
if i<j-1 & j>i+1
A(i,j)=0;
end
if i==j
A(i,j)=-2;
A(i,j-1)=1;
A(i,j+1)=1;
end
end
end