On Nov 15, 2007 7:49 AM, William Stein <[EMAIL PROTECTED]> wrote: > > f := x -> integrate(exp(t^2), t=0..x); (Maple) > > evalf(f(2)); ... > > You can define that function and compute the integral in Sage as > follows: > > sage: assume(x > 0) > sage: f(x,t) = integrate(exp(t^2), t, 0, x) > sage: f > (x, t) |--> -sqrt(pi)*I*erf(I*x)/2 > sage: a = f(x,0) > sage: a > -sqrt(pi)*I*erf(I*x)/2
By the way, one can compute the numerical value at 2 in Sage of this function as follows: sage: val, err_bound = numerical_integral(lambda t: exp(t^2), 0, 2) sage: val 16.452627765507231 sage: err_bound 1.826608616085327e-13 William --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~----------~----~----~----~------~----~------~--~---
