cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A197035 Decimal expansion of the shortest distance from the x axis through (3,1) to the line y=x.

Original entry on oeis.org

2, 6, 0, 8, 1, 9, 4, 0, 2, 4, 9, 6, 1, 0, 1, 8, 9, 0, 1, 9, 9, 0, 1, 4, 4, 5, 4, 2, 8, 3, 5, 2, 2, 3, 9, 5, 9, 0, 8, 3, 5, 8, 9, 1, 1, 5, 8, 7, 9, 5, 9, 7, 6, 7, 4, 4, 9, 4, 9, 1, 7, 7, 5, 6, 3, 8, 5, 7, 7, 4, 4, 9, 2, 8, 8, 4, 1, 8, 9, 9, 6, 8, 0, 3, 9, 1, 0, 4, 3, 1, 5, 5, 9, 0, 5, 1, 4, 5, 0
Offset: 1

Views

Author

Clark Kimberling, Oct 10 2011

Keywords

Comments

The shortest segment from one side of an angle T through a point P inside T is called the Philo line of P in T. For discussions and guides to related sequences, see A197008 and A195284.

Examples

			length of Philo line:  2.60819402496101...
endpoint on x axis:   (3.47797, 0); see A197034
endpoint on line y=x: (2.35321, 2.35321)
		

Crossrefs

Programs

  • Mathematica
    f[t_] := (t - k*t/(k + m*t - m*h))^2 + (m*k*t/(k + m*t - m*h))^2;
    g[t_] := D[f[t], t]; Factor[g[t]]
    p[t_] :=  h^2 k + k^3 - h^3 m - h k^2 m - 3 h k t + 3 h^2 m t + 2 k t^2 - 3 h m t^2 + m t^3
    m = 1; h = 3; k = 1;  (* slope m; point (h,k) *)
    t = t1 /. FindRoot[p[t1] == 0, {t1, 1, 2}, WorkingPrecision -> 100]
    RealDigits[t]  (* A197034 *)
    {N[t], 0} (* endpoint on x axis *)
    {N[k*t/(k + m*t - m*h)], N[m*k*t/(k + m*t - m*h)]} (* upper endpoint *)
    d = N[Sqrt[f[t]], 100]
    RealDigits[d] (* A197035 *)
    Show[Plot[{k*(x - t)/(h - t), m*x}, {x, 0, 4}],
    ContourPlot[(x - h)^2 + (y - k)^2 == .002, {x, 0, 3.5}, {y, 0, 3}], PlotRange -> {0, 3}, AspectRatio -> Automatic]