function y = spl(x) global m n a y = log(1+exp(a*x))/a; return