MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Gauss-Kronrod algorithm

  • To: mathgroup at smc.vnet.net
  • Subject: [mg77006] Re: Gauss-Kronrod algorithm
  • From: dimitris <dimmechan at yahoo.com>
  • Date: Thu, 31 May 2007 03:34:57 -0400 (EDT)
  • References: <f04nlc$8kr$1@smc.vnet.net><f3jhhf$dgk$1@smc.vnet.net>

Hello the following might be interesting for you

http://library.wolfram.com/infocenter/Conferences/6485/

A talk from Anton Antonov of WRI for NIntegrate (version 6)

Dimitris


 /  acmilan       :
> Dear Dimitris and Peter,
> Thanks so much for your help.
> Regards,
> Milan
> On May 4, 4:31 pm, dimitris <dimmec... at yahoo.com> wrote:
> > Hello again.
> >
> > You should have a look here.
> >
> > http://reference.wolfram.com/mathematica/tutorial/NIntegrateIntegrati...
> >
> > Dimitris
> >
> > =CF/=C7 dimitris =DD=E3=F1=E1=F8=E5:
> >
> >
> >
> > > Although this is from the documentation of the upcoming version of
> > > Mathematica
> > > I suggest you to have a look!
> >
> > >http://reference.wolfram.com/mathematica/ref/NIntegrate.html
> >
> > > Now let's see some things.
> >
> > > 1) Some useful links
> >
> > >http://en.wikipedia.org/wiki/Gaussian_quadrature
> > >http://numericalmethods.eng.usf.edu/topics/gauss_quadrature.html
> > > (amazing site!)
> > >http://www.sitmo.com/eqcat/13
> > >http://scholar.google.com/scholar?hl=en&lr=&q=%22Calculation+of+Gau=
> > ss=
> > > -Kronrod+Quadrature+Rules%22
> > > (the first article is of special notice)
> >
> > > 2) Some Mathematica code
> >
> > > In[64]:=
> > > oo[x_] = Exp[-5*x^4] + (1/2)*x;
> >
> > > In[71]:=
> > > Integrate[oo[x], {x, -3, 5}]
> > > {N[%], ({#1, NIntegrate[oo[x], {x, -3, 5}, Method -> #1]} & ) /@
> > > {Trapezoidal, DoubleExponential, GaussKronrod}}
> >
> > > Out[71]=
> > > 4 + (8*Gamma[5/4] - Gamma[1/4, 405] - Gamma[1/4, 3125])/(4*5^(1/4))
> >
> > > Out[72]=
> > > {5.212295737874928, {{Trapezoidal, 5.212295737874928},
> > > {DoubleExponential, 5.212295737874928},
> > >    {GaussKronrod, 5.212295737874985}}}
> >
> > > In[73]:=
> > > (ListPlot[Reap[NIntegrate[oo[x], {x, -3, 3}, Method -> #1,
> > > EvaluationMonitor :> Sow[x]]][[2,1]], PlotLabel -> #1] & ) /@
> > > {Trapezoidal, DoubleExponential, GaussKronrod};
> >
> > > I hope I help you a bit!
> >
> > > Dimitris
> >
> > > =CF/=C7 acmilan =DD=E3=F1=E1=F8=E5:
> > > > On Apr 24, 3:21 pm, dimitris <dimmec... at yahoo.com> wrote:
> > > > > Hello.
> >
> > > > > As a first hint I suggest you searching in the archives about posts
> > > > > of V. Antonov of WRI.
> >
> > > > > Dimitris
> >
> > > > > =CF/=C7 acmilan =DD=E3=F1=E1=F8=E5:
> >
> > > > > > On Apr 18, 5:16 pm, acmilan <fanvn_Acmi... at yahoo.com> wrote:
> > > > > > > Dear All,
> > > > > > > I am looking for the algorithm to calculate the general Gauss-Kro=
> > nr=
> > > od
> > > > > > > points and weights  in Mathematica( which is similar with the bui=
> > lt=
> > > -in
> > > > > > > function
> > > > > > > GaussianQuadratureWeights[n, a, b, prec]  for Gauss points and
> > > > > > > weights)
> > > > > > > Thank you very much for your help.
> > > > > > > Best regards,
> > > > > > > Milan
> >
> > > > > > Hi all,
> > > > > > Anybody helps me on this Mathematica Gauss-Kronrod points and weigh=
> > ts
> > > > > > code?
> > > > > > Thanks and regards,
> > > > > > Khoi Khoa- Hide quoted text -
> >
> > > > > - Show quoted text -
> >
> > > > Dear Dimitris,
> > > > I did try to search all the archives but could not find what I want.
> > > > Do you have any other suggestions?
> > > > Thanks and regards,
> > > > Milan- Hide quoted text -
> >
> > - Show quoted text -



  • Prev by Date: Re: Re: Ellipse equation simplification on Mathematica:
  • Next by Date: Re: LogLinearListPlot
  • Previous by thread: Re: Gauss-Kronrod algorithm
  • Next by thread: Re: Simplification