On Mar 30, 9:48 am, Minh Nguyen <[email protected]> wrote: > Hi Jared, > > On Wed, Mar 31, 2010 at 12:31 AM, Jared Schlieper > > <[email protected]> wrote: > > <SNIP> > > > Am I missing something? > > > Mathematica gives 3954.63 using the NIntegrate command. > > Is the following what you want? > > [mv...@sage ~]$ sage > ---------------------------------------------------------------------- > | Sage Version 4.3.3, Release Date: 2010-02-21 | > | Type notebook() for the GUI, and license() for information. | > ---------------------------------------------------------------------- > sage: x, y = var("x, y") > sage: density = e^(sqrt(x)) + e^(y^2/2) > sage: g1 = e^(2*x) - 1 > sage: g2 = 5 - 5*(x-1)^2 > sage: a = (g2 - g1).find_root(-0.1, 0.2, x) > sage: b = (g2 - g1).find_root(.7, 1.0, x) > sage: f = integral(density, y, g1, g2) > sage: f.nintegral(x, a, b) > (3954.6341726670939, 2.7847910517020398e-06, 105, 0) >
Thanks. That works for what I need. Jared > -- > Regards > Minh Van Nguyen -- To post to this group, send email to [email protected] To unsubscribe from this group, send email to [email protected] For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org To unsubscribe, reply using "remove me" as the subject.
