here are some subjectively interesting functions that are not polynomials, accompanied with a range over which the function does not do anything weird (no singularities).
motivation was to test numerical integration. many quadrature schemes are exact for low-degree polynomials (future post), so we avoid polynomials.
we give the functions in Mathematica syntax for convenience of integration, followed by notes. some integrals Mathematica cannot do. all of these can be evaluated quickly to high precision by some way other than numerical integration, sometimes a special function, given in notes.
do some quadrature schemes work poorly if the first derivative is infinite (i.e., function is vertical) at an endpoint?
Integrate[Exp[-x^2/2]/Sqrt[2*Pi],{x,0,2}] : standard normal distribution, Erf, inflection at 1
Integrate[Sin[x],{x,0,Pi}] : integrating sine over a much larger range would even more strongly avoid polynomiality
Integrate[Cos[x],{x,0,Pi/2}]
Integrate[Cos[n*tau-x*Sin[tau]]/Pi,{tau,0,Pi}] : BesselJ[n,x], limits fixed, let n=1 x=10
Integrate[x^x,{x,0,1}] : sophomore's dream, limits fixed, defined as limit at 0, vertical at 0
Integrate[x^-x,{x,0,1}] : sophomore's dream, limits fixed, defined as limit at 0, vertical at 0
Integrate[Sqrt[x],{x,0,1}] : vertical at 0
Integrate[x^((1+Sqrt[5])/2-1),{x,0,1}] : vertical at 0, exponent is golden ratio
Integrate[Log[x],{x,1,E}]
Integrate[1/x,{x,1,E}]
Integrate[Exp[x],{x,0,1}]
Integrate[Exp[x]/x,{x,1/4,2}] : exponential integral
Integrate[Exp[-1/x],{x,0,1}] : exponential integral, defined as limit at 0, flat at zero
Integrate[1/Log[x],{x,0,1/2}] : logarithmic integral, defined as limit at 0
Integrate[Sqrt[1-x^2],{x,-1,1}] : semicircle area, vertical at endpoints
Integrate[Sqrt[1/(1-x^2)],{x,-1/2,1/2}] : 1/6 circle perimeter, we avoid integrating from -1 and 1 because function is undefined at endpoints
Integrate[1/(1+x^2),{x,0,1}] : atan(1) = pi/4
Integrate[Sqrt[1-k^2*Sin[x]^2], {x,0,Pi/2}] : complete elliptic integral of the second kind, arithmetic-geometric mean, limits fixed, ellipse perimeter, let k=255/256, we have used trigonometric form to avoid the function being undefined at endpoints
Integrate[1/Sqrt[1-k^2*Sin[x]^2], {x,0,Pi/2}] : complete elliptic integral of the first kind, limits fixed, let k=15/16, arithmetic-geometric mean, we have used trigonometric form to avoid the function being undefined at endpoints
Integrate[x^a*(1-x)^b, {x,0,1}] : beta function, limits fixed, let a=1/3 b=1/3, vertical at endpoints
future work: rescale, compute Taylor expansions around the midpoint of the interval.
No comments :
Post a Comment