On $U$, one has $t=\sigma (r)$. We first construct two triples of coefficients, positive on the middle band, with product $f_*(\sigma )(\sigma ')^2$. At each endpoint of the band, exactly one of the coefficients of $P$ and $Q$ will vanish.
Put $d_*(r)=f_*(\sigma (r))(\sigma '(r))^2$. Since $a=\sigma (-1)$, $b=\sigma (1)$, define the positive smooth divided differences
\[ k_-(r)=\int _0^1\sigma '(-1+s(r+1))\, \mathrm{d}s,\qquad k_+(r)=\int _0^1\sigma '(1+s(r-1))\, \mathrm{d}s, \]
and $H_*(r)=k_-(r)k_+(r)\sigma '(r)>0$. The fundamental theorem of calculus gives
3FD3DD269BA8
14.29.5
\begin{equation} \sigma (r)-a=(r+1)k_-(r),\quad \sigma (r)-b=(r-1)k_+(r),\quad d_*=(1-r^2)^2H_*^2. \end{equation}
These formulas show that $H_*$ is smooth and positive also at the endpoints, without any division by a vanishing function. The symmetry $\sigma (-r)=1-\sigma (r)$ makes $d_*$, and hence $H_*$, even.
Choose the function
\[ p(r)= \begin{cases} \mathrm{e}^{-1/(1-r^2)},& |r|<1,\\ 0,& |r|\geq 1, \end{cases} \qquad q=p'. \]
All its derivatives vanish at $r=\pm 1$: each derivative is a finite sum of polynomials in $r$ times powers of $(1-r^2)^{-1}$, multiplied by the exponential, and the exponential factor dominates every such power. On $|r|<1$, $q=-2rp/(1-r^2)^2$. Define
14762A8A84D0
14.29.6
\begin{align} z_ b& =\sqrt{q^2+2H_*^2},\qquad S_* =\sqrt{r^2+4p^2+2z_ b(1-r^2)^2},\\ x_ b& =(S_*-r)/2,\qquad y_ b=(S_*+r)/2.\tag{14.29.7} \end{align}
Both radicands are strictly positive everywhere, so these formulas define smooth functions on $\mathbb R$, including at the endpoints. They satisfy
E916650EE31A
14.29.8
\begin{equation} z_ b>|q|,\qquad x_ b-y_ b=-r,\qquad x_ by_ b=p^2+\frac{z_ b}2(1-r^2)^2. \end{equation}
The coefficients $x_ b,y_ b$ are nonnegative, positive away from the endpoints, with values $(1,0)$ at $-1$ and $(0,1)$ at $1$. On the open band,
\[ S_*^2-(2p+|r|)^2=2(1-r^2)^2(z_ b-|q|)>0, \]
so $x_ b>p$, $y_ b>p$. Consequently both triples
025CE4264442
14.29.9
\begin{equation} (x_ b-p,y_ b+p,z_ b-q),\qquad (x_ b+p,y_ b-p,z_ b+q) \end{equation}
are strictly positive there.
We now verify that the two triples have the same product. Set $D_*=x_ by_ b-p^2$, $E_*=p(x_ b-y_ b)=-rp$. Equations (14.29.8) and the formula for $q$ give
\[ D_*=(z_ b/2)(1-r^2)^2,\qquad qD_*=z_ bE_*. \]
The two products in (14.29.9) therefore equal
17DDFD1D68B2
14.29.10
\begin{equation} (D_*+E_*)(z_ b-q)=(D_*-E_*)(z_ b+q) =\frac{D_*}{z_ b}(z_ b^2-q^2)=d_*. \end{equation}
For $|r|\geq 1$, $p=q=0$ and the same identity reads $x_ by_ bz_ b=d_*$. Thus both triples have exactly the prescribed product, and at the endpoints exactly one of their entries vanishes.
Comments
Corrections, questions, and mathematical discussion are welcome. Sign in with GitHub to comment. Comments are public.