معادله پخش به صورت زیر می باشد:

روش ترکیبی b :

ریخمایر و مورتن الگوریتم کلی زیر را برای روش ضمنی با سه سطح زمانی ارائه دادند:

 کد مورد نظر (tarkibi_B.m) ،  معادله گرما را در حالت دو بعدی و به روش ترکیبی B حل می نماید.

پارامترهایی که پس از اجرای برنامه از کاربر خواسته می شوند:

طول،عرض، تعداد تقسیمات در امتداد x و y ، تعداد گام زمانی، ضریب پخش، شرایط مرزی، شرط اولیه و…

……………..

شرایط مرزی برای  مرزها به صورت زیر می باشند:

مرز بالا،پایین، چپ: نوع دیریکله

مرز چپ:نیومن

پس از اجرای کد ، توزیع دما به صورت سه بعدی به نمایش در می آید.