clcclear?all?close?all
ua?=?10;%賦初值
i?=?131;
j?=?51;
A_chafen（ij）?=?0;
Temp（ij）?=?0;
derta?=?0.01;%精度
womiga?=?1.0;%超松弛法系數(shù)
epise?=?8.85E-12;%介電常數(shù)
R1?=?10;%賦初值
R2?=?35;
R3?=?50;
C1?=?2*3.14*70*epise/log（R2/R1）;%計算C1
C2?=?2*3.14*60*epise/log（R3/R1）;%計算C2
for?j?=?12:35???????????????????	%B邊界
????A_chafen（1j）?=?10?-?（ua/log（R2/R1）*log（（j-1）/R1））;
end
for?j?=?12:50???????????????????????%C邊界
????A_chafen（131j）?=?10?-?（ua/log（R3/R1）*log（（j-1）/R1））;
end
A_chafen（1:1311:11）?=?10;??%內(nèi)邊界
Temp?=?A_chafen;
?for?n=1:50000
????for?i?=?2:71
????????for?j?=?12:35
????????????A_chafen（ij）?=?0.25*（A_chafen（i+1j）?+?A_chafen（i-1j）?+?（1+1/（2*（j-1）））*A_chafen（ij+1）?+?（1-1/（2*（j-1）））*A_chafen（