Hogy kell ezt megoldani Matlabban?

Szeretnék egy olyan "scriptet" írni, ami Newton-féle hűlési törvényt felhasználva számolja ki a "k" értékét, és ebből a "k" értékből adott idő múlva kiszámolja a várható hőmérsékletet.

Newton-féle törvény: To+(To-Tk)*exp(-k*t)

Ennyit hoztam össze eddig:

function [k,T] = coollaw(a,b,d,e,f)

k=real(-(log((e-a)/(a-b)))/d);

T=a+(a-b)*exp(-k*f);

k értékét jól számolja ki, de utána számolt hőmérséklet nem az, amit az ellenőrzés során kapok....

Mi itt a hiba?