يك صفحه به طول L متر و عرض H متر مفروض است. اگر صفحه را به قسمت هاي كوچك تقسيم نماييم، مي توان دماي هر نقطه در حالت پايا (steady state) را متوسط دماي چهار قسمت اطراف آن دانست. اگر دماي چهار طرف صفحه اصلي مشخص باشد، برنامه اي بنويسيد كه با فرض گفته شده، دماي نقاط مياني را محاسبه نموده و پروفايل دما داخل صفحه را رسم نمايد.

برنامه را به گونه اي بنويسيد كه بتوان تعداد قسمت ها را در هر جهت تغيير داد و اثر آن را در پروفايل دما خروجي مشاهده نمود.

function NumericalSolution
t1=200;
t2=300;
t3=400;
t4=500;
x=30; % length grid
y=30; %width grid

for i=1:x
s(y,i)=t3;
end
s(:,1)=t4;
s(:,x)=t2;
k=s.*0;
err=s.*5; %error value
while max(max(err))>0
k=s;
for i=2:y-1
for j=2:x-1
s(i,j)=(s(i+1,j)+s(i-1,j)+s(i,j+1)+s(i,j-1))/4;
end
end
err=abs(k-s);
end
mesh(s)
xlabel('length')
ylabel('width')
zlabel('temp')