S.put(k, exp(logA - N * logwn - B * exp(-M * logwn)))