function[y]=expfunc(b,t_vec)

y = b(1)-exp(-b(2)*t_vec);