- fixed bug in normal(): normal(x^a) became (x^(-a))^(-1)