unsimplified result from integrate
- Subject: unsimplified result from integrate
- From: Barton Willis
- Date: Sun, 13 May 2007 23:04:32 -0500
Consider (a test from rtest5):
(%i1) e : integrate(%e^x*erf(x),x);
(%o1) %e^x*erf(x)-%e^(1/4)*erf(x-1/2)
Yeechs--look at ((mtimes simp) ((mtimes simp) ...)) -- that's not
simplified:
(agree / disagree ?)
(%i2) ?print(e);
((MPLUS SIMP)
((MTIMES SIMP) ((MTIMES SIMP) -1 ((MEXPT SIMP) $%E ((RAT SIMP) 1 4)))
((%ERF SIMP) ((MPLUS SIMP) ((RAT SIMP) -1 2) $X)))
((MTIMES SIMP) ((MEXPT SIMP) $%E $X) ((%ERF SIMP) $X)))
The function resimplify crunches e:
(%i3) :lisp(resimplify $e));
((MPLUS SIMP)
((MTIMES SIMP) -1 ((MEXPT SIMP) $%E ((RAT SIMP) 1 4))
((%ERF SIMP) ((MPLUS SIMP) ((RAT SIMP) -1 2) $X)))
((MTIMES SIMP) ((MEXPT SIMP) $%E $X) ((%ERF SIMP) $X)))
Barton