Re: Problem with evaluating Integral

*To*: mathgroup at yoda.physics.unc.edu*Subject*: Re: Problem with evaluating Integral*From*: David Withoff <withoff>*Date*: Mon, 13 Apr 1992 15:44:27 -0500

John Yearwood writes: > Can anyone tell me why Mathematica v2.0 (for windows or Mac) does not > make any progress with the following: > > f[x_]:= Sin[x]+Exp[x] > j[x_]:=Evaluate[Normal[Series[f[x],{x,a,4}]]] > err[a_]:=Integrate[(f[x] - j[x])^2,{x,-6,2}] > Evaluate[err[a]] The expanded integrand here has 233 terms. For various reasons (cancelling non-elementary summands, etc.) the internal algorithm handles all of these terms at once, and ends up manipulating enormous intermediate expressions. In this example, however, it is perfectly safe to apply the idea that the integral of the sum is the sum of the integrals, and give the integration code a far simpler task: In[4]:= {time, result} = Timing[ Map[Integrate[#, {x,-6,2}] &, Expand[(f[x] - j[x])^2] ] ]; In[5]:= time Out[5]= 175.467 Second ... a little less than a second for each integral. Your timing may vary -- I did this on a SPARC. Dave Withoff withoff at wri.com