برنامه متلب حل معادله گرمای دوبعدی روی مربع واحد با روش Peaceman-rachford با جواب دقیق مشخص P-R Method for 2D heat equation on unit square % % u_t = u_{xx} + u_{yy} + f(x,t) % % Test problme: % % Exact solution: u(t,x,y) = exp(-pi^2*t) sin(pi*x+pi*y) % % Source term: f(t,x,y) = pi^2*exp(-pi^2