我們求解的問題是
上述問題的精確解為 u=sin(3*pi*x+pi/4)*sin(2*pi*y+pi/4);
接下來用Matlab工具箱求解
首先打開matlab，命令窗口處輸入：pdetool
將看到如下界面：
顯示柵格：Options->Grid
選擇求解區域：Draw (這個地方我們的求解區域為[0 1]X[0 1])->Rectange/square 在面板上隨意拖動，雙擊出現的區域，輸入如下數值
進入邊界條件設置：
上邊界與下邊界 sin(3*pi*x+pi/4)*sin(pi/4)
(注意：此處不能為'x',后同)
左邊界 3*pi*sin(2*pi*y+pi/4)*cos(pi/4)
右邊界 -3*pi*sin(2*pi*y+pi/4)*cos(pi/4)
方程設置：13*pi*pi.*sin(3*pi*x+pi/4).*sin(2*pi*y+pi/4)
網格
求解