Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Eqs. (29)-(30): gain rate and relative values of a two-state Markov-renewal process

Proved
JewellMRP.GainRate.two_state_gain_rate

by mikedeng1 · Sep 28, 2026 · Mathlib 0df444a (Lean v4.33.1)

average-rewardmarkov-renewalp2o-batch-p200bp2o-gran-per-chapterp2o-plan-paperp2o-v1

Consider a Markov-renewal program with two states and a stationary policy zzz whose chain is

P=(1−p12p12p211−p21),p12>0, p21>0,P = \begin{pmatrix} 1 - p_{12} & p_{12} \\ p_{21} & 1 - p_{21} \end{pmatrix}, \qquad p_{12} > 0,\ p_{21} > 0,P=(1−p12​p21​​p12​1−p21​​),p12​>0, p21​>0,

with mean sojourn times ν1,ν2>0\nu_1, \nu_2 > 0ν1​,ν2​>0 and expected rewards ρ1,ρ2\rho_1, \rho_2ρ1​,ρ2​. Let

g=p21ρ1+p12ρ2ν1p21+ν2p12,v1=ν2ρ1−ν1ρ2ν1p21+ν2p12,v2=0.g = \frac{p_{21}\rho_1 + p_{12}\rho_2}{\nu_1 p_{21} + \nu_2 p_{12}}, \qquad v_1 = \frac{\nu_2\rho_1 - \nu_1\rho_2}{\nu_1 p_{21} + \nu_2 p_{12}}, \qquad v_2 = 0 .g=ν1​p21​+ν2​p12​p21​ρ1​+p12​ρ2​​,v1​=ν1​p21​+ν2​p12​ν2​ρ1​−ν1​ρ2​​,v2​=0.

Then (g,v1,v2)(g, v_1, v_2)(g,v1​,v2​) solves the value-determination equations (13) of zzz, and ggg equals the gain rate (B 7) of zzz computed from any stationary probability vector of PPP.

These closed forms let the optimal policy of a two-state program be found by direct evaluation, and they give an explicit, nonvacuous instance of the objects of the algorithm of Fig. 2.

Formalization Note States 1,21, 21,2 are indices 0,10, 10,1 of Fin 2. The hypotheses p12,p21>0p_{12}, p_{21} > 0p12​,p21​>0 make the chain ergodic, as the paper notes below (21); no ergodicity of other policies is assumed.

Preamble
import Mathlib
import Definitions.Def_JewellMRP_GainRate_MRP
import Definitions.Def_JewellMRP_GainRate_PolicyIteration
Formal statement
namespace JewellMRP.GainRate
theorem two_state_gain_rate {α : Type*} (M : MRP 2 α) (z : Fin 2 → α)
    (h12 : 0 < M.p 0 (z 0) 1) (h21 : 0 < M.p 1 (z 1) 0) :
    let p12 := M.p 0 (z 0) 1
    let p21 := M.p 1 (z 1) 0
    let ν1 := M.ν 0 (z 0)
    let ν2 := M.ν 1 (z 1)
    let ρ1 := M.ρ 0 (z 0)
    let ρ2 := M.ρ 1 (z 1)
    let g := (p21 * ρ1 + p12 * ρ2) / (ν1 * p21 + ν2 * p12)
    M.SolvesValueDetermination z g ![(ν2 * ρ1 - ν1 * ρ2) / (ν1 * p21 + ν2 * p12), 0] ∧
      ∀ π : Fin 2 → ℝ, IsStationaryDist (M.policyMatrix z) π → M.gainRate z π = g := by sorry
end JewellMRP.GainRate
Source
Jewell, Markov-Renewal Programming. II: Infinite Return Models, Example, Oper. Res. 11 (1963), p. 961, Eqs. (29)-(30), with (21) (p. 960)
Human review
  • Endorsed by Shuze Chen · Oct 1, 2026

    Confirmed by the moderator at approval.

  • Endorsed by mikedeng1 · Oct 1, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me