南开大学计算物理讲义9.ppt

b）隐格式 利用关系 可得三层隐式格式： 截断误差 : 绝对稳定 9.3.4 Model for general initial conditions Our analysis of the Euler method for solving the diffusion equation in section 9.3.3 assumed initial conditions of the form sin(kπx/Lx) sin(lπy/Ly) where k, l are integers and Lx, Ly are the dimensions of the domain. In addition to satisfying the boundary conditions, these initial conditions represent a set of orthogonal functions which may be used to construct any arbitrary initial conditions as a Fourier series. Now, since the diffusion equation is linear, and as our stability analysis of the previous section shows