eqn=function(a) 0.057564*a+0.197412*25-1.970*(0.00004707311*(a^2)-0.0002925398*a*25+0.0004771203*(25)^2)^0.5
x=seq(-100,100)
plot(x,eqn(x))
##从图上看，方程的解在a=-50附近，先大致确定好范围

> uniroot(eqn,c(-100,100))
$root
[1] -54.44617

$f.root
[1] -4.48938e-10

$iter
[1] 4

$init.it
[1] NA

$estim.prec
[1] 6.103516e-05
##方程的解为
a=-54.44617