# convert to exp(x); #toexp(a_^b_):=If(isfree(a,x), e^(log(a)*b)); #toexp(a_*b_):=If(isfree(a,x), a*toexp(b), toexp(a)*toexp(b)); #toexp(a_+b_):=If(isfree(a,x), a+toexp(b), toexp(a)+toexp(b)); #toexp(a_-b_):=If(isfree(a,x), a-toexp(b), toexp(a)-toexp(b)); #toexp(1/b_):= 1/toexp(b); #toexp(-b_):= -toexp(b); #toexp(x_):=if(iscomplex(x),toexp(tocomplex(x))); #toexp(complex(a_,b_)):= abs(a,b)*exp(i*todegree(atan2(b,a))); toexp(cos(a_)+i*sin(a_)):=exp(a*i); toexp(a_^x_):=exp(log(a)*x); toexp(x_/W(x_)):=exp(W(x)); toexp(sqrt(x_)) := x^(1/2); toexp(cbrt(x_)) := x^(1/3); toexp(sinh(x_)) := exp(x)/2-exp(-x)/2; toexp(cosh(x_)) := exp(x)/2+exp(-x)/2; toexp(tanh(x_)) := (exp(x)-exp(-x))/(exp(x)+exp(-x)); toexp(coth(x_)) := (exp(x)+exp(-x))/(exp(x)-exp(-x)); toexp(csch(x_)) := 2/(exp(x)-exp(-x)); toexp(sech(x_)) := 2/(exp(x)+exp(-x)); toexp(sin(x_)) := i*exp(-i*x)/2-i*exp(i*x)/2; toexp(cos(x_)) := exp(i*x)/2+exp(-i*x)/2; toexp(tan(x_)) := (i*exp(-i*x)-i*exp(i*x))/(exp(i*x)+exp(-i*x)); toexp(cot(x_)) := (i*exp(i*x)+i*exp(-i*x))/(exp(i*x)-exp(-i*x)); toexp(csc(x_)) := 2i/(exp(-i*x)-exp(i*x)); toexp(sec(x_)) := 2/(exp(i*x)+exp(-i*x)); toexp(gamma(n_,x_)) := If(n>0, (n-1)*toexp(gamma(n-1,x))+x^(n-1)*exp(-x), If(n <0, (toexp(gamma(n+1,x))-x^n*exp(-x))/n) ); toexp(En(n_,x_)) := If(n>0, (-x*toexp(En(n-1,x))+exp(-x))/(n-1), If(n< 0, expand(-n*toexp(En(n+1,x))/x)+exp(-x)/x )); toexp(gamma(1/2,x_)) := -sqrt(pi)*erf(sqrt(x)); toexp(gamma(1/2,-x_)) := -sqrt(pi)*i*erfi(sqrt(x)); toexp(gamma(1/2,x_^2)) := -sqrt(pi)*erf(x); toexp(Ei(0.5,x_)) := sqrt(pi)*erfc(sqrt(x))/sqrt(x); toexp(Ei(-0.5,x_)) := sqrt(pi)*erf(sqrt(x))*x^(-1.5); toexp(mittag(1/2,x_)):= exp(x^2)*erfc(-x); toexp(mittag(2,x_)):= cosh(sqrt(x)); toexp(mittag(2,-x_)):= cos(sqrt(x)); toexp(mittag(4,x_)):=cos(x^(1/4))/2+cosh(x^(1/4))/2; toexp(gauss(x_)):= exp(-x^2/2)/sqrt(2pi); toexp(polar(a_,b_)):=a*exp(b*i);