معادله قانون سرمایش نیوتن به صورت زیر می باشد:

T'(t) + k*T(t) = k*Ts

فایل cooling.m ، معادله را به روش کرنک نیکلسون حل می کند و در پایان ، تغییرات دما بر حسب زمان به نمایش در می آید
پارامترهای ورودی:

tn = زمان نهایی
k = ضریب هدایت جسم
Ts= دمای محیط
T1=دمای اولیه جسم