Phil Lucht Math & Physics Archive
Home / Math and Physics Files / Math / Integrals series sums+ GR / Study of 1 over sqrt(R) integral

more study of the Goldstine p 77 integral

DOCX · 109.7 KB
Open DOCX file

Phil's working notes dated 11.3.16 that collect and extend earlier derivations for the integral of 1/sqrt(R), R = a+bx+cx^2. They derive the arcsinh form from the log form, continue it to arcsin and arccos forms by rotating c to negative values, and list conditions on a, b, c. The notes then show a sign error in Goldstein's (3-45) that does not change the orbit equation (3-46), and list other forms from Euler substitutions.

AI-written summary; may contain errors.

Extracted text (machine-read; may contain errors)
More study of the 1/sqrt(R) integral PhL 11.3.16 Here I summarize work done in two other docs, and derive new forms as well. It was a royal pain, but I think everything is down pat. Ultimately you just differentiate to test a result, perhaps with Maple help. But I wanted to derive results, not test results handed to me. 1. Derive the arcsinh formula 1 2. Obtain the arcsin formula 2 3. Obtain the arccos formula 5 4. Summary of formulas derived above 5 5. Goldstein page 77 Application 5 6. Other forms of the integral 7 1. Derive the arcsinh formula Here is a statement of the arcsinh formula of interest from GR7 : I have derived the first form involving the log. I understand that the denominator in the log result is not significant and the basic fact is this: I = !Syntax Error, Idx/ = (1/) ln [ 2cx + b + 2] + constant where the constant can depend on a,b,c but not on x. Question: How can one use this last result to derive the sh-1 version of the result stated above, namely, I = !Syntax Error, Idx/ = (1/) sh-1 [ (2cx +b)/] + constant It must be possible to show that sh-1 [ (2cx +b)/] = ln [ 2cx + b + 2] + constant that does not depend on x In this case, differentiating both functions should give 1/. The sh-1 form also appears in Schaum p 72. One thing I know is, from Schaum p 29, sh-1y = ln(y + ) // valid for all real y Using this, I can obtain a "result to show" that involves only logs. y = (2cx +b)/ y2+1 = (2cx +b)2/(4ac-b2) + 1 = [ (2cx +b)2 + (4ac-b2) ] / (4ac-b2) = [ 4c2x2+ b2 + 4cbx + 4ac-b2 ] / (4ac-b2) = [ 4c2x2 + 4cbx + 4ac] / (4ac-b2) = 4c[ cx2 + bx + a] / (4ac-b2) = 4cR/(4ac-b2) = 2 / . Then sh-1 [ (2cx +b)/] = sh-1y = ln(y + ) = ln [ (2cx +b)/ + 2/] = ln [ 2cx + b + 2] - ln() and there you have it, case closed. Thus I have proven both these results, I = !Syntax Error, Idx/ = (1/) ln [ 2cx + b + 2] + constant I = !Syntax Error, Idx/ = (1/) sh-1 [ (2cx +b)/] + constant This last form assumes that 4ac - b2 > 0 and c > 0. Let's summarize this as I = (1/) sh-1 [ (2cx +b)/] 4ac - b2 > 0 and c > 0 [ a > 0 is forced ] This situation implies that 4ac > b2 and so a > b2/(4c) so it is implied that a > 0. 2. Obtain the arcsin formula Let us now assume that a > 0 before continuing. Then 4ac-b2 is going to go negative for c < 0. Consider this picture At θ = 0, the vector 4ac-b2 lies on the real axis and points to the right assuming 4ac-b2 > 0. But as we swing the vector c over to θ = +π, we end up with 4ac pointing to the left and 4ac-b2 also points even more to the left, so here is the status of things angle(c) = 0 angle(4ac-b2) = 0 // before swing angle(c) = π angle(4ac-b2) = π // after swing We end up with 4ac - b2 = -4a|c| - b2 . Since we have assumed a > 0, we have 4ac - b2 < 0 after doing this swing operation, so this quantity changed sign. After the swing we have c = |c| eiπ = (-c)eiπ (4ac-b2) = |4ac-b2| eiπ = |b2-4ac| eiπ = (b2-4ac)eiπ Therefore (note that the phases are determined here, you don't pick separate ± i for each one) = i = i We can then modify our integral I above using these rules to get I = (1/) sh-1 [ (2cx +b)/] = (1/i ) sh-1 [ (2cx +b)/( i)] = (-i)(1/ ) sh-1 [ -i(2cx +b)/] = (-i)2(1/ ) sin-1 [ (2cx +b)/] = - (1/ ) sin-1 [ (2cx +b)/] Therefore we find I = - (1/) sin-1 [ (2cx +b)/] 4ac - b2 < 0 and c < 0 [ a > 0 assumed ] and this is our first statement of the sin-1 evaluation of the I integral. Now recall from GR7 Recall our two results from above: I = (1/) sh-1 [ (2cx +b)/] 4ac - b2 > 0 and c > 0 [ a > 0 is forced ] I = - (1/) sin-1 [ (2cx +b)/] 4ac - b2 < 0 and c < 0 [ a > 0 assumed ] This agrees exactly with the lines 2 and 4 of the quotes above. I have already derived the other two lines in a separate doc, so all four forms are now derived. Now in the first case a > 0 is forced, while in the second it was only assumed. In the second case we need 4ac < b2 -4a|c| < b2 4a|c| > -b2 a > -b2/(4|c|) This in fact allows a range of negative a without impacting the sin-1 formula shown above. In the first case we have 4ac > b2 a > b2/(4c) So let's restate the above two evaluations this way I = (1/) sh-1 [ (2cx +b)/] 4ac - b2 > 0 and c > 0 [ a > b2/(4c) ] (1) I = - (1/) sin-1 [ (2cx +b)/] 4ac - b2 < 0 and c < 0 [ a > -b2/(4|c|) ] (2) 3. Obtain the arccos formula It is a fact that sin-1(x) + cos-1(x) = π/2 where x is real, this is a very obvious fact when stated in English text. There is no corresponding equation for hyperbolic functions!!!! So we can shift the constant and write the second result above as I = + (1/) cos-1 [ (2cx +b)/] 4ac - b2 < 0 and c < 0 [ a > -b2/(4|c|) ] (3) There is no simple cosh-1 formula. However, I did derive a certain cosh-1 formula with a completely different argument than that shown above. I will gather up my obscure formulas later on. There is another version of the cos-1 which follows from (2), I = (1/) sin-1 [ -(2cx +b)/] 4ac - b2 < 0 and c < 0 [ a > -b2/(4|c|) ] (4) Then I = - (1/) cos-1 [ -(2cx +b)/] 4ac - b2 < 0 and c < 0 [ a > -b2/(4|c|) ] (5) You get this same result going directly from the previous cos result, using Schaum p 18. 4. Summary of formulas derived above I = (1/) sh-1 [ (2cx +b)/] 4ac - b2 > 0 and c > 0 [ a > b2/(4c) ] (1) I = - (1/) sin-1 [ (2cx +b)/] 4ac - b2 < 0 and c < 0 [ a > -b2/(4|c|) ] (2) I = + (1/) cos-1 [ (2cx +b)/] 4ac - b2 < 0 and c < 0 [ a > -b2/(4|c|) ] (3) I = (1/) sin-1 [ -(2cx +b)/] 4ac - b2 < 0 and c < 0 [ a > -b2/(4|c|) ] (4) I = - (1/) cos-1 [ -(2cx +b)/] 4ac - b2 < 0 and c < 0 [ a > -b2/(4|c|) ] (5) 5. Goldstein page 77 Application I claim that Goldstein page (3-45) is item (5) above but with an overall sign error. This was my conclusion a long time ago, though I did not document it at that time (it was marked in pencil). Let us now see if this sign error affects anything of significance. He has a = 2mE/l2 and b = 2mk/l2 and c = -1. Going back to his (3.44) I would say θ = θ' - I = θ' - { - (1/) cos-1 [ -(2cx +b)/] + K} = θ' + (1/) cos-1 [ -(2cx +b)/] + K Now consider b2-4ac = [2mk/l2]2 + 4[ 2mE/l2] = [2mk/l2]2 { 1 +4[ 2mE/l2]/ [2mk/l2]2 } = [2mk/l2]2 { 1 +4[ 2mE/l2]/ [4m2k2/l4]} = [2mk/l2]2 { 1 +[ 2E]/ [mk2/l2]} = [2mk/l2]2 { 1 +2El2/mk2} = q // agrees with page 77 α and also - 2cx - b = +2u - [2mk/l2] = [2mk/l2] { 2u/ [2mk/l2] - 1 } = [2mk/l2] { ul2/ mk - 1 } The above is then θ = θ'+ K + cos-1 [ [2mk/l2] { ul2/ mk - 1 } / [2mk/l2] ] = θ' + K + cos-1 [ (ul2/ mk - 1)/ ] = θ' + cos-1 [ (ul2/ mk - 1)/ ] // set K = 0 So I claim he has a sign error as well in equation β. Continuing. θ-θ' = cos-1 [ (ul2/ mk - 1)/ ] cos(θ-θ') = (ul2/ mk - 1)/ Notice that this last equation is exactly the same with or without the sign error, because cos is even !! Continue along then cos(θ-θ') = (ul2/ mk - 1) u(l2/mk) = 1 + cos(θ-θ') 1/r = (mk/l2) [ 1 + cos(θ-θ') ] // agrees with (3-46) So Goldstein has a sign wrong in (3-45) and in β, but this does not affect (3-46) !! 6. Other forms of the integral These other proven-valid forms I found by doing Euler substitutions, I = GR7(x) = + (1/) ln [ +2cx + b + 2] I = PL(x) = - (1/) ln [ - 2cx - b + 2] I = (2/) tanh-1 [ ] I = (2/) tanh-1 [ ] I = (2/) sech-1 ( ) = (2/) sech-1 ( ) I = (2/)cosh-1( ) In general these are unpleasant formulas because they involve and worse.