# pow(a,b)=a^b; pow(1,n_):=1; pow(e,x_):=exp(x); pow(x_,1):=x; pow(-1,-0.5):= -i; pow(-1,0.5):= i; pow(-1,-1):= -1;