OFFSET
1,3
COMMENTS
See A195304 for definitions and a general discussion.
EXAMPLE
(B)=1.0684730171986100357539820259956...
MATHEMATICA
a = 1; b = Sqrt[GoldenRatio]; h = 2 a/3; k = b/3;
f[t_] := (t - a)^2 + ((t - a)^2) ((a*k - b*t)/(a*h - a*t))^2
s = NSolve[D[f[t], t] == 0, t, 150]
f1 = (f[t])^(1/2) /. Part[s, 4]
RealDigits[%, 10, 100] (* (A) A195491 *)
f[t_] := (t - a)^2 + ((t - a)^2) (k/(h - t))^2
s = NSolve[D[f[t], t] == 0, t, 150]
f2 = (f[t])^(1/2) /. Part[s, 4]
RealDigits[%, 10, 100] (* (B) A195492 *)
f[t_] := (b*t/a)^2 + ((b*t/a)^2) ((a*h - a*t)/(b*t - a*k))^2
s = NSolve[D[f[t], t] == 0, t, 150]
f3 = (f[t])^(1/2) /. Part[s, 4]
RealDigits[%, 10, 100] (* (C) A195493 *)
c = Sqrt[a^2 + b^2]; (f1 + f2 + f3)/(a + b + c)
RealDigits[%, 10, 100] (* Philo(ABC, G) A195494 *)
PROG
(PARI) polrootsreal(3464*x^8 - 16028*x^7 + 40544*x^6 + 47079*x^5 - 64883*x^4 + 50115*x^3 - 43008*x^2 - 60658*x + 32292)[4] \\ Charles R Greathouse IV, Nov 26 2024
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Clark Kimberling, Sep 19 2011
STATUS
approved
