gistfile1.txt

(%i1) integrate(x^2 * %e ^ (-x ^ 2), x);
                                                    2
                                                 - x
                          sqrt(%pi) erf(x)   x %e
(%o1)                     ---------------- - --------
                                 4              2
(%i2) plot2d(%o1, [x, -1.5, 1.5]);