programmation sur matlab

Bonjour cette équation suivante est celle de la chaleur en coordonnées cylindriques discrétisées:
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
Sa solution est de la forme: B.T_(n+1)=T_n+A.S_(n+1)
ou B,T et S sont des matrices.
Voici ce que j'ai pu faire et je ne peux plus avancer
%*************************************************
% Conduction de la chaleur dans plaque de Kapok *
%*************************************************
clear;clf;clc;hold off;
Tf=100;Ta=25;
alpha=8*10^(-6);dx=(5-1)/4;
bf=92;ba=13;T1=(-bf/(1-bf))*Tf;T5=ba*Ta/(1+ba);
p=0.13;dx=(20-1)/19;
%*******************************
% Construction de la matrice A *
%*******************************
A(1,1:2)=[1+2*p -1.25*p];A(1,3:20)=0;
A(20,1:18)=0;A(20,19:20)=[-0.75*p 1+2*p];
for i=2:19
for j=1:20
if i<j-1 & j>i+1
A(i,j)=0;
end
if i==j
A(i,j)=1+2*p;
A(i,j-1)=-0.75*p;
A(i,j+1)=-1.25*p;
end
end
end
S(1)=T1;S(20)=T5;S(2:19)=0;S=S';