login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A195296
Decimal expansion of shortest length, (C), of segment from side CA through incenter to side CB in right triangle ABC with sidelengths (a,b,c)=(8,15,17).
4
6, 9, 9, 7, 1, 4, 2, 2, 7, 3, 8, 1, 4, 3, 6, 0, 5, 6, 5, 0, 4, 8, 9, 8, 3, 4, 5, 3, 0, 5, 4, 6, 9, 9, 6, 9, 1, 8, 2, 5, 6, 7, 8, 0, 0, 1, 8, 6, 3, 1, 6, 6, 3, 4, 5, 3, 4, 0, 0, 0, 8, 0, 9, 3, 8, 4, 1, 3, 7, 2, 0, 7, 4, 3, 9, 6, 0, 5, 5, 1, 5, 3, 1, 9, 8, 2, 8, 8, 3, 9, 1, 7, 4, 3, 6, 4, 2, 4, 7
OFFSET
1,1
COMMENTS
See A195284 for definitions and a general discussion.
EXAMPLE
(C)=6.99714227381436056504898345305469969182567800...
MATHEMATICA
a = 8; b = 15; c = 17;
h = a (a + c)/(a + b + c); k = a*b/(a + b + c);
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) A195293 *)
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, 1]
RealDigits[%, 10, 100] (* (B) A195296 *)
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] (* (C) A010524 *)
(f1 + f2 + f3)/(a + b + c)
RealDigits[%, 10, 100] (* Philo(ABC, I), A195297 *)
CROSSREFS
Cf. A195284.
Sequence in context: A246709 A051496 A330564 * A273951 A307053 A100403
KEYWORD
nonn,cons
AUTHOR
Clark Kimberling, Sep 14 2011
STATUS
approved