multiplicative model modified

Previously, we rely on argument of seasonal in HoltWinters(). We found it is weak for some data (non-trend or non seasonal). 
We now use log-additive backtransform way to construct the multiplicative model. However, this way we can't get the most exactly right coefficient.