A curious integral through an algebraic equation

Calculus Level 5

If 0 1 1 + 1 + y x y d x = π \displaystyle \int_0^1 \sqrt{-1 + \sqrt{\frac {1+y}x-y}}\ dx = \pi , find y 2 78 y y^2 - 78y .


The answer is 15.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

2 solutions

Mark Hennings
Jan 23, 2020

Let us assume that y > 0 y > 0 .

The substitutions u 2 = 1 + y x y u^2 = \tfrac{1+y}{x} - y and v = u 1 v = u - 1 gives A = 0 1 1 + 1 + y x y d x = 2 ( 1 + y ) 1 u u 1 ( u 2 + y ) 2 d u = 2 ( 1 + y ) 0 ( v + 1 ) v ( ( v + 1 ) 2 + y ) 2 d v A \; = \; \int_0^1 \sqrt{-1 + \sqrt{\frac{1+y}{x} - y}}\,dx \; = \; 2(1+y)\int_1^\infty \frac{u\sqrt{u-1}}{(u^2+y)^2}\,du \; = \; 2(1+y)\int_0^\infty \frac{(v+1)\sqrt{v}}{\big((v+1)^2+y\big)^2}\,dv Defining v \sqrt{v} on the cut plane v > 0 |v| > 0 , 0 < A r g v < 2 π 0 < \mathrm{Arg}\,v < 2\pi , and integrating around a suitable keyhole contour gives A = 2 π i ( 1 + y ) ( R e s v = 1 + i y ( v + 1 ) v ( ( v + 1 ) 2 + y ) 2 + R e s v = 1 i y ( v + 1 ) v ( ( v + 1 ) 2 + y ) 2 ) A \;= \; 2\pi i(1 + y)\Big(\mathrm{Res}_{v=-1+i\sqrt{y}}\frac{(v+1)\sqrt{v}}{\big((v+1)^2 + y\big)^2} + \mathrm{Res}_{v=-1-i\sqrt{y}}\frac{(v+1)\sqrt{v}}{\big((v+1)^2 + y\big)^2} \Big) Now the function v f y ( v ) = ( v + 1 ) v ( ( v + 1 ) 2 + y ) 2 v \; \mapsto \; f_y(v) \; = \; \frac{(v+1)\sqrt{v}}{\big((v+1)^2 + y\big)^2} has a double pole at each of v = 1 ± i y v = -1 \pm i\sqrt{y} , and hence R e s y = 1 + i y ( v + 1 ) v ( ( v + 1 ) 2 + y ) 2 = d d v ( v + 1 ) v ( v + 1 + i y ) 2 v = 1 + i y = 1 + i y ( y i ) 8 y ( 1 + y ) \mathrm{Res}_{y = -1+i\sqrt{y}}\frac{(v+1)\sqrt{v}}{\big((v+1)^2 + y\big)^2} \; = \; \frac{d}{dv} \frac{(v+1)\sqrt{v}}{(v + 1 + i\sqrt{y})^2} \Big|_{v=-1+i\sqrt{y}} \; = \; -\frac{\sqrt{-1+i\sqrt{y}}(\sqrt{y}-i)}{8\sqrt{y}(1+y)} and, similarly, R e s y = 1 i y ( v + 1 ) v ( ( v + 1 ) 2 + y ) 2 = 1 i y ( y + i ) 8 y ( 1 + y ) \mathrm{Res}_{y = -1-i\sqrt{y}}\frac{(v+1)\sqrt{v}}{\big((v+1)^2 + y\big)^2} \; = \; -\frac{\sqrt{-1-i\sqrt{y}}(\sqrt{y}+i)}{8\sqrt{y}(1+y)} With this form of the cut plane, we have i + i y = a + i b \sqrt{-i+i\sqrt{y}} = a + ib and 1 i y = a + i b \sqrt{-1-i\sqrt{y}} = -a + ib for some a , b > 0 a,b > 0 , and we deduce that A = π 2 y ( b y a ) A \; = \; \frac{\pi}{2\sqrt{y}}(b\sqrt{y}-a) We now calculate that a = 1 + y 1 2 b = y 2 ( 1 + y 1 ) a \; = \; \sqrt{\frac{\sqrt{1+y}-1}{2}} \hspace{2cm} b \; = \; \sqrt{\frac{y}{2\big(\sqrt{1+y}-1\big)}} obtaining, finally, A = π ( 1 + y ) ( 1 + y 1 ) 8 y A \; = \; \pi\sqrt{\frac{(1+y)\big(\sqrt{1+y}-1\big)}{8y}} Putting t = 1 + y t = \sqrt{1+y} , we need to solve the equation ( 1 + y ) ( 1 + y 1 ) = 8 y t 2 ( t 1 ) = 8 ( t 2 1 ) ( t 1 ) ( t 2 8 t 8 ) = 0 \begin{aligned} (1+y)\big(\sqrt{1+y}-1\big) & = \; 8y \\ t^2(t-1) & = \; 8(t^2-1) \\ (t-1)(t^2 - 8t - 8) & = \; 0 \end{aligned} Thus we must have t = 1 , 4 ± 2 6 t=1\,,\, 4 \pm 2\sqrt{6} . Since y > 0 y > 0 we must have t > 1 t > 1 and hence t = 4 + 2 6 t = 4+2\sqrt{6} . Thus y = 39 + 16 6 y = 39 + 16\sqrt{6} , and hence y 2 78 y = 15 y^2 - 78y = \boxed{15} .

Note that although t = 1 t=1 is a solution of the final cubic equation, y = 0 y=0 is not a solution to this problem. The formula for A A is not valid for y = 0 y=0 , and indeed the limit of A A as y 0 y \to 0 is 1 4 π \tfrac14\pi , which is the correct value of the integral when y = 0 y=0 , as an x = sin 4 θ x = \sin^4\theta substitution readily shows). By its definition, t t must be positive, and so the 4 2 6 4 - 2\sqrt{6} value must be discarded from this analysis.


That said, there is a second solution for y y , where 1 < y < 0 -1 < y < 0 . If we put y = p 2 y = -p^2 where 0 < p < 1 0 < p < 1 , then the poles for f y ( v ) f_y(v) are ( 1 + p ) -(1+p) and ( 1 p ) -(1-p) , both real and negative. The same contour integration can be attempted, and in this case we obtain A = π 1 p 2 4 p ( 1 + p 1 p ) A \; = \; \pi\frac{\sqrt{1-p^2}}{4p}\big(\sqrt{1+p} - \sqrt{1-p}\big) and we obtain A = π A=\pi when p 4 + 78 p 2 15 = 0 p^4 + 78p^2 - 15 = 0 , or y 2 78 y 15 = 0 y^2 - 78y - 15 = 0 , and hence we obtain a solution for y y with 1 < y < 0 -1 < y < 0 , namely y = 39 16 6 y = 39 - 16\sqrt{6} . This second solution still gives y 2 78 y = 15 y^2 - 78y = 15 , of course.

Prince Zarzees
Jan 31, 2020
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
from scipy.integrate  import quad                                                           
from math import sqrt,pi                                                           
def f(x,y):                                                           
    return sqrt(-1+sqrt((1+y)/x-y))                                                           
def bisection(lower,upper):                                                           
    mid=(lower+upper)/2                                                           
    if abs(quad(f,0,1,args=(mid,))[0]-pi)<=0.000000000001:                                                           
        return mid                                                           
    if quad(f,0,1,args=(mid,))[0]>pi:                                                           
        return bisection(lower,mid)                                                           
    elif quad(f,0,1,args=(mid,))[0]<pi:                                                           
        return bisection(mid,upper)                                                           
print (bisection(0,100)**2-78*bisection(0,100))                                                           
-  #You can try this python also :) 

In each iteration, your algorithm potentially evaluates

quad(f,0,1,args=(mid,))[0]

at least twice (except in the final iteration), and potentially three times, which must slow your routine down somewhat. Surely it would be better to assign this calculated value to a variable, and then test this variable (is it "equal to" π \pi , less than π \pi , or greater than π \pi )?

Your algorithm goes from (lower,upper) to (mid,upper) if the integral is less than π \pi , and to (mid,upper) if the integral is greater than π \pi . In other words, you are assuming that the integral is a increasing function of y y . Why? Given the exact answer, this is not at all obvious, although it is true.

Mark Hennings - 1 year, 4 months ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...