Hypergeometric Differential Equation
Gradshteyn & Ryzhik
According to §9.151 (p. 1045) of Gradshteyn & Ryzhik (1965) , "… a hypergeometric series is one of the solutions of the differential equation,
0
{\displaystyle 0}
=
{\displaystyle =}
z
(
1
−
z
)
d
2
u
d
z
2
+
[
γ
−
(
α
+
β
+
1
)
z
]
d
u
d
z
−
α
β
u
,
{\displaystyle z(1-z){\frac {d^{2}u}{dz^{2}}}+[\gamma -(\alpha +\beta +1)z]{\frac {du}{dz}}-\alpha \beta u\,,}
which is called the hypergeometric equation. And, according to §9.10 (p. 1039) of Gradshteyn & Ryzhik (1965) , "A hypergeometric series is a series of the form,
F
(
α
,
β
;
γ
;
z
)
{\displaystyle F(\alpha ,\beta ;\gamma ;z)}
=
{\displaystyle =}
1
+
[
α
⋅
β
γ
⋅
1
]
z
+
[
α
(
α
+
1
)
β
(
β
+
1
)
γ
(
γ
+
1
)
⋅
1
⋅
2
]
z
2
+
[
α
(
α
+
1
)
(
α
+
2
)
β
(
β
+
1
)
(
β
+
2
)
γ
(
γ
+
1
)
(
γ
+
2
)
⋅
1
⋅
2
⋅
3
]
z
3
+
…
{\displaystyle 1+{\biggl [}{\frac {\alpha \cdot \beta }{\gamma \cdot 1}}{\biggr ]}z+{\biggl [}{\frac {\alpha (\alpha +1)\beta (\beta +1)}{\gamma (\gamma +1)\cdot 1\cdot 2}}{\biggr ]}z^{2}+{\biggl [}{\frac {\alpha (\alpha +1)(\alpha +2)\beta (\beta +1)(\beta +2)}{\gamma (\gamma +1)(\gamma +2)\cdot 1\cdot 2\cdot 3}}{\biggr ]}z^{3}+\dots }
Among other attributes, Gradshteyn & Ryzhik (1965) note that this, "… series terminates if
α
{\displaystyle \alpha }
or
β
{\displaystyle \beta }
is equal to a negative integer or to zero."
Van der Borght
General Value for b
Comment by J. E. Tohline: In Van der Borght (1970), the fourth argument of the hypergeometric function is ax2 whereas the more general power law form, axb , applies. In association with his equation (3), 📚 R. Van der Borght (1970, Proc. Astr. Soc. Australia, Vol. 1, Issue 7, pp. 325 - 326) states that a displacement function of the form,
ξ
{\displaystyle \xi }
=
{\displaystyle =}
x
c
F
(
α
,
β
,
γ
,
a
x
b
)
,
{\displaystyle x^{c}F(\alpha ,\beta ,\gamma ,ax^{b})\,,}
provides a solution to the following 2nd -order ODE:
0
{\displaystyle 0}
=
{\displaystyle =}
x
2
(
a
x
b
−
1
)
⋅
d
2
ξ
d
x
2
+
(
a
p
x
b
+
λ
)
x
⋅
d
ξ
d
x
+
(
a
r
x
b
+
s
)
⋅
ξ
.
{\displaystyle x^{2}(ax^{b}-1)\cdot {\frac {d^{2}\xi }{dx^{2}}}+(apx^{b}+\lambda )x\cdot {\frac {d\xi }{dx}}+(arx^{b}+s)\cdot \xi \,.}
Is this ODE essentially the same as the above-defined hypergeometric equation ?
A mapping between the two differential equations requires,
a
x
b
{\displaystyle ax^{b}}
↔
{\displaystyle \leftrightarrow }
z
,
{\displaystyle z\,,}
and,
ξ
x
c
{\displaystyle {\frac {\xi }{x^{c}}}}
↔
{\displaystyle \leftrightarrow }
u
,
{\displaystyle u\,,}
⇒
d
x
d
z
{\displaystyle \Rightarrow ~~~{\frac {dx}{dz}}}
=
{\displaystyle =}
d
d
z
(
z
a
)
1
/
b
=
a
−
1
/
b
⋅
d
z
1
/
b
d
z
=
1
b
⋅
a
−
1
/
b
z
−
1
+
1
/
b
=
1
b
z
⋅
(
z
a
)
1
/
b
=
x
b
z
=
x
1
−
b
a
b
,
{\displaystyle {\frac {d}{dz}}{\biggl (}{\frac {z}{a}}{\biggr )}^{1/b}=a^{-1/b}\cdot {\frac {dz^{1/b}}{dz}}={\frac {1}{b}}\cdot a^{-1/b}z^{-1+1/b}={\frac {1}{bz}}\cdot {\biggl (}{\frac {z}{a}}{\biggr )}^{1/b}={\frac {x}{bz}}={\frac {x^{1-b}}{ab}}\,,}
in which case,
0
{\displaystyle 0}
=
{\displaystyle =}
z
(
1
−
z
)
d
2
u
d
z
2
+
[
γ
−
(
α
+
β
+
1
)
z
]
d
u
d
z
−
α
β
u
{\displaystyle z(1-z){\frac {d^{2}u}{dz^{2}}}+[\gamma -(\alpha +\beta +1)z]{\frac {du}{dz}}-\alpha \beta u}
=
{\displaystyle =}
a
x
b
(
1
−
a
x
b
)
{
d
x
d
z
⋅
d
d
x
[
d
x
d
z
⋅
d
d
x
(
ξ
x
−
c
)
]
}
+
[
γ
−
(
α
+
β
+
1
)
a
x
b
]
d
x
d
z
⋅
d
d
x
⋅
(
ξ
x
−
c
)
−
α
β
(
ξ
x
−
c
)
{\displaystyle ax^{b}(1-ax^{b}){\biggl \{}{\frac {dx}{dz}}\cdot {\frac {d}{dx}}{\biggl [}{\frac {dx}{dz}}\cdot {\frac {d}{dx}}{\biggl (}\xi x^{-c}{\biggr )}{\biggr ]}{\biggr \}}+[\gamma -(\alpha +\beta +1)ax^{b}]{\frac {dx}{dz}}\cdot {\frac {d}{dx}}\cdot {\biggl (}\xi x^{-c}{\biggr )}-\alpha \beta {\biggl (}\xi x^{-c}{\biggr )}}
=
{\displaystyle =}
a
x
b
(
1
−
a
x
b
)
{
x
1
−
b
a
b
⋅
d
d
x
[
x
1
−
b
a
b
⋅
(
x
−
c
d
ξ
d
x
−
c
ξ
x
−
1
−
c
)
]
}
+
[
γ
−
(
α
+
β
+
1
)
a
x
b
]
⋅
x
1
−
b
a
b
⋅
(
x
−
c
d
ξ
d
x
−
c
ξ
x
−
1
−
c
)
−
α
β
(
ξ
x
−
c
)
{\displaystyle ax^{b}(1-ax^{b}){\biggl \{}{\frac {x^{1-b}}{ab}}\cdot {\frac {d}{dx}}{\biggl [}{\frac {x^{1-b}}{ab}}\cdot {\biggl (}x^{-c}{\frac {d\xi }{dx}}-c\xi x^{-1-c}{\biggr )}{\biggr ]}{\biggr \}}+[\gamma -(\alpha +\beta +1)ax^{b}]\cdot {\frac {x^{1-b}}{ab}}\cdot {\biggl (}x^{-c}{\frac {d\xi }{dx}}-c\xi x^{-1-c}{\biggr )}-\alpha \beta {\biggl (}\xi x^{-c}{\biggr )}}
=
{\displaystyle =}
a
x
b
(
1
−
a
x
b
)
{
x
1
−
b
a
2
b
2
⋅
d
d
x
[
x
1
−
b
−
c
d
ξ
d
x
−
c
ξ
x
−
b
−
c
]
}
+
[
γ
−
(
α
+
β
+
1
)
a
x
b
]
a
b
⋅
(
x
1
−
b
−
c
d
ξ
d
x
−
c
ξ
x
−
b
−
c
)
−
α
β
(
ξ
x
−
c
)
{\displaystyle ax^{b}(1-ax^{b}){\biggl \{}{\frac {x^{1-b}}{a^{2}b^{2}}}\cdot {\frac {d}{dx}}{\biggl [}x^{1-b-c}{\frac {d\xi }{dx}}-c\xi x^{-b-c}{\biggr ]}{\biggr \}}+{\frac {[\gamma -(\alpha +\beta +1)ax^{b}]}{ab}}\cdot {\biggl (}x^{1-b-c}{\frac {d\xi }{dx}}-c\xi x^{-b-c}{\biggr )}-\alpha \beta {\biggl (}\xi x^{-c}{\biggr )}}
=
{\displaystyle =}
a
x
b
(
1
−
a
x
b
)
⋅
x
1
−
b
a
2
b
2
⋅
[
x
1
−
b
−
c
d
2
ξ
d
x
2
+
(
1
−
b
−
c
)
x
−
b
−
c
d
ξ
d
x
−
c
x
−
b
−
c
d
ξ
d
x
−
c
(
−
b
−
c
)
ξ
x
−
1
−
b
−
c
]
{\displaystyle ax^{b}(1-ax^{b})\cdot {\frac {x^{1-b}}{a^{2}b^{2}}}\cdot {\biggl [}x^{1-b-c}{\frac {d^{2}\xi }{dx^{2}}}+(1-b-c)x^{-b-c}{\frac {d\xi }{dx}}-cx^{-b-c}{\frac {d\xi }{dx}}-c(-b-c)\xi x^{-1-b-c}{\biggr ]}}
+
[
γ
−
(
α
+
β
+
1
)
a
x
b
]
a
b
⋅
[
x
1
−
b
−
c
d
ξ
d
x
−
c
ξ
x
−
b
−
c
]
−
(
α
β
x
−
c
)
ξ
{\displaystyle +{\frac {[\gamma -(\alpha +\beta +1)ax^{b}]}{ab}}\cdot {\biggl [}x^{1-b-c}{\frac {d\xi }{dx}}-c\xi x^{-b-c}{\biggr ]}-{\biggl (}\alpha \beta x^{-c}{\biggr )}\xi }
=
{\displaystyle =}
(
1
−
a
x
b
)
⋅
x
a
b
2
⋅
[
x
1
−
b
−
c
d
2
ξ
d
x
2
]
+
(
1
−
a
x
b
)
⋅
x
a
b
2
⋅
[
(
1
−
b
−
c
)
x
−
b
−
c
d
ξ
d
x
−
c
x
−
b
−
c
d
ξ
d
x
]
+
(
1
−
a
x
b
)
⋅
x
a
b
2
⋅
[
c
(
b
+
c
)
ξ
x
−
1
−
b
−
c
]
{\displaystyle (1-ax^{b})\cdot {\frac {x}{ab^{2}}}\cdot {\biggl [}x^{1-b-c}{\frac {d^{2}\xi }{dx^{2}}}{\biggr ]}+(1-ax^{b})\cdot {\frac {x}{ab^{2}}}\cdot {\biggl [}(1-b-c)x^{-b-c}{\frac {d\xi }{dx}}-cx^{-b-c}{\frac {d\xi }{dx}}{\biggr ]}+(1-ax^{b})\cdot {\frac {x}{ab^{2}}}\cdot {\biggl [}c(b+c)\xi x^{-1-b-c}{\biggr ]}}
+
[
γ
−
(
α
+
β
+
1
)
a
x
b
]
a
b
⋅
[
x
1
−
b
−
c
d
ξ
d
x
]
−
[
γ
−
(
α
+
β
+
1
)
a
x
b
]
a
b
⋅
[
c
ξ
x
−
b
−
c
]
−
(
α
β
x
−
c
)
ξ
{\displaystyle +{\frac {[\gamma -(\alpha +\beta +1)ax^{b}]}{ab}}\cdot {\biggl [}x^{1-b-c}{\frac {d\xi }{dx}}{\biggr ]}-{\frac {[\gamma -(\alpha +\beta +1)ax^{b}]}{ab}}\cdot {\biggl [}c\xi x^{-b-c}{\biggr ]}-{\biggl (}\alpha \beta x^{-c}{\biggr )}\xi }
=
{\displaystyle =}
(
1
−
a
x
b
)
a
b
2
⋅
x
2
[
x
−
b
−
c
d
2
ξ
d
x
2
]
+
{
(
1
−
a
x
b
)
a
b
2
⋅
[
(
1
−
b
−
c
)
x
1
−
b
−
c
−
c
x
1
−
b
−
c
]
+
[
γ
−
(
α
+
β
+
1
)
a
x
b
]
a
b
⋅
[
x
1
−
b
−
c
]
}
d
ξ
d
x
{\displaystyle {\frac {(1-ax^{b})}{ab^{2}}}\cdot x^{2}{\biggl [}x^{-b-c}{\frac {d^{2}\xi }{dx^{2}}}{\biggr ]}+{\biggl \{}{\frac {(1-ax^{b})}{ab^{2}}}\cdot {\biggl [}(1-b-c)x^{1-b-c}-cx^{1-b-c}{\biggr ]}+{\frac {[\gamma -(\alpha +\beta +1)ax^{b}]}{ab}}\cdot {\biggl [}x^{1-b-c}{\biggr ]}{\biggr \}}{\frac {d\xi }{dx}}}
+
{
−
c
[
γ
−
(
α
+
β
+
1
)
a
x
b
]
a
b
⋅
[
x
−
b
−
c
]
−
(
α
β
x
−
b
−
c
)
x
b
+
(
1
−
a
x
b
)
a
b
2
⋅
[
c
(
b
+
c
)
x
−
b
−
c
]
}
ξ
.
{\displaystyle +{\biggl \{}-{\frac {c[\gamma -(\alpha +\beta +1)ax^{b}]}{ab}}\cdot {\biggl [}x^{-b-c}{\biggr ]}-{\biggl (}\alpha \beta x^{-b-c}{\biggr )}x^{b}+{\frac {(1-ax^{b})}{ab^{2}}}\cdot {\biggl [}c(b+c)x^{-b-c}{\biggr ]}{\biggr \}}\xi \,.}
Multiplying through by
(
−
a
b
2
x
b
+
c
)
{\displaystyle (-ab^{2}x^{b+c})}
gives,
(
−
a
b
2
x
b
+
c
)
×
0
{\displaystyle (-ab^{2}x^{b+c})\times 0}
=
{\displaystyle =}
x
2
(
a
x
b
−
1
)
d
2
ξ
d
x
2
+
{
(
a
x
b
−
1
)
⋅
[
(
1
−
b
−
c
)
−
c
]
−
b
[
γ
−
(
α
+
β
+
1
)
a
x
b
]
}
x
⋅
d
ξ
d
x
{\displaystyle x^{2}(ax^{b}-1){\frac {d^{2}\xi }{dx^{2}}}+{\biggl \{}(ax^{b}-1)\cdot {\biggl [}(1-b-c)-c{\biggr ]}-b[\gamma -(\alpha +\beta +1)ax^{b}]{\biggr \}}x\cdot {\frac {d\xi }{dx}}}
+
{
b
c
[
γ
−
(
α
+
β
+
1
)
a
x
b
]
+
(
α
β
a
b
2
)
x
b
−
(
1
−
a
x
b
)
⋅
[
c
(
b
+
c
)
]
}
ξ
{\displaystyle +{\biggl \{}bc[\gamma -(\alpha +\beta +1)ax^{b}]+{\biggl (}\alpha \beta ab^{2}{\biggr )}x^{b}-(1-ax^{b})\cdot {\biggl [}c(b+c){\biggr ]}{\biggr \}}\xi }
=
{\displaystyle =}
x
2
(
a
x
b
−
1
)
d
2
ξ
d
x
2
+
{
−
(
1
+
b
γ
−
b
−
2
c
)
+
a
x
b
(
1
−
b
−
2
c
)
+
b
(
α
+
β
+
1
)
a
x
b
}
x
⋅
d
ξ
d
x
{\displaystyle x^{2}(ax^{b}-1){\frac {d^{2}\xi }{dx^{2}}}+{\biggl \{}-(1+b\gamma -b-2c)+ax^{b}(1-b-2c)+b(\alpha +\beta +1)ax^{b}{\biggr \}}x\cdot {\frac {d\xi }{dx}}}
+
{
b
c
γ
−
c
(
b
+
c
)
−
b
c
(
α
+
β
+
1
)
a
x
b
+
α
β
a
b
2
x
b
+
c
(
b
+
c
)
a
x
b
}
ξ
{\displaystyle +{\biggl \{}bc\gamma -c(b+c)-bc(\alpha +\beta +1)ax^{b}+\alpha \beta ab^{2}x^{b}+c(b+c)ax^{b}{\biggr \}}\xi }
=
{\displaystyle =}
x
2
(
a
x
b
−
1
)
d
2
ξ
d
x
2
+
(
a
p
x
b
+
λ
)
x
⋅
d
ξ
d
x
+
(
a
r
x
b
+
s
)
ξ
,
{\displaystyle x^{2}(ax^{b}-1){\frac {d^{2}\xi }{dx^{2}}}+(apx^{b}+\lambda )x\cdot {\frac {d\xi }{dx}}+(arx^{b}+s)\xi \,,}
which matches equation (3) of 📚 Van der Borght (1970) if the expressions for the four new scalar coefficients are,
1st :
p
{\displaystyle p}
=
{\displaystyle =}
(
1
−
b
−
2
c
)
+
b
(
α
+
β
+
1
)
=
1
−
2
c
+
b
(
α
+
β
)
,
{\displaystyle (1-b-2c)+b(\alpha +\beta +1)=1-2c+b(\alpha +\beta )\,,}
2nd :
λ
{\displaystyle \lambda }
=
{\displaystyle =}
−
(
1
+
b
γ
−
b
−
2
c
)
,
{\displaystyle -(1+b\gamma -b-2c)\,,}
3rd :
s
{\displaystyle s}
=
{\displaystyle =}
b
c
γ
−
c
(
b
+
c
)
=
c
(
b
γ
−
b
−
c
)
,
{\displaystyle bc\gamma -c(b+c)=c(b\gamma -b-c)\,,}
4th :
r
{\displaystyle r}
=
{\displaystyle =}
−
b
c
(
α
+
β
+
1
)
+
α
β
b
2
+
c
(
b
+
c
)
.
{\displaystyle -bc(\alpha +\beta +1)+\alpha \beta b^{2}+c(b+c)\,.}
If, for any given problem, we are given the values of these four scalar coefficients along with a choice of the exponent,
b
{\displaystyle b}
, that appears in the fourth argument of the hypergeometric series, we can determine values the other three arguments of the hypergeometric series —
α
,
β
,
γ
{\displaystyle \alpha ,\beta ,\gamma }
— and the exponent,
c
{\displaystyle c}
. In what follows we show how this is done.
Determining the Value of c
i.) Determining the value of the exponent,
c
{\displaystyle c}
, from
λ
{\displaystyle \lambda }
and
s
{\displaystyle s}
Equating
b
γ
{\displaystyle b\gamma }
in the 2nd and 3rd of these expressions, gives,
−
(
1
−
b
−
2
c
+
λ
)
{\displaystyle -(1-b-2c+\lambda )}
=
{\displaystyle =}
s
c
+
b
+
c
{\displaystyle {\frac {s}{c}}+b+c}
⇒
(
1
+
λ
)
{\displaystyle \Rightarrow ~~~(1+\lambda )}
=
{\displaystyle =}
c
−
s
c
{\displaystyle c-{\frac {s}{c}}}
⇒
0
{\displaystyle \Rightarrow ~~~0}
=
{\displaystyle =}
c
2
−
c
(
1
+
λ
)
−
s
.
{\displaystyle c^{2}-c(1+\lambda )-s\,.}
The pair of roots,
c
±
{\displaystyle c_{\pm }}
, of this quadratic equation are then obtained from the relation,
2
c
±
{\displaystyle 2c_{\pm }}
=
{\displaystyle =}
(
1
+
λ
)
±
[
(
1
+
λ
)
2
+
4
s
]
1
/
2
.
{\displaystyle (1+\lambda )\pm [(1+\lambda )^{2}+4s]^{1/2}\,.}
Note for further use below that,
2
(
c
+
−
c
−
)
{\displaystyle 2(c_{+}-c_{-})}
=
{\displaystyle =}
{
(
1
+
λ
)
+
[
(
1
+
λ
)
2
+
4
s
]
1
/
2
}
−
{
(
1
+
λ
)
−
[
(
1
+
λ
)
2
+
4
s
]
1
/
2
}
.
{\displaystyle \{(1+\lambda )+[(1+\lambda )^{2}+4s]^{1/2}\}-\{(1+\lambda )-[(1+\lambda )^{2}+4s]^{1/2}\}\,.}
=
{\displaystyle =}
2
[
(
1
+
λ
)
2
+
4
s
]
1
/
2
.
{\displaystyle 2[(1+\lambda )^{2}+4s]^{1/2}\,.}
Consistent with our derivation, 📚 Van der Borght (1970) states, "… if
A
1
{\displaystyle A_{1}}
,
A
2
{\displaystyle A_{2}}
are the solutions of
A
2
−
(
λ
+
1
)
A
−
s
=
0
{\displaystyle A^{2}-(\lambda +1)A-s=0}
… then
c
=
A
1
{\displaystyle c=A_{1}}
…"
Determining the Value of γ
ii.) Determining the value of,
b
(
1
−
γ
)
{\displaystyle b(1-\gamma )}
, from
λ
{\displaystyle \lambda }
and
s
{\displaystyle s}
Combining our 3rd expression with the quadratic equation for
c
±
{\displaystyle c_{\pm }}
in such a way as to eliminate
s
{\displaystyle s}
, we find,
c
b
(
γ
−
1
)
−
c
2
{\displaystyle cb(\gamma -1)-c^{2}}
=
{\displaystyle =}
c
2
−
c
(
1
+
λ
)
{\displaystyle c^{2}-c(1+\lambda )}
⇒
b
(
γ
−
1
)
{\displaystyle \Rightarrow ~~~b(\gamma -1)}
=
{\displaystyle =}
2
c
±
−
(
1
+
λ
)
{\displaystyle 2c_{\pm }-(1+\lambda )}
Adopting the superior sign, we find that,
b
(
γ
−
1
)
{\displaystyle b(\gamma -1)}
=
{\displaystyle =}
2
c
+
−
(
1
+
λ
)
=
[
(
1
+
λ
)
2
+
4
s
]
1
/
2
{\displaystyle 2c_{+}-(1+\lambda )=[(1+\lambda )^{2}+4s]^{1/2}}
=
{\displaystyle =}
(
c
+
−
c
−
)
,
{\displaystyle (c_{+}-c_{-})\,,}
where, in order to make this last step we have drawn from the relation derived in box "i." immediately above.
Consistent with this derivation, 📚 Van der Borght (1970) states, "…
(
1
−
γ
)
b
=
A
2
−
A
1
{\displaystyle (1-\gamma )b=A_{2}-A_{1}}
…"
Determining the Value of β
iii.) Determining
b
β
{\displaystyle b\beta }
Combining our 1st and 4th expressions in such a way as to cancel terms involving
b
c
(
α
+
β
)
{\displaystyle bc(\alpha +\beta )}
, we find,
(
1
−
p
)
c
−
2
c
2
+
b
c
(
α
+
β
)
{\displaystyle (1-p)c-2c^{2}+bc(\alpha +\beta )}
=
{\displaystyle =}
b
c
(
α
+
β
)
+
b
c
−
α
β
b
2
−
c
(
b
+
c
)
+
r
{\displaystyle bc(\alpha +\beta )+bc-\alpha \beta b^{2}-c(b+c)+r}
⇒
(
1
−
p
)
c
−
c
2
{\displaystyle \Rightarrow ~~~(1-p)c-c^{2}}
=
{\displaystyle =}
−
α
β
b
2
+
r
{\displaystyle -\alpha \beta b^{2}+r}
⇒
(
b
α
)
(
b
β
)
{\displaystyle \Rightarrow ~~~(b\alpha )(b\beta )}
=
{\displaystyle =}
−
(
1
−
p
)
c
+
c
2
+
r
.
{\displaystyle -(1-p)c+c^{2}+r\,.}
Also, from the 1st expression alone we can write,
(
b
α
)
{\displaystyle (b\alpha )}
=
{\displaystyle =}
(
p
−
1
)
+
2
c
−
b
β
.
{\displaystyle (p-1)+2c-b\beta \,.}
Together, then, we have,
(
b
β
)
[
(
p
−
1
)
+
2
c
−
b
β
]
{\displaystyle (b\beta ){\biggl [}(p-1)+2c-b\beta {\biggr ]}}
=
{\displaystyle =}
−
(
1
−
p
)
c
+
c
2
+
r
{\displaystyle -(1-p)c+c^{2}+r}
⇒
0
{\displaystyle \Rightarrow ~~~0}
=
{\displaystyle =}
(
b
β
)
2
−
(
b
β
)
[
(
p
−
1
)
+
2
c
]
⏟
(
b
α
+
b
β
)
+
[
c
2
+
(
p
−
1
)
c
+
r
]
⏞
(
b
α
)
⋅
(
b
β
)
{\displaystyle (b\beta )^{2}-(b\beta )\underbrace {[(p-1)+2c]} _{(b\alpha +b\beta )}+\overbrace {[c^{2}+(p-1)c+r]} ^{(b\alpha )\cdot (b\beta )}}
The pair of roots,
(
b
β
)
±
{\displaystyle (b\beta )_{\pm }}
, of this quadratic equation are then obtained from the relation,
2
(
b
β
)
±
{\displaystyle 2(b\beta )_{\pm }}
=
{\displaystyle =}
[
(
p
−
1
)
+
2
c
]
±
{
[
(
p
−
1
)
+
2
c
]
2
−
4
[
c
2
+
(
p
−
1
)
c
+
r
]
}
1
/
2
{\displaystyle [(p-1)+2c]\pm {\biggl \{}[(p-1)+2c]^{2}-4[c^{2}+(p-1)c+r]{\biggr \}}^{1/2}}
📚 Van der Borght (1970) states that if, "…
B
1
,
B
2
{\displaystyle B_{1},B_{2}}
are solutions of
B
2
−
(
p
−
1
)
B
+
r
=
0
{\displaystyle B^{2}-(p-1)B+r=0}
, then …
b
β
=
A
1
+
B
2
.
{\displaystyle b\beta =A_{1}+B_{2}.}
" Let's see if our derivation leads to this same conclusion. First note that the roots,
B
±
{\displaystyle B_{\pm }}
, of this Van der Borght quadratic equation are,
2
B
±
{\displaystyle 2B_{\pm }}
=
{\displaystyle =}
(
p
−
1
)
±
[
(
p
−
1
)
2
−
4
r
]
1
/
2
.
{\displaystyle (p-1)\pm [(p-1)^{2}-4r]^{1/2}\,.}
If we assign the inferior root with Van der Borght's notation,
B
2
{\displaystyle B_{2}}
, then,
2
(
A
1
+
B
2
)
=
2
(
c
+
+
B
−
)
{\displaystyle 2(A_{1}+B_{2})=2(c_{+}+B_{-})}
=
{\displaystyle =}
{
(
1
+
λ
)
+
[
(
1
+
λ
)
2
+
4
s
]
1
/
2
}
+
{
(
p
−
1
)
−
[
(
p
−
1
)
2
−
4
r
]
1
/
2
}
{\displaystyle \{(1+\lambda )+[(1+\lambda )^{2}+4s]^{1/2}\}+\{(p-1)-[(p-1)^{2}-4r]^{1/2}\}}
Completing the Squares to Determine the Value of α
First, Complete the square in the quadratic equation for
B
2
{\displaystyle B^{2}}
:
B
2
−
(
p
−
1
)
B
{\displaystyle B^{2}-(p-1)B}
=
{\displaystyle =}
−
r
{\displaystyle -r}
⇒
B
2
−
(
p
−
1
)
B
+
[
(
p
−
1
)
2
]
2
{\displaystyle \Rightarrow ~~~B^{2}-(p-1)B+{\biggl [}{\frac {(p-1)}{2}}{\biggr ]}^{2}}
=
{\displaystyle =}
[
(
p
−
1
)
2
]
2
−
r
{\displaystyle {\biggl [}{\frac {(p-1)}{2}}{\biggr ]}^{2}-r}
⇒
[
B
−
(
p
−
1
)
2
]
2
{\displaystyle \Rightarrow ~~~{\biggl [}B-{\frac {(p-1)}{2}}{\biggr ]}^{2}}
=
{\displaystyle =}
[
(
p
−
1
)
2
]
2
−
r
{\displaystyle {\biggl [}{\frac {(p-1)}{2}}{\biggr ]}^{2}-r}
Second, complete the square in the quadratic equation for
A
2
{\displaystyle A^{2}}
— which also completes the square for
c
2
{\displaystyle c^{2}}
:
A
2
−
(
λ
+
1
)
A
{\displaystyle A^{2}-(\lambda +1)A}
=
{\displaystyle =}
s
{\displaystyle s}
⇒
A
2
−
(
λ
+
1
)
A
+
[
(
λ
+
1
)
2
]
2
{\displaystyle \Rightarrow ~~~A^{2}-(\lambda +1)A+{\biggl [}{\frac {(\lambda +1)}{2}}{\biggr ]}^{2}}
=
{\displaystyle =}
s
+
[
(
λ
+
1
)
2
]
2
{\displaystyle s+{\biggl [}{\frac {(\lambda +1)}{2}}{\biggr ]}^{2}}
⇒
[
A
−
(
λ
+
1
)
2
]
2
{\displaystyle \Rightarrow ~~~{\biggl [}A-{\frac {(\lambda +1)}{2}}{\biggr ]}^{2}}
=
{\displaystyle =}
s
+
[
(
λ
+
1
)
2
]
2
{\displaystyle s+{\biggl [}{\frac {(\lambda +1)}{2}}{\biggr ]}^{2}}
Third, complete the square in the quadratic equation for
(
b
β
)
2
{\displaystyle (b\beta )^{2}}
:
(
b
β
)
2
−
(
b
β
)
[
(
p
−
1
)
+
2
c
]
{\displaystyle (b\beta )^{2}-(b\beta )[(p-1)+2c]}
=
{\displaystyle =}
−
[
c
2
+
(
p
−
1
)
c
+
r
]
{\displaystyle -[c^{2}+(p-1)c+r]}
⇒
(
b
β
)
2
−
(
b
β
)
[
(
p
−
1
)
+
2
c
]
+
[
(
p
−
1
)
+
2
c
2
]
2
{\displaystyle \Rightarrow ~~~(b\beta )^{2}-(b\beta )[(p-1)+2c]+{\biggl [}{\frac {(p-1)+2c}{2}}{\biggr ]}^{2}}
=
{\displaystyle =}
[
(
p
−
1
)
+
2
c
2
]
2
−
[
c
2
+
(
p
−
1
)
c
+
r
]
{\displaystyle {\biggl [}{\frac {(p-1)+2c}{2}}{\biggr ]}^{2}-[c^{2}+(p-1)c+r]}
⇒
[
(
b
β
)
−
(
p
−
1
)
+
2
c
2
]
2
{\displaystyle \Rightarrow ~~~{\biggl [}(b\beta )-{\frac {(p-1)+2c}{2}}{\biggr ]}^{2}}
=
{\displaystyle =}
[
(
p
−
1
)
+
2
c
2
]
2
−
[
c
2
+
(
p
−
1
)
c
+
(
p
−
1
)
2
2
2
]
−
r
+
(
p
−
1
)
2
2
2
{\displaystyle {\biggl [}{\frac {(p-1)+2c}{2}}{\biggr ]}^{2}-{\biggl [}c^{2}+(p-1)c+{\frac {(p-1)^{2}}{2^{2}}}{\biggr ]}-r+{\frac {(p-1)^{2}}{2^{2}}}}
=
{\displaystyle =}
[
(
p
−
1
)
2
+
c
]
2
−
[
c
+
(
p
−
1
)
2
]
2
−
r
+
(
p
−
1
)
2
2
2
{\displaystyle {\biggl [}{\frac {(p-1)}{2}}+c{\biggr ]}^{2}-{\biggl [}c+{\frac {(p-1)}{2}}{\biggr ]}^{2}-r+{\frac {(p-1)^{2}}{2^{2}}}}
=
{\displaystyle =}
−
r
+
(
p
−
1
)
2
2
2
{\displaystyle -r+{\frac {(p-1)^{2}}{2^{2}}}}
=
{\displaystyle =}
[
B
−
(
p
−
1
)
2
]
2
.
{\displaystyle {\biggl [}B-{\frac {(p-1)}{2}}{\biggr ]}^{2}\,.}
Taking the positive root of both sides of this expression, we find that,
(
b
β
)
−
(
p
−
1
)
2
−
c
+
{\displaystyle (b\beta )-{\frac {(p-1)}{2}}-c_{+}}
=
{\displaystyle =}
B
+
−
(
p
−
1
)
2
{\displaystyle B_{+}-{\frac {(p-1)}{2}}}
⇒
(
b
β
)
{\displaystyle \Rightarrow ~~~(b\beta )}
=
{\displaystyle =}
B
+
+
c
+
.
{\displaystyle B_{+}+c_{+}\,.}
But,
c
±
=
A
±
{\displaystyle c_{\pm }=A_{\pm }}
. So we conclude, as did 📚 Van der Borght (1970) , that,
(
b
β
)
{\displaystyle (b\beta )}
=
{\displaystyle =}
B
+
+
A
+
.
{\displaystyle B_{+}+A_{+}\,.}
Alternatively, taking the negative root of the RHS of this expression, we find that,
(
b
β
)
−
(
p
−
1
)
2
−
c
+
{\displaystyle (b\beta )-{\frac {(p-1)}{2}}-c_{+}}
=
{\displaystyle =}
−
B
−
+
(
p
−
1
)
2
{\displaystyle -B_{-}+{\frac {(p-1)}{2}}}
⇒
(
b
β
)
{\displaystyle \Rightarrow ~~~(b\beta )}
=
{\displaystyle =}
(
p
−
1
)
2
+
c
+
−
B
−
+
(
p
−
1
)
2
{\displaystyle {\frac {(p-1)}{2}}+c_{+}-B_{-}+{\frac {(p-1)}{2}}}
=
{\displaystyle =}
c
+
−
B
−
+
(
p
−
1
)
.
{\displaystyle c_{+}-B_{-}+(p-1)\,.}
Also, given that,
1
−
2
c
+
b
(
α
+
β
)
{\displaystyle 1-2c+b(\alpha +\beta )}
=
{\displaystyle =}
p
{\displaystyle p}
⇒
(
b
α
)
{\displaystyle \Rightarrow ~~~(b\alpha )}
=
{\displaystyle =}
p
−
1
+
2
c
−
(
b
β
)
{\displaystyle p-1+2c-(b\beta )}
=
{\displaystyle =}
(
p
−
1
)
+
2
c
−
[
c
+
−
B
−
+
(
p
−
1
)
]
{\displaystyle (p-1)+2c-{\biggl [}c_{+}-B_{-}+(p-1){\biggr ]}}
=
{\displaystyle =}
2
c
−
c
+
+
B
−
.
{\displaystyle 2c-c_{+}+B_{-}\,.}
As long as we assume that
c
=
c
+
{\displaystyle c=c_{+}}
in this expression, we also obtain the 📚 Van der Borght (1970) expression for
(
b
α
)
{\displaystyle (b\alpha )}
, namely,
(
b
α
)
{\displaystyle (b\alpha )}
=
{\displaystyle =}
B
−
+
A
+
.
{\displaystyle B_{-}+A_{+}\,.}
If b = 2
0
{\displaystyle 0}
=
{\displaystyle =}
z
(
1
−
z
)
d
2
u
d
z
2
+
[
γ
−
(
α
+
β
+
1
)
z
]
d
u
d
z
−
α
β
u
{\displaystyle z(1-z){\frac {d^{2}u}{dz^{2}}}+[\gamma -(\alpha +\beta +1)z]{\frac {du}{dz}}-\alpha \beta u}
=
{\displaystyle =}
a
x
2
(
1
−
a
x
2
)
{
d
x
d
z
⋅
d
d
x
[
d
x
d
z
⋅
d
d
x
(
ξ
x
−
c
)
]
}
+
[
γ
−
(
α
+
β
+
1
)
a
x
2
]
d
x
d
z
⋅
d
d
x
⋅
(
ξ
x
−
c
)
−
α
β
(
ξ
x
−
c
)
{\displaystyle ax^{2}(1-ax^{2}){\biggl \{}{\frac {dx}{dz}}\cdot {\frac {d}{dx}}{\biggl [}{\frac {dx}{dz}}\cdot {\frac {d}{dx}}{\biggl (}\xi x^{-c}{\biggr )}{\biggr ]}{\biggr \}}+[\gamma -(\alpha +\beta +1)ax^{2}]{\frac {dx}{dz}}\cdot {\frac {d}{dx}}\cdot {\biggl (}\xi x^{-c}{\biggr )}-\alpha \beta {\biggl (}\xi x^{-c}{\biggr )}}
=
{\displaystyle =}
a
x
2
(
1
−
a
x
2
)
{
1
2
a
x
⋅
d
d
x
[
1
2
a
x
⋅
(
x
−
c
d
ξ
d
x
−
c
ξ
x
−
1
−
c
)
]
}
+
[
γ
−
(
α
+
β
+
1
)
a
x
2
]
⋅
1
2
a
x
⋅
(
x
−
c
d
ξ
d
x
−
c
ξ
x
−
1
−
c
)
−
α
β
(
ξ
x
−
c
)
{\displaystyle ax^{2}(1-ax^{2}){\biggl \{}{\frac {1}{2ax}}\cdot {\frac {d}{dx}}{\biggl [}{\frac {1}{2ax}}\cdot {\biggl (}x^{-c}{\frac {d\xi }{dx}}-c\xi x^{-1-c}{\biggr )}{\biggr ]}{\biggr \}}+[\gamma -(\alpha +\beta +1)ax^{2}]\cdot {\frac {1}{2ax}}\cdot {\biggl (}x^{-c}{\frac {d\xi }{dx}}-c\xi x^{-1-c}{\biggr )}-\alpha \beta {\biggl (}\xi x^{-c}{\biggr )}}
=
{\displaystyle =}
a
x
2
(
1
−
a
x
2
)
{
1
4
a
2
x
⋅
d
d
x
[
(
x
−
1
−
c
d
ξ
d
x
−
c
ξ
x
−
2
−
c
)
]
}
+
[
γ
−
(
α
+
β
+
1
)
a
x
2
]
⋅
1
2
a
⋅
(
x
−
1
−
c
d
ξ
d
x
−
c
ξ
x
−
2
−
c
)
−
α
β
(
ξ
x
−
c
)
{\displaystyle ax^{2}(1-ax^{2}){\biggl \{}{\frac {1}{4a^{2}x}}\cdot {\frac {d}{dx}}{\biggl [}{\biggl (}x^{-1-c}{\frac {d\xi }{dx}}-c\xi x^{-2-c}{\biggr )}{\biggr ]}{\biggr \}}+[\gamma -(\alpha +\beta +1)ax^{2}]\cdot {\frac {1}{2a}}\cdot {\biggl (}x^{-1-c}{\frac {d\xi }{dx}}-c\xi x^{-2-c}{\biggr )}-\alpha \beta {\biggl (}\xi x^{-c}{\biggr )}}
=
{\displaystyle =}
x
(
1
−
a
x
2
)
4
a
{
x
−
1
−
c
d
2
ξ
d
x
2
−
(
1
+
c
)
x
−
2
−
c
d
ξ
d
x
−
c
x
−
2
−
c
d
ξ
d
x
+
(
2
c
+
c
2
)
ξ
x
−
3
−
c
}
{\displaystyle {\frac {x(1-ax^{2})}{4a}}{\biggl \{}x^{-1-c}{\frac {d^{2}\xi }{dx^{2}}}-(1+c)x^{-2-c}{\frac {d\xi }{dx}}-cx^{-2-c}{\frac {d\xi }{dx}}+(2c+c^{2})\xi x^{-3-c}{\biggr \}}}
+
[
γ
−
(
α
+
β
+
1
)
a
x
2
2
a
]
⋅
(
x
−
1
−
c
d
ξ
d
x
)
−
{
[
γ
−
(
α
+
β
+
1
)
a
x
2
2
a
]
⋅
c
x
−
2
−
c
+
α
β
x
−
c
}
ξ
{\displaystyle +{\biggl [}{\frac {\gamma -(\alpha +\beta +1)ax^{2}}{2a}}{\biggr ]}\cdot {\biggl (}x^{-1-c}{\frac {d\xi }{dx}}{\biggr )}-{\biggl \{}{\biggl [}{\frac {\gamma -(\alpha +\beta +1)ax^{2}}{2a}}{\biggr ]}\cdot cx^{-2-c}+\alpha \beta x^{-c}{\biggr \}}\xi }
=
{\displaystyle =}
x
(
1
−
a
x
2
)
4
a
{
x
−
1
−
c
d
2
ξ
d
x
2
}
+
{
x
(
1
−
a
x
2
)
4
a
[
−
(
1
+
c
)
x
−
2
−
c
−
c
x
−
2
−
c
]
+
[
γ
−
(
α
+
β
+
1
)
a
x
2
2
a
]
⋅
x
−
1
−
c
}
d
ξ
d
x
{\displaystyle {\frac {x(1-ax^{2})}{4a}}{\biggl \{}x^{-1-c}{\frac {d^{2}\xi }{dx^{2}}}{\biggr \}}+{\biggl \{}{\frac {x(1-ax^{2})}{4a}}{\biggl [}-(1+c)x^{-2-c}-cx^{-2-c}{\biggr ]}+{\biggl [}{\frac {\gamma -(\alpha +\beta +1)ax^{2}}{2a}}{\biggr ]}\cdot x^{-1-c}{\biggr \}}{\frac {d\xi }{dx}}}
+
{
[
x
(
1
−
a
x
2
)
4
a
]
(
2
c
+
c
2
)
x
−
3
−
c
−
[
γ
−
(
α
+
β
+
1
)
a
x
2
2
a
]
⋅
c
x
−
2
−
c
−
α
β
x
−
c
}
ξ
.
{\displaystyle +{\biggl \{}{\biggl [}{\frac {x(1-ax^{2})}{4a}}{\biggr ]}(2c+c^{2})x^{-3-c}-{\biggl [}{\frac {\gamma -(\alpha +\beta +1)ax^{2}}{2a}}{\biggr ]}\cdot cx^{-2-c}-\alpha \beta x^{-c}{\biggr \}}\xi \,.}
Multiplying through by a term proportional to
x
2
+
c
{\displaystyle x^{2+c}}
gives us,
[
−
4
a
x
2
+
c
]
×
[
0
]
{\displaystyle {\biggl [}-4ax^{2+c}{\biggr ]}\times {\biggl [}0{\biggr ]}}
=
{\displaystyle =}
x
2
(
a
x
2
−
1
)
{
d
2
ξ
d
x
2
}
+
{
x
(
1
−
a
x
2
)
[
(
1
+
c
)
+
c
]
−
2
[
γ
−
(
α
+
β
+
1
)
a
x
2
]
⋅
x
}
d
ξ
d
x
{\displaystyle x^{2}(ax^{2}-1){\biggl \{}{\frac {d^{2}\xi }{dx^{2}}}{\biggr \}}+{\biggl \{}x(1-ax^{2}){\biggl [}(1+c)+c{\biggr ]}-2[\gamma -(\alpha +\beta +1)ax^{2}]\cdot x{\biggr \}}{\frac {d\xi }{dx}}}
+
{
(
a
x
2
−
1
)
(
2
c
+
c
2
)
−
2
c
[
γ
−
(
α
+
β
+
1
)
a
x
2
]
+
4
a
α
β
x
2
}
ξ
{\displaystyle +{\biggl \{}(ax^{2}-1)(2c+c^{2})-2c[\gamma -(\alpha +\beta +1)ax^{2}]+4a\alpha \beta x^{2}{\biggr \}}\xi }
=
{\displaystyle =}
x
2
(
a
x
2
−
1
)
{
d
2
ξ
d
x
2
}
+
{
(
1
+
2
c
−
2
γ
)
−
a
x
2
(
1
+
2
c
)
+
2
(
α
+
β
+
1
)
a
x
2
]
}
x
⋅
d
ξ
d
x
{\displaystyle x^{2}(ax^{2}-1){\biggl \{}{\frac {d^{2}\xi }{dx^{2}}}{\biggr \}}+{\biggl \{}(1+2c-2\gamma )-ax^{2}(1+2c)+2(\alpha +\beta +1)ax^{2}]{\biggr \}}x\cdot {\frac {d\xi }{dx}}}
+
{
(
a
x
2
−
1
)
(
2
c
+
c
2
)
−
2
c
[
γ
−
(
α
+
β
+
1
)
a
x
2
]
+
4
a
α
β
x
2
}
ξ
{\displaystyle +{\biggl \{}(ax^{2}-1)(2c+c^{2})-2c[\gamma -(\alpha +\beta +1)ax^{2}]+4a\alpha \beta x^{2}{\biggr \}}\xi }
LAWE
Familiar Foundation
Drawing from an accompanying discussion , we have the,
LAWE: Linear Adiabatic Wave (or Radial Pulsation ) Equation
d
2
x
d
r
0
2
+
[
4
r
0
−
(
g
0
ρ
0
P
0
)
]
d
x
d
r
0
+
(
ρ
0
γ
g
P
0
)
[
ω
2
+
(
4
−
3
γ
g
)
g
0
r
0
]
x
=
0
{\displaystyle {\frac {d^{2}x}{dr_{0}^{2}}}+{\biggl [}{\frac {4}{r_{0}}}-{\biggl (}{\frac {g_{0}\rho _{0}}{P_{0}}}{\biggr )}{\biggr ]}{\frac {dx}{dr_{0}}}+{\biggl (}{\frac {\rho _{0}}{\gamma _{\mathrm {g} }P_{0}}}{\biggr )}{\biggl [}\omega ^{2}+(4-3\gamma _{\mathrm {g} }){\frac {g_{0}}{r_{0}}}{\biggr ]}x=0}
where,
g
0
{\displaystyle g_{0}}
=
{\displaystyle =}
−
1
ρ
0
d
P
0
d
r
0
.
{\displaystyle -{\frac {1}{\rho _{0}}}{\frac {dP_{0}}{dr_{0}}}\,.}
Multiplying through by
R
2
{\displaystyle R^{2}}
, and making the variable substitutions,
x
{\displaystyle x}
→
{\displaystyle \rightarrow }
f
,
{\displaystyle f\,,}
r
0
R
{\displaystyle {\frac {r_{0}}{R}}}
→
{\displaystyle \rightarrow }
x
,
{\displaystyle x\,,}
(
4
−
3
γ
g
)
{\displaystyle (4-3\gamma _{g})}
→
{\displaystyle \rightarrow }
−
α
γ
g
,
{\displaystyle -\alpha \gamma _{g}\,,}
the LAWE may be rewritten as,
0
{\displaystyle 0}
=
{\displaystyle =}
d
2
f
d
x
2
+
[
4
x
−
(
g
0
ρ
0
R
P
0
)
]
d
f
d
x
+
(
ρ
0
R
2
γ
g
P
0
)
[
ω
2
−
α
γ
g
g
0
r
0
]
f
{\displaystyle {\frac {d^{2}f}{dx^{2}}}+{\biggl [}{\frac {4}{x}}-{\biggl (}{\frac {g_{0}\rho _{0}R}{P_{0}}}{\biggr )}{\biggr ]}{\frac {df}{dx}}+{\biggl (}{\frac {\rho _{0}R^{2}}{\gamma _{\mathrm {g} }P_{0}}}{\biggr )}{\biggl [}\omega ^{2}-{\frac {\alpha \gamma _{g}g_{0}}{r_{0}}}{\biggr ]}f}
=
{\displaystyle =}
d
2
f
d
x
2
+
1
x
[
4
−
(
g
0
ρ
0
r
0
P
0
)
]
d
f
d
x
+
[
(
ω
2
ρ
0
R
2
γ
g
P
0
)
−
α
γ
g
g
0
r
0
(
ρ
0
R
2
γ
g
P
0
)
]
f
{\displaystyle {\frac {d^{2}f}{dx^{2}}}+{\frac {1}{x}}{\biggl [}4-{\biggl (}{\frac {g_{0}\rho _{0}r_{0}}{P_{0}}}{\biggr )}{\biggr ]}{\frac {df}{dx}}+{\biggl [}{\biggl (}{\frac {\omega ^{2}\rho _{0}R^{2}}{\gamma _{\mathrm {g} }P_{0}}}{\biggr )}-{\frac {\alpha \gamma _{g}g_{0}}{r_{0}}}{\biggl (}{\frac {\rho _{0}R^{2}}{\gamma _{\mathrm {g} }P_{0}}}{\biggr )}{\biggr ]}f}
=
{\displaystyle =}
d
2
f
d
x
2
+
1
x
[
4
−
(
g
0
ρ
0
r
0
P
0
)
]
d
f
d
x
+
[
(
ω
2
ρ
0
R
2
γ
g
P
0
)
−
α
x
2
(
g
0
r
0
ρ
0
P
0
)
]
f
.
{\displaystyle {\frac {d^{2}f}{dx^{2}}}+{\frac {1}{x}}{\biggl [}4-{\biggl (}{\frac {g_{0}\rho _{0}r_{0}}{P_{0}}}{\biggr )}{\biggr ]}{\frac {df}{dx}}+{\biggl [}{\biggl (}{\frac {\omega ^{2}\rho _{0}R^{2}}{\gamma _{\mathrm {g} }P_{0}}}{\biggr )}-{\frac {\alpha }{x^{2}}}{\biggl (}{\frac {g_{0}r_{0}\rho _{0}}{P_{0}}}{\biggr )}{\biggr ]}f\,.}
If we furthermore adopt the variable definition,
μ
{\displaystyle \mu }
≡
{\displaystyle \equiv }
(
g
0
ρ
0
r
0
P
0
)
=
−
d
ln
P
0
d
ln
r
0
,
{\displaystyle {\biggl (}{\frac {g_{0}\rho _{0}r_{0}}{P_{0}}}{\biggr )}=-{\frac {d\ln P_{0}}{d\ln r_{0}}}\,,}
we obtain what we will refer to as the,
Kopal (1948) LAWE
0
{\displaystyle 0}
=
{\displaystyle =}
d
2
f
d
x
2
+
(
4
−
μ
)
x
⋅
d
f
d
x
+
[
(
ω
2
ρ
0
R
2
γ
g
P
0
)
−
α
μ
x
2
]
f
.
{\displaystyle {\frac {d^{2}f}{dx^{2}}}+{\frac {(4-\mu )}{x}}\cdot {\frac {df}{dx}}+{\biggl [}{\biggl (}{\frac {\omega ^{2}\rho _{0}R^{2}}{\gamma _{\mathrm {g} }P_{0}}}{\biggr )}-{\frac {\alpha \mu }{x^{2}}}{\biggr ]}f\,.}
📚 Kopal (1948) , p. 378, Eq. (6)
📚 Van der Borght (1970) , p. 325, Eq. (1)
Specifically for Polytropes
Let's look at the expression for the function,
μ
{\displaystyle \mu }
, that arises in the context of polytropic spheres.
General Expression for the Function μ
First, we note that,
r
0
{\displaystyle r_{0}}
=
{\displaystyle =}
a
ξ
,
{\displaystyle a\xi \,,}
ρ
0
{\displaystyle \rho _{0}}
=
{\displaystyle =}
ρ
c
θ
n
,
{\displaystyle \rho _{c}\theta ^{n}\,,}
P
0
{\displaystyle P_{0}}
=
{\displaystyle =}
K
[
ρ
c
θ
n
]
(
n
+
1
)
/
n
,
{\displaystyle K{\biggl [}\rho _{c}\theta ^{n}{\biggr ]}^{(n+1)/n}\,,}
M
r
{\displaystyle M_{r}}
=
{\displaystyle =}
4
π
a
3
ρ
c
(
−
ξ
2
d
θ
d
ξ
)
,
{\displaystyle 4\pi a^{3}\rho _{c}{\biggl (}-\xi ^{2}{\frac {d\theta }{d\xi }}{\biggr )}\,,}
g
0
≡
G
M
r
r
0
2
{\displaystyle g_{0}\equiv {\frac {GM_{r}}{r_{0}^{2}}}}
=
{\displaystyle =}
4
π
G
a
ρ
c
(
−
d
θ
d
ξ
)
,
{\displaystyle 4\pi Ga\rho _{c}{\biggl (}-{\frac {d\theta }{d\xi }}{\biggr )}\,,}
where,
a
{\displaystyle a}
≡
{\displaystyle \equiv }
[
(
n
+
1
)
K
4
π
G
]
1
/
2
ρ
c
(
1
−
n
)
/
2
n
.
{\displaystyle {\biggl [}{\frac {(n+1)K}{4\pi G}}{\biggr ]}^{1/2}\rho _{c}^{(1-n)/2n}\,.}
⇒
K
{\displaystyle \Rightarrow ~~~K}
=
{\displaystyle =}
[
4
π
G
(
n
+
1
)
]
a
2
ρ
c
(
n
−
1
)
/
n
.
{\displaystyle {\biggl [}{\frac {4\pi G}{(n+1)}}{\biggr ]}a^{2}\rho _{c}^{(n-1)/n}\,.}
Hence,
μ
=
(
g
0
ρ
0
r
0
P
0
)
{\displaystyle \mu ={\biggl (}{\frac {g_{0}\rho _{0}r_{0}}{P_{0}}}{\biggr )}}
=
{\displaystyle =}
4
π
G
a
ρ
c
(
−
d
θ
d
ξ
)
ρ
c
θ
n
a
ξ
[
ρ
c
θ
n
]
−
(
n
+
1
)
/
n
[
(
n
+
1
)
4
π
G
]
a
−
2
ρ
c
−
(
n
−
1
)
/
n
{\displaystyle 4\pi Ga\rho _{c}{\biggl (}-{\frac {d\theta }{d\xi }}{\biggr )}\rho _{c}\theta ^{n}a\xi {\biggl [}\rho _{c}\theta ^{n}{\biggr ]}^{-(n+1)/n}{\biggl [}{\frac {(n+1)}{4\pi G}}{\biggr ]}a^{-2}\rho _{c}^{-(n-1)/n}}
=
{\displaystyle =}
(
n
+
1
)
(
−
d
θ
d
ξ
)
θ
n
ξ
θ
−
(
n
+
1
)
ρ
c
2
−
(
n
+
1
)
/
n
−
(
n
−
1
)
/
n
{\displaystyle (n+1){\biggl (}-{\frac {d\theta }{d\xi }}{\biggr )}\theta ^{n}\xi \theta ^{-(n+1)}\rho _{c}^{2-(n+1)/n-(n-1)/n}}
=
{\displaystyle =}
(
n
+
1
)
(
−
ξ
θ
⋅
d
θ
d
ξ
)
=
(
n
+
1
)
(
−
d
ln
θ
d
ln
ξ
)
.
{\displaystyle (n+1){\biggl (}-{\frac {\xi }{\theta }}\cdot {\frac {d\theta }{d\xi }}{\biggr )}=(n+1){\biggl (}-{\frac {d\ln \theta }{d\ln \xi }}{\biggr )}\,.}
Alternatively,
μ
=
−
d
ln
P
0
d
ln
r
0
{\displaystyle \mu =-{\frac {d\ln P_{0}}{d\ln r_{0}}}}
=
{\displaystyle =}
−
r
0
P
0
⋅
d
P
0
d
r
0
{\displaystyle -{\frac {r_{0}}{P_{0}}}\cdot {\frac {dP_{0}}{dr_{0}}}}
=
{\displaystyle =}
−
ξ
θ
−
(
n
+
1
)
⋅
d
d
ξ
[
θ
(
n
+
1
)
]
{\displaystyle -\xi \theta ^{-(n+1)}\cdot {\frac {d}{d\xi }}{\biggl [}\theta ^{(n+1)}{\biggr ]}}
=
{\displaystyle =}
−
(
n
+
1
)
(
ξ
θ
⋅
d
θ
d
ξ
)
=
(
n
+
1
)
(
−
d
ln
θ
d
ln
ξ
)
.
{\displaystyle -(n+1){\biggl (}{\frac {\xi }{\theta }}\cdot {\frac {d\theta }{d\xi }}{\biggr )}=(n+1){\biggl (}-{\frac {d\ln \theta }{d\ln \xi }}{\biggr )}\,.}
Yes!
Trial Displacement Function
Now, building on an accompanying discussion , let's guess,
f
t
r
i
a
l
{\displaystyle f_{\mathrm {trial} }}
=
{\displaystyle =}
3
(
n
−
1
)
2
n
[
1
+
(
n
−
3
n
−
1
)
(
1
ξ
θ
n
)
d
θ
d
ξ
]
{\displaystyle {\frac {3(n-1)}{2n}}{\biggl [}1+{\biggl (}{\frac {n-3}{n-1}}{\biggr )}{\biggl (}{\frac {1}{\xi \theta ^{n}}}{\biggr )}{\frac {d\theta }{d\xi }}{\biggr ]}}
⇒
[
2
n
3
(
n
−
1
)
]
f
t
r
i
a
l
{\displaystyle \Rightarrow ~~~{\biggl [}{\frac {2n}{3(n-1)}}{\biggr ]}f_{\mathrm {trial} }}
=
{\displaystyle =}
1
+
(
n
−
3
n
−
1
)
ξ
−
1
θ
−
n
[
θ
ξ
⋅
d
ln
θ
d
ln
ξ
]
{\displaystyle 1+{\biggl (}{\frac {n-3}{n-1}}{\biggr )}\xi ^{-1}\theta ^{-n}{\biggl [}{\frac {\theta }{\xi }}\cdot {\frac {d\ln \theta }{d\ln \xi }}{\biggr ]}}
=
{\displaystyle =}
1
+
(
3
−
n
n
−
1
)
ξ
−
2
θ
(
1
−
n
)
⋅
(
−
d
ln
θ
d
ln
ξ
)
{\displaystyle 1+{\biggl (}{\frac {3-n}{n-1}}{\biggr )}\xi ^{-2}\theta ^{(1-n)}\cdot {\biggl (}-{\frac {d\ln \theta }{d\ln \xi }}{\biggr )}}
=
{\displaystyle =}
1
+
[
3
−
n
(
n
+
1
)
(
n
−
1
)
]
ξ
−
2
θ
(
1
−
n
)
⋅
μ
{\displaystyle 1+{\biggl [}{\frac {3-n}{(n+1)(n-1)}}{\biggr ]}\xi ^{-2}\theta ^{(1-n)}\cdot \mu }
Flipping it around, we have alternatively,
μ
{\displaystyle \mu }
=
{\displaystyle =}
[
(
n
+
1
)
(
n
−
1
)
3
−
n
]
{
[
2
n
3
(
n
−
1
)
]
f
t
r
i
a
l
−
1
}
ξ
2
θ
(
n
−
1
)
{\displaystyle {\biggl [}{\frac {(n+1)(n-1)}{3-n}}{\biggr ]}{\biggl \{}{\biggl [}{\frac {2n}{3(n-1)}}{\biggr ]}f_{\mathrm {trial} }-1{\biggr \}}\xi ^{2}\theta ^{(n-1)}}
Plug into Kopal (1948) LAWE
Replace ftrial by μ
Plugging this trial function into the Kopal (1948) LAWE and recognizing that
x
=
ξ
/
ξ
1
{\displaystyle x=\xi /\xi _{1}}
, we find,
ξ
1
−
2
⋅
L
A
W
E
{\displaystyle \xi _{1}^{-2}\cdot ~\mathrm {LAWE} }
=
{\displaystyle =}
d
2
f
t
r
i
a
l
d
ξ
2
+
(
4
−
μ
)
ξ
⋅
d
f
t
r
i
a
l
d
ξ
+
[
(
ω
2
ρ
0
R
2
γ
g
P
0
ξ
1
2
)
−
α
μ
ξ
2
]
f
t
r
i
a
l
{\displaystyle {\frac {d^{2}f_{\mathrm {trial} }}{d\xi ^{2}}}+{\frac {(4-\mu )}{\xi }}\cdot {\frac {df_{\mathrm {trial} }}{d\xi }}+{\biggl [}{\biggl (}{\frac {\omega ^{2}\rho _{0}R^{2}}{\gamma _{\mathrm {g} }P_{0}\xi _{1}^{2}}}{\biggr )}-{\frac {\alpha \mu }{\xi ^{2}}}{\biggr ]}f_{\mathrm {trial} }}
⇒
[
2
n
3
(
n
−
1
)
]
ξ
1
−
2
⋅
L
A
W
E
{\displaystyle \Rightarrow ~~~{\biggl [}{\frac {2n}{3(n-1)}}{\biggr ]}\xi _{1}^{-2}\cdot ~\mathrm {LAWE} }
=
{\displaystyle =}
d
2
d
ξ
2
{
1
+
[
3
−
n
(
n
+
1
)
(
n
−
1
)
]
ξ
−
2
θ
(
1
−
n
)
⋅
μ
}
+
(
4
−
μ
)
ξ
⋅
d
d
ξ
{
1
+
[
3
−
n
(
n
+
1
)
(
n
−
1
)
]
ξ
−
2
θ
(
1
−
n
)
⋅
μ
}
{\displaystyle {\frac {d^{2}}{d\xi ^{2}}}{\biggl \{}1+{\biggl [}{\frac {3-n}{(n+1)(n-1)}}{\biggr ]}\xi ^{-2}\theta ^{(1-n)}\cdot \mu {\biggr \}}+{\frac {(4-\mu )}{\xi }}\cdot {\frac {d}{d\xi }}{\biggl \{}1+{\biggl [}{\frac {3-n}{(n+1)(n-1)}}{\biggr ]}\xi ^{-2}\theta ^{(1-n)}\cdot \mu {\biggr \}}}
+
[
(
ω
2
ρ
0
R
2
γ
g
P
0
ξ
1
2
)
−
α
μ
ξ
2
]
{
1
+
[
3
−
n
(
n
+
1
)
(
n
−
1
)
]
ξ
−
2
θ
(
1
−
n
)
⋅
μ
}
{\displaystyle +{\biggl [}{\biggl (}{\frac {\omega ^{2}\rho _{0}R^{2}}{\gamma _{\mathrm {g} }P_{0}\xi _{1}^{2}}}{\biggr )}-{\frac {\alpha \mu }{\xi ^{2}}}{\biggr ]}{\biggl \{}1+{\biggl [}{\frac {3-n}{(n+1)(n-1)}}{\biggr ]}\xi ^{-2}\theta ^{(1-n)}\cdot \mu {\biggr \}}}
⇒
[
2
n
(
n
+
1
)
3
]
ξ
1
−
2
⋅
L
A
W
E
{\displaystyle \Rightarrow ~~~{\biggl [}{\frac {2n(n+1)}{3}}{\biggr ]}\xi _{1}^{-2}\cdot ~\mathrm {LAWE} }
=
{\displaystyle =}
d
2
d
ξ
2
{
(
3
−
n
)
ξ
−
2
θ
(
1
−
n
)
⋅
μ
}
+
(
4
−
μ
)
ξ
⋅
d
d
ξ
{
(
3
−
n
)
ξ
−
2
θ
(
1
−
n
)
⋅
μ
}
{\displaystyle {\frac {d^{2}}{d\xi ^{2}}}{\biggl \{}(3-n)\xi ^{-2}\theta ^{(1-n)}\cdot \mu {\biggr \}}+{\frac {(4-\mu )}{\xi }}\cdot {\frac {d}{d\xi }}{\biggl \{}(3-n)\xi ^{-2}\theta ^{(1-n)}\cdot \mu {\biggr \}}}
+
[
(
ω
2
ρ
0
R
2
γ
g
P
0
ξ
1
2
)
−
α
μ
ξ
2
]
{
(
n
+
1
)
(
n
−
1
)
+
(
3
−
n
)
ξ
−
2
θ
(
1
−
n
)
⋅
μ
}
{\displaystyle +{\biggl [}{\biggl (}{\frac {\omega ^{2}\rho _{0}R^{2}}{\gamma _{\mathrm {g} }P_{0}\xi _{1}^{2}}}{\biggr )}-{\frac {\alpha \mu }{\xi ^{2}}}{\biggr ]}{\biggl \{}(n+1)(n-1)+(3-n)\xi ^{-2}\theta ^{(1-n)}\cdot \mu {\biggr \}}}
⇒
[
2
n
(
n
+
1
)
3
(
3
−
n
)
]
ξ
1
−
2
⋅
L
A
W
E
{\displaystyle \Rightarrow ~~~{\biggl [}{\frac {2n(n+1)}{3(3-n)}}{\biggr ]}\xi _{1}^{-2}\cdot ~\mathrm {LAWE} }
=
{\displaystyle =}
d
2
d
ξ
2
{
ξ
−
2
θ
(
1
−
n
)
⋅
μ
}
+
(
4
−
μ
)
ξ
⋅
d
d
ξ
{
ξ
−
2
θ
(
1
−
n
)
⋅
μ
}
+
[
(
ω
2
ρ
0
R
2
γ
g
P
0
ξ
1
2
)
−
α
μ
ξ
2
]
{
(
n
+
1
)
(
n
−
1
)
(
3
−
n
)
+
ξ
−
2
θ
(
1
−
n
)
⋅
μ
}
.
{\displaystyle {\frac {d^{2}}{d\xi ^{2}}}{\biggl \{}\xi ^{-2}\theta ^{(1-n)}\cdot \mu {\biggr \}}+{\frac {(4-\mu )}{\xi }}\cdot {\frac {d}{d\xi }}{\biggl \{}\xi ^{-2}\theta ^{(1-n)}\cdot \mu {\biggr \}}+{\biggl [}{\biggl (}{\frac {\omega ^{2}\rho _{0}R^{2}}{\gamma _{\mathrm {g} }P_{0}\xi _{1}^{2}}}{\biggr )}-{\frac {\alpha \mu }{\xi ^{2}}}{\biggr ]}{\biggl \{}{\frac {(n+1)(n-1)}{(3-n)}}+\xi ^{-2}\theta ^{(1-n)}\cdot \mu {\biggr \}}\,.}
Noting that,
R
/
ξ
1
=
a
{\displaystyle R/\xi _{1}=a}
and
ρ
0
P
0
{\displaystyle {\frac {\rho _{0}}{P_{0}}}}
=
{\displaystyle =}
ρ
c
θ
n
K
−
1
ρ
c
−
(
n
+
1
)
/
n
θ
−
(
n
+
1
)
=
K
−
1
ρ
c
−
1
/
n
θ
−
1
,
{\displaystyle \rho _{c}\theta ^{n}K^{-1}\rho _{c}^{-(n+1)/n}\theta ^{-(n+1)}=K^{-1}\rho _{c}^{-1/n}\theta ^{-1}\,,}
the frequency-squared term may be rewritten as,
ω
2
γ
g
(
ρ
0
a
2
P
0
)
{\displaystyle {\frac {\omega ^{2}}{\gamma _{g}}}{\biggl (}{\frac {\rho _{0}a^{2}}{P_{0}}}{\biggr )}}
=
{\displaystyle =}
(
ω
2
γ
g
)
K
−
1
ρ
c
−
1
/
n
θ
−
1
[
(
n
+
1
)
K
4
π
G
]
ρ
c
(
1
−
n
)
/
n
=
(
ω
2
γ
g
)
[
(
n
+
1
)
4
π
ρ
c
G
]
θ
−
1
.
{\displaystyle {\biggl (}{\frac {\omega ^{2}}{\gamma _{g}}}{\biggr )}K^{-1}\rho _{c}^{-1/n}\theta ^{-1}{\biggl [}{\frac {(n+1)K}{4\pi G}}{\biggr ]}\rho _{c}^{(1-n)/n}={\biggl (}{\frac {\omega ^{2}}{\gamma _{g}}}{\biggr )}{\biggl [}{\frac {(n+1)}{4\pi \rho _{c}G}}{\biggr ]}\theta ^{-1}\,.}
Replace μ by ftrial
Making instead the alternate substitution, namely,
μ
{\displaystyle \mu }
=
{\displaystyle =}
[
(
n
+
1
)
(
n
−
1
)
3
−
n
]
{
[
2
n
3
(
n
−
1
)
]
f
t
r
i
a
l
−
1
}
ξ
2
θ
(
n
−
1
)
{\displaystyle {\biggl [}{\frac {(n+1)(n-1)}{3-n}}{\biggr ]}{\biggl \{}{\biggl [}{\frac {2n}{3(n-1)}}{\biggr ]}f_{\mathrm {trial} }-1{\biggr \}}\xi ^{2}\theta ^{(n-1)}}
=
{\displaystyle =}
{
2
n
(
n
+
1
)
3
(
3
−
n
)
⋅
f
t
r
i
a
l
−
(
n
+
1
)
(
n
−
1
)
3
−
n
}
ξ
2
θ
(
n
−
1
)
{\displaystyle {\biggl \{}{\frac {2n(n+1)}{3(3-n)}}\cdot f_{\mathrm {trial} }-{\frac {(n+1)(n-1)}{3-n}}{\biggr \}}\xi ^{2}\theta ^{(n-1)}}
=
{\displaystyle =}
1
3
(
3
−
n
)
[
2
n
(
n
+
1
)
⏟
A
f
t
r
i
a
l
−
3
(
n
+
1
)
(
n
−
1
)
⏟
B
]
ξ
2
θ
(
n
−
1
)
{\displaystyle {\frac {1}{3(3-n)}}{\biggl [}\underbrace {2n(n+1)} _{A}f_{\mathrm {trial} }-\underbrace {3(n+1)(n-1)} _{B}{\biggr ]}\xi ^{2}\theta ^{(n-1)}}
we have,
ξ
1
−
2
⋅
L
A
W
E
{\displaystyle \xi _{1}^{-2}\cdot ~\mathrm {LAWE} }
=
{\displaystyle =}
d
2
f
t
r
i
a
l
d
ξ
2
+
{
4
ξ
}
d
f
t
r
i
a
l
d
ξ
−
{
μ
ξ
}
d
f
t
r
i
a
l
d
ξ
+
(
ω
2
γ
g
)
[
(
n
+
1
)
4
π
ρ
c
G
]
f
t
r
i
a
l
θ
−
(
α
ξ
2
)
μ
f
t
r
i
a
l
{\displaystyle {\frac {d^{2}f_{\mathrm {trial} }}{d\xi ^{2}}}+{\biggl \{}{\frac {4}{\xi }}{\biggr \}}{\frac {df_{\mathrm {trial} }}{d\xi }}-{\biggl \{}{\frac {\mu }{\xi }}{\biggr \}}{\frac {df_{\mathrm {trial} }}{d\xi }}+{\biggl (}{\frac {\omega ^{2}}{\gamma _{g}}}{\biggr )}{\biggl [}{\frac {(n+1)}{4\pi \rho _{c}G}}{\biggr ]}{\frac {f_{\mathrm {trial} }}{\theta }}-{\biggl (}{\frac {\alpha }{\xi ^{2}}}{\biggr )}\mu f_{\mathrm {trial} }}
⇒
[
3
(
3
−
n
)
ξ
1
2
]
L
A
W
E
{\displaystyle \Rightarrow ~~~{\biggl [}{\frac {3(3-n)}{\xi _{1}^{2}}}{\biggr ]}~\mathrm {LAWE} }
=
{\displaystyle =}
3
(
3
−
n
)
d
2
f
t
r
i
a
l
d
ξ
2
+
{
12
(
3
−
n
)
ξ
}
d
f
t
r
i
a
l
d
ξ
−
[
A
f
t
r
i
a
l
−
B
]
ξ
θ
(
n
−
1
)
d
f
t
r
i
a
l
d
ξ
{\displaystyle 3(3-n){\frac {d^{2}f_{\mathrm {trial} }}{d\xi ^{2}}}+{\biggl \{}{\frac {12(3-n)}{\xi }}{\biggr \}}{\frac {df_{\mathrm {trial} }}{d\xi }}-{\biggl [}Af_{\mathrm {trial} }-B{\biggr ]}\xi \theta ^{(n-1)}{\frac {df_{\mathrm {trial} }}{d\xi }}}
+
3
(
3
−
n
)
(
ω
2
γ
g
)
[
(
n
+
1
)
4
π
ρ
c
G
]
f
t
r
i
a
l
θ
−
α
[
A
f
t
r
i
a
l
−
B
]
θ
(
n
−
1
)
f
t
r
i
a
l
.
{\displaystyle +3(3-n){\biggl (}{\frac {\omega ^{2}}{\gamma _{g}}}{\biggr )}{\biggl [}{\frac {(n+1)}{4\pi \rho _{c}G}}{\biggr ]}{\frac {f_{\mathrm {trial} }}{\theta }}-\alpha {\biggl [}Af_{\mathrm {trial} }-B{\biggr ]}\theta ^{(n-1)}f_{\mathrm {trial} }\,.}
Noting that,
d
d
ξ
[
ξ
θ
(
n
−
1
)
f
t
r
i
a
l
]
{\displaystyle {\frac {d}{d\xi }}{\biggl [}\xi \theta ^{(n-1)}f_{\mathrm {trial} }{\biggr ]}}
=
{\displaystyle =}
{
θ
(
n
−
1
)
f
t
r
i
a
l
+
ξ
f
t
r
i
a
l
(
n
−
1
)
θ
(
n
−
2
)
d
θ
d
ξ
+
ξ
θ
(
n
−
1
)
d
f
t
r
i
a
l
d
ξ
}
{\displaystyle {\biggl \{}\theta ^{(n-1)}f_{\mathrm {trial} }+\xi f_{\mathrm {trial} }(n-1)\theta ^{(n-2)}{\frac {d\theta }{d\xi }}+\xi \theta ^{(n-1)}{\frac {df_{\mathrm {trial} }}{d\xi }}{\biggr \}}}
⇒
ξ
θ
(
n
−
1
)
d
f
t
r
i
a
l
d
ξ
{\displaystyle \Rightarrow ~~~\xi \theta ^{(n-1)}{\frac {df_{\mathrm {trial} }}{d\xi }}}
=
{\displaystyle =}
d
d
ξ
[
ξ
θ
(
n
−
1
)
f
t
r
i
a
l
]
−
[
θ
(
n
−
1
)
f
t
r
i
a
l
+
ξ
f
t
r
i
a
l
(
n
−
1
)
θ
(
n
−
2
)
d
θ
d
ξ
]
,
{\displaystyle {\frac {d}{d\xi }}{\biggl [}\xi \theta ^{(n-1)}f_{\mathrm {trial} }{\biggr ]}-{\biggl [}\theta ^{(n-1)}f_{\mathrm {trial} }+\xi f_{\mathrm {trial} }(n-1)\theta ^{(n-2)}{\frac {d\theta }{d\xi }}{\biggr ]}\,,}
we furthermore can write,
⇒
[
3
(
3
−
n
)
ξ
1
2
]
L
A
W
E
{\displaystyle \Rightarrow ~~~{\biggl [}{\frac {3(3-n)}{\xi _{1}^{2}}}{\biggr ]}~\mathrm {LAWE} }
=
{\displaystyle =}
3
(
3
−
n
)
d
2
f
t
r
i
a
l
d
ξ
2
+
{
12
(
3
−
n
)
ξ
}
d
f
t
r
i
a
l
d
ξ
+
3
(
3
−
n
)
(
ω
2
γ
g
)
[
(
n
+
1
)
4
π
ρ
c
G
]
f
t
r
i
a
l
θ
{\displaystyle 3(3-n){\frac {d^{2}f_{\mathrm {trial} }}{d\xi ^{2}}}+{\biggl \{}{\frac {12(3-n)}{\xi }}{\biggr \}}{\frac {df_{\mathrm {trial} }}{d\xi }}+3(3-n){\biggl (}{\frac {\omega ^{2}}{\gamma _{g}}}{\biggr )}{\biggl [}{\frac {(n+1)}{4\pi \rho _{c}G}}{\biggr ]}{\frac {f_{\mathrm {trial} }}{\theta }}}
−
α
[
A
f
t
r
i
a
l
−
B
]
θ
(
n
−
1
)
f
t
r
i
a
l
−
[
A
f
t
r
i
a
l
−
B
]
{
d
d
ξ
[
ξ
θ
(
n
−
1
)
f
t
r
i
a
l
]
−
[
θ
(
n
−
1
)
f
t
r
i
a
l
+
ξ
f
t
r
i
a
l
(
n
−
1
)
θ
(
n
−
2
)
d
θ
d
ξ
]
}
{\displaystyle -\alpha {\biggl [}Af_{\mathrm {trial} }-B{\biggr ]}\theta ^{(n-1)}f_{\mathrm {trial} }-{\biggl [}Af_{\mathrm {trial} }-B{\biggr ]}{\biggl \{}{\frac {d}{d\xi }}{\biggl [}\xi \theta ^{(n-1)}f_{\mathrm {trial} }{\biggr ]}-{\biggl [}\theta ^{(n-1)}f_{\mathrm {trial} }+\xi f_{\mathrm {trial} }(n-1)\theta ^{(n-2)}{\frac {d\theta }{d\xi }}{\biggr ]}{\biggr \}}}
=
{\displaystyle =}
3
(
3
−
n
)
d
2
f
t
r
i
a
l
d
ξ
2
+
{
12
(
3
−
n
)
ξ
}
d
f
t
r
i
a
l
d
ξ
+
3
(
3
−
n
)
(
ω
2
γ
g
)
[
(
n
+
1
)
4
π
ρ
c
G
]
f
t
r
i
a
l
θ
{\displaystyle 3(3-n){\frac {d^{2}f_{\mathrm {trial} }}{d\xi ^{2}}}+{\biggl \{}{\frac {12(3-n)}{\xi }}{\biggr \}}{\frac {df_{\mathrm {trial} }}{d\xi }}+3(3-n){\biggl (}{\frac {\omega ^{2}}{\gamma _{g}}}{\biggr )}{\biggl [}{\frac {(n+1)}{4\pi \rho _{c}G}}{\biggr ]}{\frac {f_{\mathrm {trial} }}{\theta }}}
−
[
A
f
t
r
i
a
l
−
B
]
{
d
d
ξ
[
ξ
θ
(
n
−
1
)
f
t
r
i
a
l
]
−
[
θ
(
n
−
1
)
f
t
r
i
a
l
+
ξ
f
t
r
i
a
l
(
n
−
1
)
θ
(
n
−
2
)
d
θ
d
ξ
]
+
α
θ
(
n
−
1
)
f
t
r
i
a
l
}
{\displaystyle -{\biggl [}Af_{\mathrm {trial} }-B{\biggr ]}{\biggl \{}{\frac {d}{d\xi }}{\biggl [}\xi \theta ^{(n-1)}f_{\mathrm {trial} }{\biggr ]}-{\biggl [}\theta ^{(n-1)}f_{\mathrm {trial} }+\xi f_{\mathrm {trial} }(n-1)\theta ^{(n-2)}{\frac {d\theta }{d\xi }}{\biggr ]}+\alpha \theta ^{(n-1)}f_{\mathrm {trial} }{\biggr \}}}
Seek Hypergeometric Form
Start with the standard LAWE, namely,
0
{\displaystyle 0}
=
{\displaystyle =}
d
2
f
d
ξ
2
+
(
4
−
μ
)
ξ
⋅
d
f
d
ξ
+
[
(
ω
2
ρ
0
R
2
γ
g
P
0
)
−
α
μ
ξ
2
]
f
{\displaystyle {\frac {d^{2}f}{d\xi ^{2}}}+{\frac {(4-\mu )}{\xi }}\cdot {\frac {df}{d\xi }}+{\biggl [}{\biggl (}{\frac {\omega ^{2}\rho _{0}R^{2}}{\gamma _{\mathrm {g} }P_{0}}}{\biggr )}-{\frac {\alpha \mu }{\xi ^{2}}}{\biggr ]}f}
=
{\displaystyle =}
d
2
f
d
ξ
2
+
1
ξ
[
4
+
(
n
+
1
)
ξ
θ
′
θ
]
⋅
d
f
d
ξ
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
1
θ
+
α
θ
′
ξ
θ
]
f
.
{\displaystyle {\frac {d^{2}f}{d\xi ^{2}}}+{\frac {1}{\xi }}{\biggl [}4+(n+1){\frac {\xi \theta ^{'}}{\theta }}{\biggr ]}\cdot {\frac {df}{d\xi }}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}{\frac {1}{\theta }}+{\frac {\alpha \theta ^{'}}{\xi \theta }}{\biggr ]}f\,.}
Part I
Try switching the independent variable from
ξ
{\displaystyle \xi }
to
z
{\displaystyle z}
such that,
z
{\displaystyle z}
=
{\displaystyle =}
ξ
−
1
θ
−
n
(
θ
′
)
⇒
(
θ
′
)
=
ξ
θ
n
z
{\displaystyle \xi ^{-1}\theta ^{-n}(\theta ^{'})~~~\Rightarrow ~~~(\theta ^{'})=\xi \theta ^{n}z}
⇒
d
z
d
ξ
{\displaystyle \Rightarrow ~~~{\frac {dz}{d\xi }}}
=
{\displaystyle =}
−
ξ
−
2
θ
−
n
(
θ
′
)
−
n
ξ
−
1
θ
−
n
−
1
(
θ
′
)
2
+
ξ
−
1
θ
−
n
(
θ
″
)
{\displaystyle -\xi ^{-2}\theta ^{-n}(\theta ^{'})-n\xi ^{-1}\theta ^{-n-1}(\theta ^{'})^{2}+\xi ^{-1}\theta ^{-n}(\theta ^{''})}
=
{\displaystyle =}
−
[
ξ
−
2
θ
−
n
(
θ
′
)
+
n
ξ
−
1
θ
−
n
−
1
(
θ
′
)
2
+
ξ
−
1
θ
−
n
(
θ
n
+
2
θ
′
ξ
)
]
{\displaystyle -{\biggl [}\xi ^{-2}\theta ^{-n}(\theta ^{'})+n\xi ^{-1}\theta ^{-n-1}(\theta ^{'})^{2}+\xi ^{-1}\theta ^{-n}{\biggl (}\theta ^{n}+{\frac {2\theta '}{\xi }}{\biggr )}{\biggr ]}}
=
{\displaystyle =}
−
[
ξ
−
2
θ
−
n
(
ξ
θ
n
z
)
+
n
ξ
−
1
θ
−
n
−
1
(
ξ
θ
n
z
)
2
+
ξ
−
1
θ
−
n
(
θ
n
)
+
ξ
−
1
θ
−
n
(
2
θ
′
ξ
)
]
{\displaystyle -{\biggl [}\xi ^{-2}\theta ^{-n}(\xi \theta ^{n}z)+n\xi ^{-1}\theta ^{-n-1}(\xi \theta ^{n}z)^{2}+\xi ^{-1}\theta ^{-n}{\biggl (}\theta ^{n}{\biggr )}+\xi ^{-1}\theta ^{-n}{\biggl (}{\frac {2\theta '}{\xi }}{\biggr )}{\biggr ]}}
=
{\displaystyle =}
−
[
ξ
−
1
z
+
n
ξ
θ
n
−
1
z
2
+
ξ
−
1
+
2
ξ
−
1
z
]
{\displaystyle -{\biggl [}\xi ^{-1}z+n\xi \theta ^{n-1}z^{2}+\xi ^{-1}+2\xi ^{-1}z{\biggr ]}}
=
{\displaystyle =}
−
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
;
{\displaystyle -{\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}\,;}
and,
d
2
z
d
ξ
2
{\displaystyle {\frac {d^{2}z}{d\xi ^{2}}}}
=
{\displaystyle =}
−
{
−
ξ
−
2
(
1
+
3
z
)
+
3
ξ
−
1
⋅
d
z
d
ξ
+
n
θ
n
−
1
z
2
+
n
(
n
−
1
)
ξ
θ
n
−
2
(
θ
′
)
z
2
+
2
n
ξ
θ
n
−
1
z
⋅
d
z
d
ξ
}
{\displaystyle -{\biggl \{}-\xi ^{-2}(1+3z)+3\xi ^{-1}\cdot {\frac {dz}{d\xi }}+n\theta ^{n-1}z^{2}+n(n-1)\xi \theta ^{n-2}(\theta ^{'})z^{2}+2n\xi \theta ^{n-1}z\cdot {\frac {dz}{d\xi }}{\biggr \}}}
=
{\displaystyle =}
−
{
−
ξ
−
2
(
1
+
3
z
)
−
3
ξ
−
1
⋅
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
+
n
θ
n
−
1
z
2
+
n
(
n
−
1
)
ξ
θ
n
−
2
(
ξ
θ
n
z
)
z
2
−
2
n
ξ
θ
n
−
1
z
⋅
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
}
{\displaystyle -{\biggl \{}-\xi ^{-2}(1+3z)-3\xi ^{-1}\cdot {\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}+n\theta ^{n-1}z^{2}+n(n-1)\xi \theta ^{n-2}(\xi \theta ^{n}z)z^{2}-2n\xi \theta ^{n-1}z\cdot {\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}{\biggr \}}}
=
{\displaystyle =}
{
ξ
−
2
(
1
+
3
z
)
+
3
ξ
−
1
⋅
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
−
n
θ
n
−
1
z
2
−
n
(
n
−
1
)
ξ
θ
n
−
2
(
ξ
θ
n
z
)
z
2
+
2
n
ξ
θ
n
−
1
z
⋅
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
}
{\displaystyle {\biggl \{}\xi ^{-2}(1+3z)+3\xi ^{-1}\cdot {\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}-n\theta ^{n-1}z^{2}-n(n-1)\xi \theta ^{n-2}(\xi \theta ^{n}z)z^{2}+2n\xi \theta ^{n-1}z\cdot {\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}{\biggr \}}}
=
{\displaystyle =}
{
ξ
−
2
(
1
+
3
z
)
+
[
3
ξ
−
2
(
1
+
3
z
)
+
3
n
θ
n
−
1
z
2
]
−
n
θ
n
−
1
z
2
−
n
(
n
−
1
)
ξ
2
θ
2
(
n
−
1
)
z
3
+
[
2
n
θ
n
−
1
z
(
1
+
3
z
)
+
2
n
2
ξ
2
θ
2
(
n
−
1
)
z
3
]
}
{\displaystyle {\biggl \{}\xi ^{-2}(1+3z)+{\biggl [}3\xi ^{-2}(1+3z)+3n\theta ^{n-1}z^{2}{\biggr ]}-n\theta ^{n-1}z^{2}-n(n-1)\xi ^{2}\theta ^{2(n-1)}z^{3}+{\biggl [}2n\theta ^{n-1}z(1+3z)+2n^{2}\xi ^{2}\theta ^{2(n-1)}z^{3}{\biggr ]}{\biggr \}}}
=
{\displaystyle =}
4
ξ
−
2
(
1
+
3
z
)
+
2
n
θ
n
−
1
[
z
(
1
+
3
z
)
+
z
2
]
+
[
2
n
2
−
n
(
n
−
1
)
]
ξ
2
θ
2
(
n
−
1
)
z
3
{\displaystyle 4\xi ^{-2}(1+3z)+2n\theta ^{n-1}[z(1+3z)+z^{2}]+[2n^{2}-n(n-1)]\xi ^{2}\theta ^{2(n-1)}z^{3}}
=
{\displaystyle =}
4
ξ
−
2
(
1
+
3
z
)
+
2
n
θ
n
−
1
[
1
+
4
z
]
z
+
n
(
n
+
1
)
ξ
2
θ
2
(
n
−
1
)
z
3
.
{\displaystyle 4\xi ^{-2}(1+3z)+2n\theta ^{n-1}[1+4z]z+n(n+1)\xi ^{2}\theta ^{2(n-1)}z^{3}\,.}
Part II
Part I Summary …
(
θ
′
)
{\displaystyle (\theta ^{'})}
=
{\displaystyle =}
ξ
θ
n
z
,
{\displaystyle \xi \theta ^{n}z\,,}
d
z
d
ξ
{\displaystyle {\frac {dz}{d\xi }}}
=
{\displaystyle =}
−
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
,
{\displaystyle -{\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}\,,}
d
2
z
d
ξ
2
{\displaystyle {\frac {d^{2}z}{d\xi ^{2}}}}
=
{\displaystyle =}
4
ξ
−
2
(
1
+
3
z
)
+
2
n
θ
n
−
1
[
1
+
4
z
]
z
+
n
(
n
+
1
)
ξ
2
θ
2
(
n
−
1
)
z
3
.
{\displaystyle 4\xi ^{-2}(1+3z)+2n\theta ^{n-1}[1+4z]z+n(n+1)\xi ^{2}\theta ^{2(n-1)}z^{3}\,.}
Also,
d
f
d
ξ
{\displaystyle {\frac {df}{d\xi }}}
→
{\displaystyle \rightarrow }
d
z
d
ξ
⋅
d
f
d
z
;
{\displaystyle {\frac {dz}{d\xi }}\cdot {\frac {df}{dz}}\,;}
d
2
f
d
ξ
2
=
d
d
ξ
[
d
z
d
ξ
⋅
d
f
d
z
]
{\displaystyle {\frac {d^{2}f}{d\xi ^{2}}}={\frac {d}{d\xi }}{\biggl [}{\frac {dz}{d\xi }}\cdot {\frac {df}{dz}}{\biggr ]}}
→
{\displaystyle \rightarrow }
d
2
z
d
ξ
2
⋅
d
f
d
z
+
[
d
z
d
ξ
]
2
⋅
d
2
f
d
z
2
{\displaystyle {\frac {d^{2}z}{d\xi ^{2}}}\cdot {\frac {df}{dz}}+{\biggl [}{\frac {dz}{d\xi }}{\biggr ]}^{2}\cdot {\frac {d^{2}f}{dz^{2}}}}
As a result,
LAWE
=
{\displaystyle =}
d
2
f
d
ξ
2
+
1
ξ
[
4
+
(
n
+
1
)
ξ
θ
′
θ
]
⋅
d
f
d
ξ
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
1
θ
+
α
θ
′
ξ
θ
]
f
{\displaystyle {\frac {d^{2}f}{d\xi ^{2}}}+{\frac {1}{\xi }}{\biggl [}4+(n+1){\frac {\xi \theta ^{'}}{\theta }}{\biggr ]}\cdot {\frac {df}{d\xi }}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}{\frac {1}{\theta }}+{\frac {\alpha \theta ^{'}}{\xi \theta }}{\biggr ]}f}
=
{\displaystyle =}
d
2
f
d
ξ
2
{\displaystyle {\frac {d^{2}f}{d\xi ^{2}}}}
+
1
ξ
[
4
+
(
n
+
1
)
ξ
2
θ
n
−
1
z
]
⋅
d
f
d
ξ
{\displaystyle +{\frac {1}{\xi }}{\biggl [}4+(n+1)\xi ^{2}\theta ^{n-1}z{\biggr ]}\cdot {\frac {df}{d\xi }}}
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
1
θ
+
α
θ
n
−
1
z
]
f
{\displaystyle +(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}{\frac {1}{\theta }}+\alpha \theta ^{n-1}z{\biggr ]}f}
=
{\displaystyle =}
[
d
z
d
ξ
]
2
⋅
d
2
f
d
z
2
+
d
2
z
d
ξ
2
⋅
d
f
d
z
{\displaystyle {\biggl [}{\frac {dz}{d\xi }}{\biggr ]}^{2}\cdot {\frac {d^{2}f}{dz^{2}}}+{\frac {d^{2}z}{d\xi ^{2}}}\cdot {\frac {df}{dz}}}
+
[
4
ξ
−
1
+
(
n
+
1
)
ξ
θ
n
−
1
z
]
d
z
d
ξ
⋅
d
f
d
z
{\displaystyle +{\biggl [}4\xi ^{-1}+(n+1)\xi \theta ^{n-1}z{\biggr ]}{\frac {dz}{d\xi }}\cdot {\frac {df}{dz}}}
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
+
α
θ
n
−
1
z
]
f
{\displaystyle +(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}+\alpha \theta ^{n-1}z{\biggr ]}f}
=
{\displaystyle =}
{
d
z
d
ξ
}
2
⋅
d
2
f
d
z
2
{\displaystyle {\biggl \{}{\frac {dz}{d\xi }}{\biggr \}}^{2}\cdot {\frac {d^{2}f}{dz^{2}}}}
+
{
d
2
z
d
ξ
2
+
[
4
ξ
−
1
+
(
n
+
1
)
ξ
θ
n
−
1
z
]
d
z
d
ξ
}
d
f
d
z
{\displaystyle +{\biggl \{}{\frac {d^{2}z}{d\xi ^{2}}}+{\biggl [}4\xi ^{-1}+(n+1)\xi \theta ^{n-1}z{\biggr ]}{\frac {dz}{d\xi }}{\biggr \}}{\frac {df}{dz}}}
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
+
α
θ
n
−
1
z
]
f
{\displaystyle +(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}+\alpha \theta ^{n-1}z{\biggr ]}f}
=
{\displaystyle =}
{
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
}
2
⋅
d
2
f
d
z
2
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
+
α
θ
n
−
1
z
]
f
{\displaystyle {\biggl \{}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr \}}^{2}\cdot {\frac {d^{2}f}{dz^{2}}}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}+\alpha \theta ^{n-1}z{\biggr ]}f}
+
{
[
4
ξ
−
2
(
1
+
3
z
)
+
2
n
θ
n
−
1
[
1
+
4
z
]
z
+
n
(
n
+
1
)
ξ
2
θ
2
(
n
−
1
)
z
3
]
−
[
4
ξ
−
1
+
(
n
+
1
)
ξ
θ
n
−
1
z
]
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
}
d
f
d
z
{\displaystyle +{\biggl \{}{\biggl [}4\xi ^{-2}(1+3z)+2n\theta ^{n-1}[1+4z]z+n(n+1)\xi ^{2}\theta ^{2(n-1)}z^{3}{\biggr ]}-{\biggl [}4\xi ^{-1}+(n+1)\xi \theta ^{n-1}z{\biggr ]}{\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}{\biggr \}}{\frac {df}{dz}}}
=
{\displaystyle =}
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
⋅
d
2
f
d
z
2
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
+
α
θ
n
−
1
z
]
f
{\displaystyle {\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}{\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}\cdot {\frac {d^{2}f}{dz^{2}}}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}+\alpha \theta ^{n-1}z{\biggr ]}f}
+
[
4
ξ
−
2
(
1
+
3
z
)
+
2
n
θ
n
−
1
[
1
+
4
z
]
z
+
n
(
n
+
1
)
ξ
2
θ
2
(
n
−
1
)
z
3
]
d
f
d
z
{\displaystyle +{\biggl [}4\xi ^{-2}(1+3z)+2n\theta ^{n-1}[1+4z]z+n(n+1)\xi ^{2}\theta ^{2(n-1)}z^{3}{\biggr ]}{\frac {df}{dz}}}
−
[
4
ξ
−
1
]
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
d
f
d
z
−
[
(
n
+
1
)
ξ
θ
n
−
1
z
]
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
d
f
d
z
{\displaystyle -{\biggl [}4\xi ^{-1}{\biggr ]}{\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}{\frac {df}{dz}}-{\biggl [}(n+1)\xi \theta ^{n-1}z{\biggr ]}{\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}{\frac {df}{dz}}}
=
{\displaystyle =}
[
ξ
−
2
(
1
+
3
z
)
2
+
2
n
(
1
+
3
z
)
z
2
θ
n
−
1
+
n
2
ξ
2
θ
2
(
n
−
1
)
z
4
]
⋅
d
2
f
d
z
2
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
+
α
θ
n
−
1
z
]
f
{\displaystyle {\biggl [}\xi ^{-2}(1+3z)^{2}+2n(1+3z)z^{2}\theta ^{n-1}+n^{2}\xi ^{2}\theta ^{2(n-1)}z^{4}{\biggr ]}\cdot {\frac {d^{2}f}{dz^{2}}}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}+\alpha \theta ^{n-1}z{\biggr ]}f}
+
[
4
ξ
−
2
(
1
+
3
z
)
+
2
n
θ
n
−
1
[
1
+
4
z
]
z
+
n
(
n
+
1
)
ξ
2
θ
2
(
n
−
1
)
z
3
]
d
f
d
z
{\displaystyle +{\biggl [}4\xi ^{-2}(1+3z)+2n\theta ^{n-1}[1+4z]z+n(n+1)\xi ^{2}\theta ^{2(n-1)}z^{3}{\biggr ]}{\frac {df}{dz}}}
−
[
4
ξ
−
2
(
1
+
3
z
)
+
4
n
θ
n
−
1
z
2
]
d
f
d
z
−
[
(
n
+
1
)
θ
n
−
1
(
1
+
3
z
)
z
+
n
(
n
+
1
)
ξ
2
θ
2
(
n
−
1
)
z
3
]
d
f
d
z
{\displaystyle -{\biggl [}4\xi ^{-2}(1+3z)+4n\theta ^{n-1}z^{2}{\biggr ]}{\frac {df}{dz}}-{\biggl [}(n+1)\theta ^{n-1}(1+3z)z+n(n+1)\xi ^{2}\theta ^{2(n-1)}z^{3}{\biggr ]}{\frac {df}{dz}}}
=
{\displaystyle =}
[
ξ
−
2
(
1
+
3
z
)
2
+
2
n
(
1
+
3
z
)
z
2
θ
n
−
1
+
n
2
ξ
2
θ
2
(
n
−
1
)
z
4
]
⋅
d
2
f
d
z
2
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
+
α
θ
n
−
1
z
]
f
{\displaystyle {\biggl [}\xi ^{-2}(1+3z)^{2}+2n(1+3z)z^{2}\theta ^{n-1}+n^{2}\xi ^{2}\theta ^{2(n-1)}z^{4}{\biggr ]}\cdot {\frac {d^{2}f}{dz^{2}}}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}+\alpha \theta ^{n-1}z{\biggr ]}f}
+
{
[
2
n
θ
n
−
1
(
1
+
4
z
)
]
−
[
4
n
θ
n
−
1
z
]
−
[
(
n
+
1
)
θ
n
−
1
(
1
+
3
z
)
]
}
z
⋅
d
f
d
z
{\displaystyle +{\biggl \{}{\biggl [}2n\theta ^{n-1}(1+4z){\biggr ]}-{\biggl [}4n\theta ^{n-1}z{\biggr ]}-{\biggl [}(n+1)\theta ^{n-1}(1+3z){\biggr ]}{\biggr \}}z\cdot {\frac {df}{dz}}}
=
{\displaystyle =}
[
ξ
−
2
(
1
+
3
z
)
2
+
2
n
(
1
+
3
z
)
z
2
θ
n
−
1
+
n
2
ξ
2
θ
2
(
n
−
1
)
z
4
]
⋅
d
2
f
d
z
2
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
+
α
θ
n
−
1
z
]
f
{\displaystyle {\biggl [}\xi ^{-2}(1+3z)^{2}+2n(1+3z)z^{2}\theta ^{n-1}+n^{2}\xi ^{2}\theta ^{2(n-1)}z^{4}{\biggr ]}\cdot {\frac {d^{2}f}{dz^{2}}}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}+\alpha \theta ^{n-1}z{\biggr ]}f}
+
{
[
2
n
+
8
n
z
]
−
[
4
n
z
]
−
[
(
n
+
1
)
+
3
(
n
+
1
)
z
]
}
θ
n
−
1
z
⋅
d
f
d
z
{\displaystyle +{\biggl \{}{\biggl [}2n+8nz{\biggr ]}-{\biggl [}4nz{\biggr ]}-{\biggl [}(n+1)+3(n+1)z{\biggr ]}{\biggr \}}\theta ^{n-1}z\cdot {\frac {df}{dz}}}
=
{\displaystyle =}
[
ξ
−
2
(
1
+
3
z
)
2
+
2
n
(
1
+
3
z
)
z
2
θ
n
−
1
+
n
2
ξ
2
θ
2
(
n
−
1
)
z
4
]
⋅
d
2
f
d
z
2
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
+
α
θ
n
−
1
z
]
f
{\displaystyle {\biggl [}\xi ^{-2}(1+3z)^{2}+2n(1+3z)z^{2}\theta ^{n-1}+n^{2}\xi ^{2}\theta ^{2(n-1)}z^{4}{\biggr ]}\cdot {\frac {d^{2}f}{dz^{2}}}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}+\alpha \theta ^{n-1}z{\biggr ]}f}
+
{
(
n
−
1
)
+
(
n
−
3
)
z
}
θ
n
−
1
z
⋅
d
f
d
z
.
{\displaystyle +{\biggl \{}(n-1)+(n-3)z{\biggr \}}\theta ^{n-1}z\cdot {\frac {df}{dz}}\,.}
Part III
Now, suppose that
f
=
(
a
0
+
b
0
z
)
{\displaystyle f=(a_{0}+b_{0}z)}
. We have,
LAWE
=
{\displaystyle =}
[
ξ
−
2
(
1
+
3
z
)
2
+
2
n
(
1
+
3
z
)
z
2
θ
n
−
1
+
n
2
ξ
2
θ
2
(
n
−
1
)
z
4
]
⋅
d
2
f
d
z
2
0
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
+
α
θ
n
−
1
z
]
(
a
0
+
b
0
z
)
+
{
(
n
−
1
)
+
(
n
−
3
)
z
}
θ
n
−
1
z
⋅
b
0
{\displaystyle {\biggl [}\xi ^{-2}(1+3z)^{2}+2n(1+3z)z^{2}\theta ^{n-1}+n^{2}\xi ^{2}\theta ^{2(n-1)}z^{4}{\biggr ]}\cdot {\cancelto {0}{\frac {d^{2}f}{dz^{2}}}}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}+\alpha \theta ^{n-1}z{\biggr ]}(a_{0}+b_{0}z)+{\biggl \{}(n-1)+(n-3)z{\biggr \}}\theta ^{n-1}z\cdot b_{0}}
=
{\displaystyle =}
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
]
(
a
0
+
b
0
z
)
+
{
(
n
+
1
)
α
(
a
0
+
b
0
z
)
+
b
0
(
n
−
1
)
+
b
0
(
n
−
3
)
z
}
θ
n
−
1
z
{\displaystyle (n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}{\biggr ]}(a_{0}+b_{0}z)+{\biggl \{}(n+1)\alpha (a_{0}+b_{0}z)+b_{0}(n-1)+b_{0}(n-3)z{\biggr \}}\theta ^{n-1}z}
=
{\displaystyle =}
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
]
(
a
0
+
b
0
z
)
+
{
(
n
+
1
)
α
(
a
0
)
+
b
0
(
n
−
1
)
+
[
(
n
−
3
)
+
(
n
+
1
)
α
]
b
0
z
}
θ
n
−
1
z
.
{\displaystyle (n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}{\biggr ]}(a_{0}+b_{0}z)+{\biggl \{}(n+1)\alpha (a_{0})+b_{0}(n-1)+{\biggl [}(n-3)+(n+1)\alpha {\biggr ]}b_{0}z{\biggr \}}\theta ^{n-1}z\,.}
Now, in order for the last term to be zero, we need,
0
{\displaystyle 0}
=
{\displaystyle =}
[
(
n
−
3
)
+
(
n
+
1
)
α
]
{\displaystyle {\biggl [}(n-3)+(n+1)\alpha {\biggr ]}}
⇒
α
{\displaystyle \Rightarrow ~~~\alpha }
=
{\displaystyle =}
3
−
n
n
+
1
.
{\displaystyle {\frac {3-n}{n+1}}\,.}
This is precisely the relation that results from the definition of
α
≡
(
3
−
4
/
γ
g
)
{\displaystyle \alpha \equiv (3-4/\gamma _{g})}
if the model is evolved assuming
γ
g
=
(
n
+
1
)
/
n
{\displaystyle \gamma _{g}=(n+1)/n}
. We simultaneously seek the relation,
0
{\displaystyle 0}
=
{\displaystyle =}
(
n
+
1
)
α
(
a
0
)
+
b
0
(
n
−
1
)
{\displaystyle (n+1)\alpha (a_{0})+b_{0}(n-1)}
⇒
b
0
{\displaystyle \Rightarrow ~~~b_{0}}
=
{\displaystyle =}
[
n
+
1
1
−
n
]
α
(
a
0
)
{\displaystyle {\biggl [}{\frac {n+1}{1-n}}{\biggr ]}\alpha (a_{0})}
=
{\displaystyle =}
a
0
[
3
−
n
1
−
n
]
{\displaystyle a_{0}{\biggl [}{\frac {3-n}{1-n}}{\biggr ]}}
It appears as though the leading coefficient,
a
0
{\displaystyle a_{0}}
, is arbitrary, so we will set it equal to unity. This means that the displacement function is,
f
{\displaystyle f}
=
{\displaystyle =}
1
+
[
3
−
n
1
−
n
]
z
.
{\displaystyle 1+{\biggl [}{\frac {3-n}{1-n}}{\biggr ]}z\,.}
This expression for the displacement function,
f
{\displaystyle f}
, is identical to the expression found inside the square brackets of our separately derived exact solution of the polytropic LAWE . Furthermore, given the notation,
(
σ
c
2
/
γ
g
)
=
(
F
−
2
α
)
{\displaystyle (\sigma _{c}^{2}/\gamma _{g})=({\mathfrak {F}}-2\alpha )}
, the first term on the RHS of the LAWE will go to zero when,
F
=
2
(
3
−
n
)
/
(
n
+
1
)
{\displaystyle {\mathfrak {F}}=2(3-n)/(n+1)}
.
Part IV
If we divide through by
(
θ
n
−
1
z
)
{\displaystyle (\theta ^{n-1}z)}
, the LAWE that was derived above in Part II assumes the following form,
[
θ
1
−
n
z
−
1
]
×
{\displaystyle {\biggl [}\theta ^{1-n}z^{-1}{\biggr ]}\times }
LAWE
=
{\displaystyle =}
θ
1
−
n
z
−
1
[
ξ
−
2
(
1
+
3
z
)
2
+
2
n
(
1
+
3
z
)
z
2
θ
n
−
1
+
n
2
ξ
2
θ
2
(
n
−
1
)
z
4
]
⋅
d
2
f
d
z
2
+
{
(
n
−
1
)
+
(
n
−
3
)
z
}
⋅
d
f
d
z
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
n
z
−
1
+
α
]
⋅
f
{\displaystyle \theta ^{1-n}z^{-1}{\biggl [}\xi ^{-2}(1+3z)^{2}+2n(1+3z)z^{2}\theta ^{n-1}+n^{2}\xi ^{2}\theta ^{2(n-1)}z^{4}{\biggr ]}\cdot {\frac {d^{2}f}{dz^{2}}}+{\biggl \{}(n-1)+(n-3)z{\biggr \}}\cdot {\frac {df}{dz}}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-n}z^{-1}+\alpha {\biggr ]}\cdot f}
=
{\displaystyle =}
[
ξ
−
2
θ
1
−
n
z
−
1
(
1
+
3
z
)
2
+
2
n
(
1
+
3
z
)
z
+
n
2
ξ
2
θ
(
n
−
1
)
z
3
]
⋅
d
2
f
d
z
2
+
{
(
n
−
1
)
+
(
n
−
3
)
z
}
⋅
d
f
d
z
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
n
z
−
1
+
α
]
⋅
f
,
{\displaystyle {\biggl [}\xi ^{-2}\theta ^{1-n}z^{-1}(1+3z)^{2}+2n(1+3z)z+n^{2}\xi ^{2}\theta ^{(n-1)}z^{3}{\biggr ]}\cdot {\frac {d^{2}f}{dz^{2}}}+{\biggl \{}(n-1)+(n-3)z{\biggr \}}\cdot {\frac {df}{dz}}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-n}z^{-1}+\alpha {\biggr ]}\cdot f\,,}
which resembles the above-discussed hypergeometric differential equation , namely,
0
{\displaystyle 0}
=
{\displaystyle =}
z
(
1
−
z
)
d
2
u
d
z
2
+
[
γ
−
(
α
+
β
+
1
)
z
]
d
u
d
z
−
α
β
u
.
{\displaystyle z(1-z){\frac {d^{2}u}{dz^{2}}}+[\gamma -(\alpha +\beta +1)z]{\frac {du}{dz}}-\alpha \beta u\,.}
For the record we note that the coefficient (in square brackets) of the first term on the RHS of our LAWE expression is the square of the first derivative of
z
{\displaystyle z}
with respect to
ξ
{\displaystyle \xi }
; that is,
θ
1
−
n
z
−
1
(
d
z
d
ξ
)
2
{\displaystyle \theta ^{1-n}z^{-1}{\biggl (}{\frac {dz}{d\xi }}{\biggr )}^{2}}
=
{\displaystyle =}
θ
1
−
n
z
−
1
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
2
{\displaystyle \theta ^{1-n}z^{-1}{\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}^{2}}
=
{\displaystyle =}
θ
1
−
n
z
−
1
[
ξ
−
2
(
1
+
3
z
)
2
+
2
n
(
1
+
3
z
)
θ
n
−
1
z
2
+
n
2
ξ
2
θ
2
(
n
−
1
)
z
4
]
.
{\displaystyle \theta ^{1-n}z^{-1}{\biggl [}\xi ^{-2}(1+3z)^{2}+2n(1+3z)\theta ^{n-1}z^{2}+n^{2}\xi ^{2}\theta ^{2(n-1)}z^{4}{\biggr ]}\,.}
Part V
Now suppose that,
f
=
(
a
0
+
b
0
z
+
c
0
z
2
)
{\displaystyle f=(a_{0}+b_{0}z+c_{0}z^{2})}
, where again,
z
{\displaystyle z}
=
{\displaystyle =}
ξ
−
1
θ
−
n
(
θ
′
)
.
{\displaystyle \xi ^{-1}\theta ^{-n}(\theta ^{'})\,.}
Recalling that,
Part I Summary …
(
θ
′
)
{\displaystyle (\theta ^{'})}
=
{\displaystyle =}
ξ
θ
n
z
,
{\displaystyle \xi \theta ^{n}z\,,}
d
z
d
ξ
{\displaystyle {\frac {dz}{d\xi }}}
=
{\displaystyle =}
−
[
ξ
−
1
(
1
+
3
z
)
+
n
ξ
θ
n
−
1
z
2
]
,
{\displaystyle -{\biggl [}\xi ^{-1}(1+3z)+n\xi \theta ^{n-1}z^{2}{\biggr ]}\,,}
d
2
z
d
ξ
2
{\displaystyle {\frac {d^{2}z}{d\xi ^{2}}}}
=
{\displaystyle =}
4
ξ
−
2
(
1
+
3
z
)
+
2
n
θ
n
−
1
[
1
+
4
z
]
z
+
n
(
n
+
1
)
ξ
2
θ
2
(
n
−
1
)
z
3
.
{\displaystyle 4\xi ^{-2}(1+3z)+2n\theta ^{n-1}[1+4z]z+n(n+1)\xi ^{2}\theta ^{2(n-1)}z^{3}\,.}
it may prove useful to recognize that,
θ
′
=
d
θ
d
ξ
{\displaystyle \theta ^{'}={\frac {d\theta }{d\xi }}}
→
{\displaystyle \rightarrow }
d
z
d
ξ
⋅
d
θ
d
z
{\displaystyle {\frac {dz}{d\xi }}\cdot {\frac {d\theta }{dz}}}
⇒
d
θ
d
z
{\displaystyle \Rightarrow ~~~{\frac {d\theta }{dz}}}
=
{\displaystyle =}
(
d
z
d
ξ
)
−
1
ξ
θ
n
z
{\displaystyle {\biggl (}{\frac {dz}{d\xi }}{\biggr )}^{-1}\xi \theta ^{n}z}
⇒
d
ln
θ
d
ln
z
{\displaystyle \Rightarrow ~~~{\frac {d\ln \theta }{d\ln z}}}
=
{\displaystyle =}
(
d
z
d
ξ
)
−
1
ξ
θ
n
−
1
z
2
{\displaystyle {\biggl (}{\frac {dz}{d\xi }}{\biggr )}^{-1}\xi \theta ^{n-1}z^{2}}
⇒
d
z
d
ξ
{\displaystyle \Rightarrow ~~~{\frac {dz}{d\xi }}}
=
{\displaystyle =}
(
d
ln
θ
d
ln
z
)
−
1
ξ
θ
n
−
1
z
2
{\displaystyle {\biggl (}{\frac {d\ln \theta }{d\ln z}}{\biggr )}^{-1}\xi \theta ^{n-1}z^{2}}
In this case we have,
LAWE
=
{\displaystyle =}
[
−
(
d
ln
θ
d
ln
z
)
−
1
ξ
θ
n
−
1
z
2
]
2
⋅
(
2
c
0
)
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
1
+
α
θ
n
−
1
z
]
(
a
0
+
b
0
z
+
c
0
z
2
)
{\displaystyle {\biggl [}-{\biggl (}{\frac {d\ln \theta }{d\ln z}}{\biggr )}^{-1}\xi \theta ^{n-1}z^{2}{\biggr ]}^{2}\cdot (2c_{0})+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-1}+\alpha \theta ^{n-1}z{\biggr ]}(a_{0}+b_{0}z+c_{0}z^{2})}
+
{
(
n
−
1
)
+
(
n
−
3
)
z
}
θ
n
−
1
z
⋅
(
b
0
+
2
c
0
z
)
{\displaystyle +{\biggl \{}(n-1)+(n-3)z{\biggr \}}\theta ^{n-1}z\cdot (b_{0}+2c_{0}z)}
Useful ?????
Try again …
[
θ
1
−
n
z
−
1
]
×
{\displaystyle {\biggl [}\theta ^{1-n}z^{-1}{\biggr ]}\times }
LAWE
=
{\displaystyle =}
θ
1
−
n
z
−
1
[
d
z
d
ξ
]
2
⋅
d
2
(
a
0
+
b
0
z
+
c
0
z
2
)
d
z
2
+
{
(
n
−
1
)
+
(
n
−
3
)
z
}
⋅
d
(
a
0
+
b
0
z
+
c
0
z
2
)
d
z
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
n
z
−
1
+
α
]
⋅
(
a
0
+
b
0
z
+
c
0
z
2
)
{\displaystyle \theta ^{1-n}z^{-1}{\biggl [}{\frac {dz}{d\xi }}{\biggr ]}^{2}\cdot {\frac {d^{2}(a_{0}+b_{0}z+c_{0}z^{2})}{dz^{2}}}+{\biggl \{}(n-1)+(n-3)z{\biggr \}}\cdot {\frac {d(a_{0}+b_{0}z+c_{0}z^{2})}{dz}}+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-n}z^{-1}+\alpha {\biggr ]}\cdot (a_{0}+b_{0}z+c_{0}z^{2})}
=
{\displaystyle =}
[
ξ
−
2
(
1
+
3
z
)
2
+
2
n
(
1
+
3
z
)
z
2
θ
n
−
1
+
n
2
ξ
2
θ
2
(
n
−
1
)
z
4
]
⋅
(
2
c
0
)
+
{
(
n
−
1
)
+
(
n
−
3
)
z
}
⋅
(
b
0
+
2
c
0
z
)
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
n
z
−
1
+
α
]
⋅
(
a
0
+
b
0
z
+
c
0
z
2
)
{\displaystyle {\biggl [}\xi ^{-2}(1+3z)^{2}+2n(1+3z)z^{2}\theta ^{n-1}+n^{2}\xi ^{2}\theta ^{2(n-1)}z^{4}{\biggr ]}\cdot (2c_{0})+{\biggl \{}(n-1)+(n-3)z{\biggr \}}\cdot (b_{0}+2c_{0}z)+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-n}z^{-1}+\alpha {\biggr ]}\cdot (a_{0}+b_{0}z+c_{0}z^{2})}
=
{\displaystyle =}
[
(
2
c
0
)
ξ
−
2
(
1
+
6
z
+
9
z
2
)
+
4
c
0
n
(
z
2
+
3
z
3
)
θ
n
−
1
+
2
c
0
n
2
ξ
2
θ
2
(
n
−
1
)
z
4
]
+
(
n
−
1
)
(
b
0
+
2
c
0
z
)
{\displaystyle {\biggl [}(2c_{0})\xi ^{-2}(1+6z+9z^{2})+4c_{0}n(z^{2}+3z^{3})\theta ^{n-1}+2c_{0}n^{2}\xi ^{2}\theta ^{2(n-1)}z^{4}{\biggr ]}+(n-1)(b_{0}+2c_{0}z)}
+
(
n
−
3
)
(
b
0
z
+
2
c
0
z
2
)
+
(
n
+
1
)
[
(
σ
c
2
6
γ
g
)
θ
−
n
z
−
1
]
⋅
(
a
0
+
b
0
z
+
c
0
z
2
)
+
(
n
+
1
)
α
(
a
0
+
b
0
z
+
c
0
z
2
)
.
{\displaystyle +(n-3)(b_{0}z+2c_{0}z^{2})+(n+1){\biggl [}{\biggl (}{\frac {\sigma _{c}^{2}}{6\gamma _{\mathrm {g} }}}{\biggr )}\theta ^{-n}z^{-1}{\biggr ]}\cdot (a_{0}+b_{0}z+c_{0}z^{2})+(n+1)\alpha (a_{0}+b_{0}z+c_{0}z^{2})\,.}
Looks pretty hopeless!
Example Density- and Pressure-Profiles
Properties of Analytically Defined, Spherically Symmetric, Equilibrium Structures
Model
ρ
(
x
)
{\displaystyle \rho (x)}
P
(
x
)
{\displaystyle P(x)}
P
′
(
x
)
{\displaystyle P^{'}(x)}
μ
(
x
)
{\displaystyle \mu (x)}
ρ
(
x
)
P
(
x
)
{\displaystyle {\frac {\rho (x)}{P(x)}}}
Uniform-density
1
{\displaystyle 1}
1
−
x
2
{\displaystyle 1-x^{2}}
−
2
x
{\displaystyle -2x}
2
x
2
(
1
−
x
2
)
{\displaystyle {\frac {2x^{2}}{(1-x^{2})}}}
1
(
1
−
x
2
)
{\displaystyle {\frac {1}{(1-x^{2})}}}
Linear
1
−
x
{\displaystyle 1-x}
(
1
−
x
)
2
(
1
+
2
x
−
9
5
x
2
)
{\displaystyle (1-x)^{2}(1+2x-{\tfrac {9}{5}}x^{2})}
−
12
5
x
(
1
−
x
)
(
4
−
3
x
)
{\displaystyle -{\tfrac {12}{5}}x(1-x)(4-3x)}
12
5
x
2
(
4
−
3
x
)
(
1
−
x
)
(
1
+
2
x
−
9
5
x
2
)
{\displaystyle {\frac {{\tfrac {12}{5}}x^{2}(4-3x)}{(1-x)(1+2x-{\tfrac {9}{5}}x^{2})}}}
1
(
1
−
x
)
(
1
+
2
x
−
9
5
x
2
)
{\displaystyle {\frac {1}{(1-x)(1+2x-{\tfrac {9}{5}}x^{2})}}}
Parabolic
1
−
x
2
{\displaystyle 1-x^{2}}
(
1
−
x
2
)
2
(
1
−
1
2
x
2
)
{\displaystyle (1-x^{2})^{2}(1-{\tfrac {1}{2}}x^{2})}
−
x
(
1
−
x
2
)
(
5
−
3
x
2
)
{\displaystyle -x(1-x^{2})(5-3x^{2})}
x
2
(
5
−
3
x
2
)
(
1
−
x
2
)
(
1
−
1
2
x
2
)
{\displaystyle {\frac {x^{2}(5-3x^{2})}{(1-x^{2})(1-{\tfrac {1}{2}}x^{2})}}}
1
(
1
−
x
2
)
(
1
−
1
2
x
2
)
{\displaystyle {\frac {1}{(1-x^{2})(1-{\tfrac {1}{2}}x^{2})}}}
n
=
1
{\displaystyle ~n=1}
Polytrope
sin
x
x
{\displaystyle {\frac {\sin x}{x}}}
(
sin
x
x
)
2
{\displaystyle {\biggl (}{\frac {\sin x}{x}}{\biggr )}^{2}}
2
x
[
cos
x
−
sin
x
x
]
sin
x
x
{\displaystyle {\frac {2}{x}}{\biggl [}\cos x-{\frac {\sin x}{x}}{\biggr ]}{\frac {\sin x}{x}}}
2
(
1
−
x
cot
x
)
{\displaystyle 2(1-x\cot x)}
x
sin
x
{\displaystyle {\frac {x}{\sin x}}}
Uniform Density
In the case of a uniform-density, incompressible configuration, the Kopal (1948) LAWE becomes,
0
{\displaystyle 0}
=
{\displaystyle =}
d
2
f
d
x
2
+
(
4
−
μ
)
x
⋅
d
f
d
x
+
[
(
ω
2
ρ
0
R
2
γ
g
P
0
)
−
α
μ
x
2
]
f
{\displaystyle {\frac {d^{2}f}{dx^{2}}}+{\frac {(4-\mu )}{x}}\cdot {\frac {df}{dx}}+{\biggl [}{\biggl (}{\frac {\omega ^{2}\rho _{0}R^{2}}{\gamma _{\mathrm {g} }P_{0}}}{\biggr )}-{\frac {\alpha \mu }{x^{2}}}{\biggr ]}f}
=
{\displaystyle =}
d
2
f
d
x
2
+
1
x
[
4
−
2
x
2
(
1
−
x
2
)
]
d
f
d
x
+
[
(
ω
2
ρ
c
R
2
γ
g
P
c
)
1
(
1
−
x
2
)
−
(
2
α
1
−
x
2
)
]
f
{\displaystyle {\frac {d^{2}f}{dx^{2}}}+{\frac {1}{x}}{\biggl [}4-{\frac {2x^{2}}{(1-x^{2})}}{\biggr ]}{\frac {df}{dx}}+{\biggl [}{\biggl (}{\frac {\omega ^{2}\rho _{c}R^{2}}{\gamma _{\mathrm {g} }P_{c}}}{\biggr )}{\frac {1}{(1-x^{2})}}-{\biggl (}{\frac {2\alpha }{1-x^{2}}}{\biggr )}{\biggr ]}f}
=
{\displaystyle =}
(
1
−
x
2
)
⋅
d
2
f
d
x
2
+
1
x
[
4
−
6
x
2
]
d
f
d
x
+
[
(
ω
2
ρ
c
R
2
γ
g
P
c
)
−
2
α
]
f
.
{\displaystyle (1-x^{2})\cdot {\frac {d^{2}f}{dx^{2}}}+{\frac {1}{x}}{\biggl [}4-6x^{2}{\biggr ]}{\frac {df}{dx}}+{\biggl [}{\biggl (}{\frac {\omega ^{2}\rho _{c}R^{2}}{\gamma _{\mathrm {g} }P_{c}}}{\biggr )}-2\alpha {\biggr ]}f\,.}
Given that, in the equilibrium state ,
ρ
c
R
2
P
c
{\displaystyle {\frac {\rho _{c}R^{2}}{P_{c}}}}
=
{\displaystyle =}
6
4
π
G
ρ
c
{\displaystyle {\frac {6}{4\pi G\rho _{c}}}}
we obtain the LAWE derived by 📚 T. E. Sterne (1937, MNRAS, Vol. 97, pp. 582 - 593) — see his equation (1.91) on p. 585 — namely,
0
{\displaystyle 0}
=
{\displaystyle =}
(
1
−
x
2
)
⋅
d
2
f
d
x
2
+
1
x
[
4
−
6
x
2
]
d
f
d
x
+
[
6
(
ω
2
4
π
γ
g
G
ρ
c
)
−
2
α
]
f
{\displaystyle (1-x^{2})\cdot {\frac {d^{2}f}{dx^{2}}}+{\frac {1}{x}}{\biggl [}4-6x^{2}{\biggr ]}{\frac {df}{dx}}+{\biggl [}6{\biggl (}{\frac {\omega ^{2}}{4\pi \gamma _{\mathrm {g} }G\rho _{c}}}{\biggr )}-2\alpha {\biggr ]}f}
=
{\displaystyle =}
(
1
−
x
2
)
⋅
d
2
f
d
x
2
+
1
x
[
4
−
6
x
2
]
d
f
d
x
+
F
f
,
{\displaystyle (1-x^{2})\cdot {\frac {d^{2}f}{dx^{2}}}+{\frac {1}{x}}{\biggl [}4-6x^{2}{\biggr ]}{\frac {df}{dx}}+{\mathfrak {F}}f\,,}
where,
F
{\displaystyle {\mathfrak {F}}}
≡
{\displaystyle \equiv }
[
6
(
ω
2
4
π
γ
g
G
ρ
c
)
−
2
α
]
.
{\displaystyle {\biggl [}6{\biggl (}{\frac {\omega ^{2}}{4\pi \gamma _{\mathrm {g} }G\rho _{c}}}{\biggr )}-2\alpha {\biggr ]}\,.}
This also matches, respectively, equations (8) and (9) of 📚 Z. Kopal (1948, Proc. NAS, Vol. 34, Issue 8, pp.377-384) , aside from what, we presume, is a type-setting error that appears in the numerator of the second term on the RHS of his equation (8):
(
4
−
x
2
)
{\displaystyle (4-x^{2})}
appears, whereas it should be
(
4
−
6
x
2
)
{\displaystyle (4-6x^{2})}
.
In order to see if this differential equation is of the same form as the hypergeometric expression , we'll make the substitution,
z
{\displaystyle z}
≡
{\displaystyle \equiv }
x
2
{\displaystyle x^{2}}
⇒
d
z
{\displaystyle \Rightarrow ~~~dz}
=
{\displaystyle =}
2
x
d
x
{\displaystyle 2xdx}
⇒
d
f
d
x
{\displaystyle \Rightarrow ~~~{\frac {df}{dx}}}
=
{\displaystyle =}
d
z
d
x
⋅
d
f
d
z
=
2
x
⋅
d
f
d
z
=
2
z
1
/
2
⋅
d
f
d
z
{\displaystyle {\frac {dz}{dx}}\cdot {\frac {df}{dz}}=2x\cdot {\frac {df}{dz}}=2z^{1/2}\cdot {\frac {df}{dz}}}
⇒
d
2
f
d
x
2
{\displaystyle \Rightarrow ~~~{\frac {d^{2}f}{dx^{2}}}}
=
{\displaystyle =}
2
z
1
/
2
⋅
d
d
z
[
2
z
1
/
2
⋅
d
f
d
z
]
=
2
z
1
/
2
[
z
−
1
/
2
⋅
d
f
d
z
+
2
z
1
/
2
⋅
d
2
f
d
z
2
]
=
[
2
⋅
d
f
d
z
+
4
z
⋅
d
2
f
d
z
2
]
,
{\displaystyle 2z^{1/2}\cdot {\frac {d}{dz}}{\biggl [}2z^{1/2}\cdot {\frac {df}{dz}}{\biggr ]}=2z^{1/2}{\biggl [}z^{-1/2}\cdot {\frac {df}{dz}}+2z^{1/2}\cdot {\frac {d^{2}f}{dz^{2}}}{\biggr ]}={\biggl [}2\cdot {\frac {df}{dz}}+4z\cdot {\frac {d^{2}f}{dz^{2}}}{\biggr ]}\,,}
in which case the 📚 Sterne (1937) LAWE may be rewritten as,
0
{\displaystyle 0}
=
{\displaystyle =}
(
1
−
z
)
[
2
⋅
d
f
d
z
+
4
z
⋅
d
2
f
d
z
2
]
+
1
z
1
/
2
[
4
−
6
z
]
2
z
1
/
2
d
f
d
z
+
F
f
{\displaystyle (1-z){\biggl [}2\cdot {\frac {df}{dz}}+4z\cdot {\frac {d^{2}f}{dz^{2}}}{\biggr ]}+{\frac {1}{z^{1/2}}}{\biggl [}4-6z{\biggr ]}2z^{1/2}{\frac {df}{dz}}+{\mathfrak {F}}f}
=
{\displaystyle =}
(
1
−
z
)
[
4
z
⋅
d
2
f
d
z
2
]
+
(
1
−
z
)
[
2
⋅
d
f
d
z
]
+
2
[
4
−
6
z
]
d
f
d
z
+
F
f
{\displaystyle (1-z){\biggl [}4z\cdot {\frac {d^{2}f}{dz^{2}}}{\biggr ]}+(1-z){\biggl [}2\cdot {\frac {df}{dz}}{\biggr ]}+2{\biggl [}4-6z{\biggr ]}{\frac {df}{dz}}+{\mathfrak {F}}f}
=
{\displaystyle =}
4
z
(
1
−
z
)
⋅
d
2
f
d
z
2
+
2
[
5
−
7
z
]
d
f
d
z
+
F
f
.
{\displaystyle 4z(1-z)\cdot {\frac {d^{2}f}{dz^{2}}}+2{\biggl [}5-7z{\biggr ]}{\frac {df}{dz}}+{\mathfrak {F}}f\,.}
This is, indeed, of the hypergeometric form if we set
(
α
,
β
;
γ
;
z
)
{\displaystyle (\alpha ,\beta ;\gamma ;z)}
γ
{\displaystyle \gamma }
=
{\displaystyle =}
5
2
,
{\displaystyle {\frac {5}{2}}\,,}
(
α
+
β
+
1
)
{\displaystyle (\alpha +\beta +1)}
=
{\displaystyle =}
7
2
,
{\displaystyle {\frac {7}{2}}\,,}
α
β
{\displaystyle \alpha \beta }
=
{\displaystyle =}
−
F
4
.
{\displaystyle -{\frac {\mathfrak {F}}{4}}\,.}
Combining this last pair of expressions gives,
0
{\displaystyle 0}
=
{\displaystyle =}
−
F
4
−
α
[
5
2
−
α
]
{\displaystyle -{\frac {\mathfrak {F}}{4}}-\alpha {\biggl [}{\frac {5}{2}}-\alpha {\biggr ]}}
=
{\displaystyle =}
α
2
−
(
5
2
)
α
−
F
4
{\displaystyle \alpha ^{2}-{\biggl (}{\frac {5}{2}}{\biggr )}\alpha -{\frac {\mathfrak {F}}{4}}}
⇒
α
{\displaystyle \Rightarrow ~~~\alpha }
=
{\displaystyle =}
1
2
{
5
2
±
[
(
5
2
)
2
−
F
]
1
/
2
}
{\displaystyle {\frac {1}{2}}{\biggl \{}{\frac {5}{2}}\pm {\biggl [}{\biggl (}{\frac {5}{2}}{\biggr )}^{2}-{\mathfrak {F}}{\biggr ]}^{1/2}{\biggr \}}}
=
{\displaystyle =}
5
4
{
1
±
[
1
−
(
4
25
)
F
]
1
/
2
}
;
{\displaystyle {\frac {5}{4}}{\biggl \{}1\pm {\biggl [}1-{\biggl (}{\frac {4}{25}}{\biggr )}{\mathfrak {F}}{\biggr ]}^{1/2}{\biggr \}}\,;}
and,
β
{\displaystyle \beta }
=
{\displaystyle =}
5
2
−
5
4
{
1
±
[
1
−
(
4
25
)
F
]
1
/
2
}
.
{\displaystyle {\frac {5}{2}}-{\frac {5}{4}}{\biggl \{}1\pm {\biggl [}1-{\biggl (}{\frac {4}{25}}{\biggr )}{\mathfrak {F}}{\biggr ]}^{1/2}{\biggr \}}\,.}
Example α = -1
If we set
α
=
−
1
{\displaystyle \alpha =-1}
, then the eigenvector is,
u
1
=
F
(
−
1
,
7
2
;
5
2
;
x
2
)
{\displaystyle u_{1}=F{\biggl (}-1,{\frac {7}{2}};{\frac {5}{2}};x^{2}{\biggr )}}
=
{\displaystyle =}
1
−
[
β
γ
]
x
2
=
1
−
(
7
5
)
x
2
;
{\displaystyle 1-{\biggl [}{\frac {\beta }{\gamma }}{\biggr ]}x^{2}=1-{\biggl (}{\frac {7}{5}}{\biggr )}x^{2}\,;}
and the corresponding eigenfrequency is obtained from the expression,
−
1
{\displaystyle -1}
=
{\displaystyle =}
5
4
{
1
±
[
1
−
(
4
25
)
F
]
1
/
2
}
{\displaystyle {\frac {5}{4}}{\biggl \{}1\pm {\biggl [}1-{\biggl (}{\frac {4}{25}}{\biggr )}{\mathfrak {F}}{\biggr ]}^{1/2}{\biggr \}}}
⇒
−
9
5
{\displaystyle \Rightarrow ~~~-{\frac {9}{5}}}
=
{\displaystyle =}
±
[
1
−
(
4
25
)
F
]
1
/
2
{\displaystyle \pm {\biggl [}1-{\biggl (}{\frac {4}{25}}{\biggr )}{\mathfrak {F}}{\biggr ]}^{1/2}}
⇒
3
4
5
2
{\displaystyle \Rightarrow ~~~{\frac {3^{4}}{5^{2}}}}
=
{\displaystyle =}
1
−
(
4
25
)
F
{\displaystyle 1-{\biggl (}{\frac {4}{25}}{\biggr )}{\mathfrak {F}}}
⇒
F
{\displaystyle \Rightarrow ~~~{\mathfrak {F}}}
=
{\displaystyle =}
(
5
2
4
)
[
1
−
3
4
5
2
]
=
1
4
[
5
2
−
3
4
]
=
14
.
{\displaystyle {\biggl (}{\frac {5^{2}}{4}}{\biggr )}{\biggl [}1-{\frac {3^{4}}{5^{2}}}{\biggr ]}={\frac {1}{4}}{\biggl [}5^{2}-3^{4}{\biggr ]}=14\,.}
As we have reviewed in a separate discussion , this is identical to the eigenvector identified by 📚 Sterne (1937) as mode "
j
=
1
{\displaystyle j=1}
".
More Generally
More generally, in agreement with 📚 Sterne (1937) , for any (positive integer) mode number,
0
≤
j
≤
∞
{\displaystyle 0\leq j\leq \infty }
, we find,
And, in terms of the hypergeometric function series , the corresponding eigenfunction is,
u
j
=
{\displaystyle u_{j}=}
=
{\displaystyle =}
F
(
α
j
,
β
j
;
5
2
;
x
2
)
.
{\displaystyle F{\biggl (}\alpha _{j},\beta _{j};{\frac {5}{2}};x^{2}{\biggr )}\,.}
See Also