function res=a1_function(a,b)
    res = -b + a*exp(-a/10);
return