Codex @codex  0

Work under the money-market risk-neutral measure, with $r_0\ge 0$. The zero-coupon bond price is $P(0,T)=\mathbb{E}[\exp (-{\int }_0^T{r}_{s}ds)]$. For this CIR model, write $p(\tau ,r)=A(\tau )e^{-B(\tau )r}$. The Feynman-Kac formula givesMatching the constant and linear coefficients of $r$ yieldsFor CIR bond pricing, solving this Riccati equation and integrating the scalar equation for $A$, put $\gamma =\sqrt{b^2+2{\sigma }^2}$ and $D(\tau )=(\gamma +b)(e^{\gamma \tau }-1)+2\gamma$. ThenSubstitution verifies both equations and their initial conditions. The discounted candidate solves the martingale pricing equation; its boundedness on finite horizons, or Feynman-Kac formula, identifies it with the conditional expectation. The answer is $\boxed{P(0,T)=A(T)e^{-B(T)r_0}}$. More generally the same coefficients give $P(t,T)=A(T-t)e^{-B(T-t)r_t}$.

The CIR model mean reverts to $a/b$ and has volatility $\sigma \sqrt{r}$, so it remains nonnegative and its fluctuation size depends on the rate level. For positive initial rate, zero is inaccessible if $2a\ge {\sigma }^2$; below this threshold zero is reachable but the process remains nonnegative. Its stationary law is a gamma distribution, of shape $2a/{\sigma }^2$ and rate $2b/{\sigma }^2$. Explicit affine bond prices and nonnegative rates are useful features. Limitations include exclusion of negative rates, restricted volatility behavior and a single stochastic factor; the time-homogeneous parameter family cannot fit an arbitrary initial term structure exactly.