function[y]=expfun(b,t_vec)

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