%****二維橢圓方程的差分解法****%

%二維橢圓方程為：uxx+uyy=0

%?x方向的結點數?nx
%?x方向的步長?hx
%?求解區間x的左端：minx
%?求解區間x的右端：maxx
%?y方向的結點數?ny
%?y方向的步長?hy
%?求解區間y的左端：miny
%?求解區間y的右端：maxy
%?求解區間上的數值解?u

%%%處理求解區域上邊界的函數%%%
clear;
minx=100;
maxx=400;
miny=0;
maxy=300;
nx=8;
ny=8;
hx=（maxx-minx）/（nx-1）;
hy=（maxy-miny）/（ny-1）;
u0=zeros（nxny）;?%u0[nxny]二維數組存放邊界條件的值
x=minx:hx:maxx;?%差分網格x向坐標
y=miny:hy:maxy;?%差分網格y向坐標
Z=zeros（nxny）;?%Z[nxny]二維數組存放所有差分網格上的數值解（包括已知邊界網格點上的值）
[XY]=meshgrid（xy）;
for?i=1:nx
????%u0（1i）=UpBorder（maxyminx+（i-1）*hx）;?%處理求解區域的上邊界
????%u0（nyi）=DownBorder（minyminx+（i-1）*hx）;?%處理求解區域的下邊界
????u0（1i）=8+（i-1）*0.1;
????u0（nyi）=6.5+（i-1）*0.1;
end
for?j=1:ny
????%u0（j1）=LeftBorder（maxyminx+（i-1）*hx）;?%處理求解區域的左邊界
????%u0（jnx）=RightBorder（）;?%處理求解區域的右邊界
????u0（j

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

?????文件???????3780??2011-08-02?22:09??difference.m

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

?????????????????3780????????????????????1