Solve for x : \begin{align} 25+2^{\log_{10} x}=x \tag{1}\label{1} \end{align}
Note that \eqref{1} is equivalent to
\begin{align}
x^{\log_{10}2}-x+25&=0
\tag{2}\label{2}
\\
\text{or }\quad
x^a-x+b&=0
\tag{3}\label{3}
\end{align}
with non-rational $a$.
It is known that such equations
don't have an algebraic solution and
can be solved only by means of numerical methods.
For example, we can use Halley's method
to iteratively find the approximation of the root as
\begin{align}
x_{n+1}&=F(x_n)
,\\
F(x)&=x-\frac{2\,f(x)\,f'(x)}{2f'(x)^2-f(x)\,f''(x)}
,\\
f(x)&=x^{\log_{10}(2)}-x+25
,\\
f'(x)&=\log_{10}(2)\cdot x^{\log_{10}(2)-1}-1
,\\
f''(x)&=\log_{10}(2)\log_{10}(\tfrac15)\cdot x^{\log_{10}(2)-2}
.
\end{align}
For example, starting with $x_0=1$, we get
\begin{align}
x_1&=6.60306336935\\
x_2&=26.5079884286\\
x_3&=27.7184046785\\
x_4&=27.7184201926\\
x_5&=27.7184201926\\
\end{align}
Edit
The rate of convergence to the root is cubic,
compare for example
to the Newton's method:
starting with the same $x_0$, the Halley's approximations would be
\begin{align}
x_0&=\color{blue}{ 27}.386363636363636363636 \\
x_1&=\color{blue}{ 27.7184}19892956254689994 \\
x_2&=\color{blue}{ 27.718420192574854316455}
\end{align}
Corresponding python code:
import decimal
decimal.getcontext().prec = 23
lg2 = decimal.Decimal(2).log10()
def f(x):
return 25+x**lg2-x
def df(x):
return lg2x*(lg2-1)-1
def ddf(x):
return lg2(lg2-1)x**(lg2-2)
def F(x):
fx=f(x)
dfx=df(x)
ddfx=ddf(x)
return x-2fxdfx/(2dfx2-fxddfx)
x=decimal.getcontext().divide(1205,44); print(x)
x=F(x); print(x)
x=F(x); print(x)
x=F(x); print(x)
27.386363636363636363636
27.718419892956254689994
27.718420192574854316455
27.718420192574854316455
$10^{log_{10} x}-2^{log_{10} x}=25$
$2^{log_{10} x}\big(5^{log_{10} x }-1\big)=25$
– sirous Nov 24 '20 at 17:38