Re: pure function to generate a list of integrals

*To*: mathgroup at smc.vnet.net*Subject*: [mg77089] Re: pure function to generate a list of integrals*From*: dimitris <dimmechan at yahoo.com>*Date*: Sun, 3 Jun 2007 06:10:06 -0400 (EDT)*References*: <f3r9er$19e$1@smc.vnet.net>

Hi. $VersionNumber 5.2 You could try In[4]:= f[g_, {x_, a_, b_}, opts___(*add options of NIntegrate ] := NInt[g, {x, a, b}, opts] (*the function you are interested in*) In[5]:= (f[2*x, {x, 0, #1}] & ) /@ Range[0.1, 0.9, 0.4] (*an application*) Out[5]= {NInt[2*x, {x, 0, 0.1}], NInt[2*x, {x, 0, 0.5}], NInt[2*x, {x, 0, 0.9}]} In[7]:= % /. NInt -> NIntegrate (*change to NIntegrate in order to get numerical results*) Out[7]= {0.01,0.25,0.81} For your second example In[11]:= (f[3*x, {x, 0, #1}] & ) /@ Range[0.1, 0.9, 0.4] Out[11]= {NInt[3 x,{x,0,0.1}],NInt[3 x,{x,0,0.5}],NInt[3 x,{x,0,0.9}]} In[12]:= % /. NInt -> NIntegrate Out[12]= {0.015,0.375,1.215} As another more advanced example consider In[38]:= (f[Sin[x], {x, 0, #1}, WorkingPrecision -> 30, MaxRecursion -> 12, PrecisionGoal -> 20] & ) /@ {10, 100, 1000} % /. NInt -> NIntegrate Out[38]= {NInt[Sin[x], {x, 0, 10}, WorkingPrecision -> 30, MaxRecursion -> 12, PrecisionGoal -> 20], NInt[Sin[x], {x, 0, 100}, WorkingPrecision -> 30, MaxRecursion -> 12, PrecisionGoal -> 20], NInt[Sin[x], {x, 0, 1000}, WorkingPrecision -> 30, MaxRecursion -> 12, PrecisionGoal -> 20]} Out[39]= {1.839071529076452452259,0.137681127712316065898,0.43762092370929700892} I hope you will find useful this workaround. Others may give you better solutions. But an important fact must be pointed out. "The problem is that NIntegrate is not listable..." you mention and at first glance you are right. In[40]:= Attributes[NIntegrate] Out[40]= {HoldAll, Protected} However NIntegrate behaves like "listable". In[46]:= g = HoldForm[NIntegrate[{x, 2*x, 3*x}, {x, 1, 3}]] Out[46]= NIntegrate[{x,2 x,3 x},{x,1,3}] In[47]:= ReleaseHold[g] Out[47]= {4.,8.,12.} It is the same situation with Plot for example In[48]:= Attributes[Plot] Out[48]= {HoldAll, Protected} But In[50]:= ReleaseHold[g /. NIntegrate -> Plot] (*plot to be displayed*) or Integrate In[51]:= Attributes[Integrate] Out[51]= {Protected, ReadProtected} but In[53]:= ReleaseHold[g /. NIntegrate -> Integrate] Out[53]= {4, 8, 12} There are now two attitudes you could face above situation for NIntegrate,Plot, Integrate. The first is the quick one: Ignore any incosistency and be happy about the fact these commands behave like these. The second is (again be happy(!) and) try to figure out what is going on. At first glance someone would say that the fact that Plot and NIntegrate behave like "listable" is astonising. However the situation is not like this. Functions that are Listable are those that have a single argument. But actually the functions can have more than one argument, and in this case equal length lists must be given for every argument. An example of this is: In[55]:= Times[{2,3,4},{x,o,e}] Out[55]= {2 x,3 o,4 e} Commands like NIntegrate and Plot DO not allow you to have a list for the second argument, the iterator. So the commands are NOT in fact Listable. Functions like Integrate and Plot are overloaded commands in the sense that they have one definition for the case of a non-List expression in the first slot and a different definition when the first slot contains a list. Indeed In[59]:= Information["Plot", LongForm -> False] "Plot[f, {x, xmin, xmax}] generates a plot of f as a function of x from xmin to xmax. Plot[{f1, f2, ... }, {x, xmin, xmax}] plots several functions fi."*Button[More..., ButtonData :> "Plot", Active -> True, ButtonStyle -> "RefGuideLink"] The two different definitions are apparent. I think something similar should have been added to the documentation of Integrate and NIntegrate (at least for version 5.2 this is not the case). Adopted this explanation one will continue his experiments with Mathematica being happy (always!) that he has understood something that look like an astonising incosistency. Since you are working with NIntegrate there are possibilities to deal with Principal Values Integrals that cannot be done analytically. In this case you should load the proper package In[63]:= Needs["NumericalMath`CauchyPrincipalValue`"] Say now that you have the list of the following functions with poles at the origin In[64]:= o = {1/x, 1/Sin[x], Tan[x]}; In view of previous explanation and considering that based on Help Brower the command CauchyPrincipalValue uses NIntegrate you could try In[65]:= CauchyPrincipalValue[o, {x, -1, {0}, 2}, AccuracyGoal -> 10] Surprisingly it doesn't work! Incosistency? Maybe yes, may be no! The truth is that this function has one definition; just for one argument in the first slot! In[67]:= CauchyPrincipalValue[Evaluate@o, {x, -1, {0}, 2}, AccuracyGoal -> 10] Not working either. In this case there is indeed no kind of "listability". Try instead In[71]:= (CauchyPrincipalValue[#1, {x, -1, {0}, 1/2}, AccuracyGoal -> 10] & ) / @ o Out[71]= {-0.6931471805602387, -0.7605693185090221, -0.4850422299423481} Best Regards Dimitris / Ruth Lazkoz Saez : > Hi everyone, > > I am trying to brush up a long code I have to make it more compliant > with the spirit of functional programming. I do not like to hear that > the kind of calculations I do should run faster in C, because I suspect > that if I managed to write good code in Mathematica it should be as > fast. So I have to go and improve my code chunk by chunk. > > My first problem is that I want to generate a pure function say f, > which, so that f[2, {0.1, 0.5, 0.9}] gives me the same output as > > {NIntegrate[2x, {x, 0, 0.1}], NIntegrate[2x, {x, 0, 0.5}], > NIntegrate[2x, {x, 0, 0.9}]} > > That is, I want to generate a list of numerical integrals of the same > function but making one of the integration limits change by taking > values from a list. > > I also want my function to admit two arguments (a number and a list) > because I want to be able to use the same definition to generate the > same output as for instance > > > {NIntegrate[3x, {x, 0, 0.1}], NIntegrate[3x, {x, 0, 0.5}], > NIntegrate[3x, {x, 0, 0.9}]} > > by evaluating f[3, {0.1, 0.5, 0.9}] this time. > > I tried for quite a while, but I failed. I suspect one of the problems > is NIntegrate is not listable. I could make some progress with Map but > I only what halfway and on top I was not satisfied with the syntax I > would have to use. > > Thanks in advance, > > Ruth Lazkoz