OFFSET
1,3
COMMENTS
Let S and C denote the graphs of y = sin x and y = cos x. For each point (u, sin u) on S, let S(u) be the line normal to S at (u, sin u), and let (snc u, cos(snc u) be the point of intersection of S(u) and C. Let d(u) be the distance from (u,sin u) to (snc u, cos(snc u)). We call d(u) the u-normal distance from S to C and note that in [0,Pi], there is a unique number u' such that d(u') > d(u) for all real numbers u except those of the form u' + k*Pi. We call d(u') the maximal normal distance between sine and cosine, and we call snc the sine-normal-to-cosine function.
For each (v, cos v) on C, let C(u) be the line normal to C at (v, cos v), and let (cns v, sin(cns v)) be the point of intersection of C(u) and S. Let e(v) be the distance from (v, cos v) to (cns v, sin(cns v)). We call d(v) the v-normal distance from C to S and note that there exists a unique number v' that maximizes e, and e(v') = d(u'). We call cns the cosine-normal-to-sine function. The numbers u' and v' are given in A332501 and A332503.
Note that the maximal normal distance (see Example) exceeds the normal distance from (Pi/2,1) in sine to (Pi/2,0) in cosine - possibly a surprise!
FORMULA
d(u') = 2*sqrt((u - 3 Pi/4)^2 + (sin u)^2).
Equals sqrt(d^2+2-2*sqrt(1-d^2)) where d = A003957. - Gleb Koloskov, Jun 16 2021
EXAMPLE
2.72573705679992524967463858129656... = the number u in [0,2 Pi] such that the line normal to S at (u, sin u) passes through the point (3 Pi/4,0); cf. A332501.
0.4039727532995172093189617400663... = sin u; cf. A086751.
1.0949989843708724286504083007155... = maximal normal distance between sine and cosine.
1.9866519235847646080193264936226... = snc u; cf A332503.
MATHEMATICA
Plot[{Sin[x], Cos[x]}, {x, -Pi, 3 Pi}, AspectRatio -> Automatic,
ImageSize -> 600, PlotLabel -> "sine and cosine"]
t = Table[x = x /. FindRoot[Cos[x] == -x Sec[u] + u Sec[u] + Sin[u], {x, 0}], {u, -2 Pi, 2 Pi, Pi/101}];
ListPlot[t, PlotLabel -> "y \[Equal] snc(x)"]
ListPlot[Cos[t], PlotLabel -> "y \[Equal] cos(snc (x))"]
t = Table[x = x /. FindRoot[Sin[x] == x Csc[u] - u Csc[u] + Cos[u], {x, 0.1}], {u, -2 Pi + .01, 2 Pi - .01, Pi/101}];
ListPlot[t, PlotLabel -> "y \[Equal] cns(x)"]
ListPlot[Sin[t], PlotLabel -> "y \[Equal] sin(cns(x)"]
u = u /. FindRoot[0 == (-3 Pi/4) Sec[u] + u Sec[u] + Sin[u], {u, 1}, WorkingPrecision ->120] (* A332501 *)
y = Sin[u] (* A086751 *)
d = 2*Sqrt[(u - 3 Pi/4)^2 + y^2] (* A332500 *)
RealDigits[u][[1]] (* A332501 *)
RealDigits[y][[1]] (* A086751 *)
RealDigits[d][[1]] (* A332500 *)
PROG
(PARI) my(d=solve(x=0, 1, cos(x)-x)); sqrt(d^2+2-2*sqrt(1-d^2)) \\ Gleb Koloskov, Jun 16 2021
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Clark Kimberling, May 05 2020
STATUS
approved