clc;?clear?all;?close?all;

%2D?WAVE?EQUATION?utt?=?c^2（uxx+uyy）?
%with?initial?condition??u（xy0）?=?sin（p*pi*x）*sin（q*pi*y）?0%?and?boundary?conditions?u（0yt）?=?u（1yt）=?u（x0t）=?u（x1t）?=?0?t>0


c?=?1;??
dx?=?0.01;
dy?=?dx;
sigma?=?1/sqrt（2）;?gamma?=?1/sqrt（2）;?%Courant-Friedrich?Stability?Condition
dt?=?sigma*（dx/c）;
t?=?0:dt:1;?x?=?0:dx:1;?y?=?0:dy:1;?
u?=?zeros（length（x）length（y）length（t））;
p?=?2;?q?=?1;

u（::1）?=?transpose（sin（p.*pi.*x））*sin（q.*pi.*y）;?%u（xy0）?=?sin（p*pi*x）*sin（q*pi*y）

%u（xydt）
for?i=2:length（x）-1?
????for?j=2:length（y）-1
?????u（ij2）=?（sigma^2）*（u（i+1j1）-2*u（ij1）+u（i-1j1））...
?????????+（gamma^2）*（u（ij+11）-2*u（ij1）+u（ij-11））+2*u（ij1）?-?u（ij1）;?
????end
end

for?n=2:length（t）-1
????for?i=2:length（x

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2123??2017-03-29?19:36??Wave_Equation_2D.m
?????文件????????1314??2017-03-29?19:36??license.txt