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

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;

endfunction
