expandlog(log(a_*b_)):=log(a)+log(b);
expandlog(log(a_^b_)):= b*log(a);