OFFSET
1,1
COMMENTS
Decimal expansion of the unique positive real number c that minimizes the integral I(c) = integral_{x=0..1} (sin((Pi/2)x)^(1/c) - sqrt(2x-x^2))^2 dx.
Geometrically, the curve y = sin((Pi/2)x)^(1/c) is a one-parameter family that interpolates between a quarter-circle (at this optimal c) and a unit square as c -> infinity.
The function sqrt(2x-x^2) is the lower-right quarter of the unit circle centered at (1,0).
The optimal c gives the best least-squares fit to the quarter-circle in this family.
The value c is the unique solution to d/dx I(x) = 0, which must be solved numerically.
It is conjectured that this constant is transcendental.
EXAMPLE
2.277552113851191588974737802127065126412089221...
MATHEMATICA
$MaxExtraPrecision = 100;
I[a_?NumericQ] := N[Integrate[(Sin[Pi*x/2]^(1/a) - Sqrt[2x - x^2])^2, {x, 0, 1}, WorkingPrecision -> 50]];
res = FindMinimum[I[a], {a, 2}, WorkingPrecision -> 50];
optimalA = a /. res[[2]];
RealDigits[N[optimalA, 80]]
PROG
(PARI) f(a) = intnum(x=0, 1, (sin((Pi/2)*x)^(1/a) - sqrt(2*x-x^2))^2);
solve(a=2.27, 2.28, f'(a)) \\ Hugo Pfoertner, Mar 16 2026
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Minyan Liu, Mar 10 2026
STATUS
approved
