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

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

%?x方向的結(jié)點(diǎn)數(shù)?nx
%?x方向的步長(zhǎng)?hx
%?求解區(qū)間x的左端：minx
%?求解區(qū)間x的右端：maxx
%?y方向的結(jié)點(diǎn)數(shù)?ny
%?y方向的步長(zhǎng)?hy
%?求解區(qū)間y的左端：miny
%?求解區(qū)間y的右端：maxy
%?求解區(qū)間上的數(shù)值解?u

%%%處理求解區(qū)域上邊界的函數(shù)%%%
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]二維數(shù)組存放邊界條件的值
x=minx:hx:maxx;?%差分網(wǎng)格x向坐標(biāo)
y=miny:hy:maxy;?%差分網(wǎng)格y向坐標(biāo)
Z=zeros（nxny）;?%Z[nxny]二維數(shù)組存放所有差分網(wǎng)格上的數(shù)值解（包括已知邊界網(wǎng)格點(diǎn)上的值）
[XY]=meshgrid（xy）;
for?i=1:nx
????%u0（1i）=UpBorder（maxyminx+（i-1）*hx）;?%處理求解區(qū)域的上邊界
????%u0（nyi）=DownBorder（minyminx+（i-1）*hx）;?%處理求解區(qū)域的下邊界
????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）;?%處理求解區(qū)域的左邊界
????%u0（jnx）=RightBorder（）;?%處理求解區(qū)域的右邊界
????u0（j

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

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

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

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