Skip to main content

calculus and analysis - Integrate gives wrong results


Integrate[a/(Sin[t]^2 + a^2), {t, 0, 2 Pi}]


∫2π0aa2+sin2(t)dt


gives 0


This cannot be true. What is going on?


If I insert a number into a, it gives a reasonable result:


NIntegrate[2/(Sin[t]^2 + 4), {t, 0, 2 Pi}]

give 2.80993



Answer



$Version


(*
Out[228]= "8.0 for Microsoft Windows (64-bit) (October 7, 2011)"
*)

The result provided by Mathematica is correct:


Integrate[a/(a^2 + Sin[t]^2), {t, 0, 2 π}]

(*
Out[213]= (2 π)/(Sqrt[1 + 1/a^2] a)

*)

Now the same procedure as "always" which "explains" the zero result. The indefinite integral is


Integrate[a/(a^2 + Sin[t]^2), t]

(*
Out[214]= ArcTan[(Sqrt[1 + a^2] Tan[t])/a]/Sqrt[1 + a^2]
*)

This result is not continuous for real a. In this graph we can see the jumps as well as the necessary additional terms to restore continuity.



The jump must be determined by taking (directed) limits. Assuming first a>0 gives


Simplify[Limit[ArcTan[(Sqrt[1 + a^2] Tan[t])/a]/Sqrt[1 + a^2], 
t -> π/2, Direction -> +1], a > 0] -
Simplify[Limit[ArcTan[(Sqrt[1 + a^2] Tan[t])/a]/Sqrt[1 + a^2],
t -> π/2, Direction -> -1], a > 0]

(*
Out[408]= π/Sqrt[1 + a^2]
*)


and therefore


With[{a = 1/2}, 
Plot[{ArcTan[(Sqrt[1 + a^2] Tan[t])/a]/Sqrt[1 + a^2],
ArcTan[(Sqrt[1 + a^2] Tan[t])/a]/Sqrt[1 + a^2] + π/Sqrt[1 + a^2],
ArcTan[(Sqrt[1 + a^2] Tan[t])/a]/Sqrt[1 + a^2] + (2 π)/Sqrt[
1 + a^2]}, {t, 0, 2 π}]]

150601_plot_int.jpg


The blue (lower) curve is the original antiderivative, the red (middle) curve has π/Sqrt[1+a^2] added, and is the continiuous continuation between π/2 and 3π/2, finally the brown (upper) curve does it for the rest adding again the same amount π/Sqrt[1+a^2].


Hence taking the difference of the antiderivative at the endpoints according to the fundamental theorem of calculus would lead to zero on the blue curve, and is only correct for the continuous version leading to the correct result given in the beginning.



EDIT #1


The case a<0 need not be discussed separately because of the (anti)symmetry of the integral.


Note that because the antiderivative vanishes at both ends the integral is equal to the "total" jump, which here is twice the amount of the jump at π/2.


This and other examples point to the (tentative) rule: if the result of Integrate[] is zero despite the fact that the integrand is positive check for jumps and calculate them using the directed Limits. The integral will then be the sum of the jumps over the whole interval.


Comments

Popular posts from this blog

functions - Get leading series expansion term?

Given a function f[x] , I would like to have a function leadingSeries that returns just the leading term in the series around x=0 . For example: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x)] x and leadingSeries[(1/x + 2 + (1 - 1/x^3)/4)/(4 + x)] -(1/(16 x^3)) Is there such a function in Mathematica? Or maybe one can implement it efficiently? EDIT I finally went with the following implementation, based on Carl Woll 's answer: lds[ex_,x_]:=( (ex/.x->(x+O[x]^2))/.SeriesData[U_,Z_,L_List,Mi_,Ma_,De_]:>SeriesData[U,Z,{L[[1]]},Mi,Mi+1,De]//Quiet//Normal) The advantage is, that this one also properly works with functions whose leading term is a constant: lds[Exp[x],x] 1 Answer Update 1 Updated to eliminate SeriesData and to not return additional terms Perhaps you could use: leadingSeries[expr_, x_] := Normal[expr /. x->(x+O[x]^2) /. a_List :> Take[a, 1]] Then for your examples: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x), x] leadingSeries[Exp[x], x] leadingSeries[(1/x + 2 + (1 - 1/x...

mathematical optimization - Minimizing using indices, error: Part::pkspec1: The expression cannot be used as a part specification

I want to use Minimize where the variables to minimize are indices pointing into an array. Here a MWE that hopefully shows what my problem is. vars = u@# & /@ Range[3]; cons = Flatten@ { Table[(u[j] != #) & /@ vars[[j + 1 ;; -1]], {j, 1, 3 - 1}], 1 vec1 = {1, 2, 3}; vec2 = {1, 2, 3}; Minimize[{Total@((vec1[[#]] - vec2[[u[#]]])^2 & /@ Range[1, 3]), cons}, vars, Integers] The error I get: Part::pkspec1: The expression u[1] cannot be used as a part specification. >> Answer Ok, it seems that one can get around Mathematica trying to evaluate vec2[[u[1]]] too early by using the function Indexed[vec2,u[1]] . The working MWE would then look like the following: vars = u@# & /@ Range[3]; cons = Flatten@{ Table[(u[j] != #) & /@ vars[[j + 1 ;; -1]], {j, 1, 3 - 1}], 1 vec1 = {1, 2, 3}; vec2 = {1, 2, 3}; NMinimize[ {Total@((vec1[[#]] - Indexed[vec2, u[#]])^2 & /@ R...

plotting - Plot 4D data with color as 4th dimension

I have a list of 4D data (x position, y position, amplitude, wavelength). I want to plot x, y, and amplitude on a 3D plot and have the color of the points correspond to the wavelength. I have seen many examples using functions to define color but my wavelength cannot be expressed by an analytic function. Is there a simple way to do this? Answer Here a another possible way to visualize 4D data: data = Flatten[Table[{x, y, x^2 + y^2, Sin[x - y]}, {x, -Pi, Pi,Pi/10}, {y,-Pi,Pi, Pi/10}], 1]; You can use the function Point along with VertexColors . Now the points are places using the first three elements and the color is determined by the fourth. In this case I used Hue, but you can use whatever you prefer. Graphics3D[ Point[data[[All, 1 ;; 3]], VertexColors -> Hue /@ data[[All, 4]]], Axes -> True, BoxRatios -> {1, 1, 1/GoldenRatio}]