function PDE1Dd_MOL
% 用MOL法求解一維動態偏微分方程：
% 先將PDE方程用有限差分法變為ODE方程，然后用龍格-庫塔法積分求解
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2002/05/21 $

global alpha dx Ta Tb
alpha = 2.0;
dx = 2.0;
Ta = 100;
Tb = 0;
tspan = [0 8];
T0 = [0 0 0 0];
[t,T] = ode45(@ODEs,tspan,T0)

% --------------------------------------------------------------------------
function dTdt = ODEs(t,T)
global alpha dx Ta Tb
dTdt1 = Ta - 2*T(1) + T(2);
dTdt2 = T(1) - 2*T(2) + T(3);
dTdt3 = T(2) - 2*T(3) + T(4);
dTdt4 = T(3) - 2*T(4) + Tb;
dTdt = alpha/dx^2 * [dTdt1; dTdt2; dTdt3; dTdt4];