T=293.15:5:323.15; k=[23.6,27.2,31,35.1,39.4,43.9,48.7]; fun = @(x,t) x(1).*exp(-x(2)./(t-x(3)));[x,r] = nlinfit(T,k,fun,[1,1,293]);plot(T,k,'*')hold onplot(T,fun(x,T))