function out = modfunc1(x,a)
out = a(1) + a(2)*exp(a(3)*x);