function S = modelo(a,dados)
%
% modelo: f(t)_k = 1 - (a1 * exp(-t/a1) - a2 * exp (-t/a2))/(a1-a2)
%

t = dados(:,1); y_k = dados(:,2);

f_k = 1 - (a(1) * exp(-t/a(1)) - a(2) * exp (-t/a(2))) / (a(1)-a(2));

e = y_k - f_k;
S = e'*e;
