%運用Crank-Nicolson差分格式求解二維擴散方程的初值問題；
clear?all
ti=0.00001;%時間步長
h=0.01;%空間步長
t=0.004;
k=t/ti;
m=1/h+1;
g=ti/h^2;?%網格比;
%初值條件及邊界條件
w=zeros（mm）;
u=zeros（mm）;
for?i=（0.4/h+1）:（0.6/h+1）
????for?j=（0.3/h+2）:（0.7/h）
????????u（ij）=1;
????end
end
for?i=（0.3/h+2）:（0.7/h）
????for?j=（0.4/h+1）:（0.6/h+1）
????????u（ij）=1;
????end
end
tic
for?l=1:k
u=Gauss_Seidel（ug）;%Gauss_Seidel迭代
end
toc
t=toc
mesh（0:h:10:h:1u）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????479??2010-11-18?22:11??程序\CN.m

?????文件????????613??2010-11-18?17:51??程序\DFF.m

?????文件????????455??2010-11-18?16:40??程序\draw\CN.m

?????文件????????619??2010-11-14?19:24??程序\draw\DFF.m

?????文件????????223??2010-11-18?16:40??程序\draw\draw.m

?????文件????????425??2010-11-14?16:31??程序\draw\Gauss_Seidel.m

?????文件????????614??2010-11-18?16:44??程序\draw\gdxs.m

?????文件???????1192??2010-11-18?10:04??程序\draw\ycjz.m

?????文件????????929??2010-11-14?19:53??程序\draw\ywjb.m

?????文件????????310??2008-04-29?17:21??程序\draw\zuigan.m

?????文件????????425??2010-11-14?16:31??程序\Gauss_Seidel.m

?????文件????????608??2010-11-18?17:52??程序\gdxs.m

?????文件???????1261??2010-11-18?17:53??程序\ycjz.m

?????文件????????970??2010-11-18?17:53??程序\ywjb.m

?????文件????????310??2008-04-29?17:21??程序\zuigan.m

?????文件?????864256??2011-01-26?10:15??第二次作業.doc

?????目錄??????????0??2010-11-18?22:13??程序\draw

?????目錄??????????0??2010-11-18?22:14??程序

-----------?---------??----------?-----??----

???????????????873689????????????????????18