A = {{a₀, a₁}, {a₁, a₂}}, b = {b₀, b₁}, x₀ = {0, 0}

Step 1: r₀ = b - Ax₀; p₀ = r₀

r[0] = b₀

r[1] = b₁

Step 2: alpha₀ = <r₀, r₀> / <p₀, Ap₀>

<p₀, Ap₀> = b₀(a₀b₀ + a₁b₁) + b₁(a₁b₀ + a₂b₁)

alpha₀ = (b₀²+b₁²) / (b₀(a₀b₀ + a₁b₁) + b₁(a₁b₀ + a₂b₁))

Step 3: r₁ = r₀ - alpha₀ Ap₀

r[0] = (-b₁(-a₁b₀² + a₀b₀b₁ - a₂b₀b₁ + a₁b₁²)) / (b₀(a₀b₀ + a₁b₁) + b₁(a₁b₀ + a₂b₁))

r[1] = (b₀(-a₁b₀² + a₀b₀b₁ - a₂b₀b₁ + a₁b₁²)) / (b₀(a₀b₀ + a₁b₁) + b₁(a₁b₀ + a₂b₁))

Step 4: x₁ = x₀ + alphap₀

x[0] = (b₀(b₀² + b₁²)) / (b₀(a₀b₀ + a₁b₁) + b₁(a₁b₀ + a₂b₁))

x[1] = (b₁(b₀² + b₁²)) / (b₀(a₀b₀ + a₁b₁) + b₁(a₁b₀ + a₂b₁))

Step 5: beta = rsnew / rsold = <r₁, r₁> / <r₀, r₀>

rsnew = ((b₀² + b₁²)(-a₁b₀² + a₀b₀b₁ - a₂b₀b₁ + a₁b₁²)²) / (b₀(a₀b₀ + a₁b₁) + b₁(a₁b₀ + a₂b₁))²

beta = (a₁b₀² - a₀b₀b₁ + a₂b₀b₁ - a₁b₁²)² / (b₀(a₀b₀ + a₁b₁) + b₁(a₁b₀ + a₂b₁))²

Step 6: p₁ = r₁ +beta p₀

p[0] = (-1)((a₁b₀ + a₂b₁)(b₀² + b₁²)(-a₁b₀² + a₀b₀b₁ - a₂b₀b₁ + a₁b₁²)) / (b₀(a₀b₀ + a₁b₁) + b₁(a₁b₀ + a₂b₁))²

p[1] = ((a₀b₀ + a₁b₁)(b₀² + b₁²)(-a₁b₀² + a₀b₀b₁ - a₂b₀b₁ + a₁b₁²)) / (b₀(a₀b₀ + a₁b₁) + b₁(a₁b₀ + a₂b₁))²

Step 7: alpha₁ = <r₁, r₁> / <p₁, Ap₁>

<p₁, Ap₁> = ((-a₁² + a₀a₂)(b₀² + b₁²)²(-a₁b₀² + a₀b₀b₁ - a₂b₀b₁ + a₁b₁²)²) / (b₀(a₀b₀ + a₁b₁) + b₁(a₁b₀ + a₂b₁))³

alpha₁ = (a₀b₀² + 2a₁b₀b₁ + a₂b₁²) / ((-a₁² + a₀a₂) (b₀² + b₁²))

Step 8: r₂ = r₁ - alpha₁ Ap₁

r[0] = 0

r[1] = 0

Step 9: x₂ = x₁ + alpha₁p₁

x[0] = (a₂b₀ - a₁b₁) / (a₀a₂ - a₁²)

x[1] = (-a₁b₀ + a₀b₁) / (a₀a₂ - a₁²)

assertion:

bncg[0] = A[0][0]x[0] + A[0][1]x[1] = a₀(a₂b₀ - a₁b₁) / (a₀a₂ - a₁²) + a₁(-a₁b₀ + a₀b₁) / (a₀a₂ - a₁²) = b₀(a₀a₂ - a₁²) / (a₀a₂ - a₁²) = b₀

bncg[1] = A[1][0]x[0] + A[1][1]x[1] = a₁(a₂b₀ - a₁b₁) / (a₀a₂ - a₁²) + a₂(-a₁b₀ + a₀b₁) / (a₀a₂ - a₁²) = b₁(a₀a₂ - a₁²) / (a₀a₂ - a₁²) = b₁

b[0] = b₀

b[1] = b₁

assert(bncg[i] == b[i])

END