A = x0, x1   b = x3[2]
      x1, x2

Step 1: r = b - Ax
  r[0] = x3[0]

  r[1] = x3[1]

Step 2: alpha = (r[i]*r[i]) / (p[i]*A[i][j]*p[j])
  alpha = (X3[0]^2^+X3[1]^2^) / ((X3[1]^2^)*X2+2*(X3[1]*X3[0]*X1)+(X3[0]^2)*X0)

Step 3: r[i] = r - alpha * A[i][j] * p[j]
  r[0] = (-1*((X3[1]^3)*X1)-1*((X3[1]^2)*X3[0]*X0)+(X3[1]^2)*X3[0]*X2+X3[1]*(X3[0]^2)*X1) / ((X3[1]^2)*X2+2*(X3[1]*X3[0]*X1)+(X3[0]^2)*X0)

  r[1] = ((X3[1]^2)*X3[0]*X1+X3[1]*(X3[0]^2)*X0-1*(X3[1]*(X3[0]^2)*X2)-1*((X3[0]^3)*X1)) / ((X3[1]^2)*X2+2*(X3[1]*X3[0]*X1)+(X3[0]^2)*X0) 

Step 4: x[i] = x + alpha*p[i]
  x[0] = ((X3[1]^2)*X3[0]+X3[0]^3) / ((X3[1]^2)*X2+2*(X3[1]*X3[0]*X1)+(X3[0]^2)*X0)

  x[1] = (X3[1]^3 + X3[1]*(X3[0]^2)) / ((X3[1]^2)*X2+2*(X3[1]*X3[0]*X1)+(X3[0]^2)*X0)

Step 5: beta = (rk[i]*rk[i]) / (r[i]*r[i])
  beta = ((X3[1]^6)*(X1^2)+2*((X3[1]^5)*X3[0]*X0*X1)-2*((X3[1]^5)*X3[0]*X1*X2)+(X3[1]^4)*(X3[0]^2)*(X0^2)-2*((X3[1]^4)*(X3[0]^2)*X0*X2)-1*((X3[1]^4)*(X3[0]^2)*(X1^2))+(X3[1]^4)*(X3[0]^2)*(X2^2)+(X3[1]^2)*(X3[0]^4)*(X0^2)-2*((X3[1]^2)*(X3[0]^4)*X0*X2)-1*((X3[1]^2)*(X3[0]^4)*(X1^2))+(X3[1]^2)*(X3[0]^4)*(X2^2)-2*(X3[1]*(X3[0]^5)*X0*X1)+2*(X3[1]*(X3[0]^5)*X1*X2)+(X3[0]^6)*(X1^2))/((X3[1]^6)*(X2^2)+4*((X3[1]^5)*X3[0]*X1*X2)+2*((X3[1]^4)*(X3[0]^2)*X0*X2)+4*((X3[1]^4)*(X3[0]^2)*(X1^2))+(X3[1]^4)*(X3[0]^2)*(X2^2)+4*((X3[1]^3)*(X3[0]^3)*X0*X1)+4*((X3[1]^3)*(X3[0]^3)*X1*X2)+(X3[1]^2)*(X3[0]^4)*(X0^2)+2*((X3[1]^2)*(X3[0]^4)*X0*X2)+4*((X3[1]^2)*(X3[0]^4)*(X1^2))+4*(X3[1]*(X3[0]^5)*X0*X1)+(X3[0]^6)*(X0^2))

Step 6: p[i] = rk[i] +beta * p
  p[0] = (-1*((X3[1]^7)*X1*X2)-1*((X3[1]^6)*X3[0]*X0*X2)-1*((X3[1]^6)*X3[0]*(X1^2))+(X3[1]^6)*X3[0]*(X2^2)-1*((X3[1]^5)*(X3[0]^2)*X0*X1)-2*((X3[1]^4)*(X3[0]^3)*X0*X2)-1*((X3[1]^4)*(X3[0]^3)*(X1^2))+2*((X3[1]^4)*(X3[0]^3)*(X2^2))-2*((X3[1]^3)*(X3[0]^4)*X0*X1)+3*((X3[1]^3)*(X3[0]^4)*X1*X2)-1*((X3[1]^2)*(X3[0]^5)*X0*X2)+(X3[1]^2)*(X3[0]^5)*(X1^2)+(X3[1]^2)*(X3[0]^5)*(X2^2)-1*(X3[1]*(X3[0]^6)*X0*X1)+2*(X3[1]*(X3[0]^6)*X1*X2)+(X3[0]^7)*(X1^2))/((X3[1]^6)*(X2^2)+4*((X3[1]^5)*X3[0]*X1*X2)+2*((X3[1]^4)*(X3[0]^2)*X0*X2)+4*((X3[1]^4)*(X3[0]^2)*(X1^2))+(X3[1]^4)*(X3[0]^2)*(X2^2)+4*((X3[1]^3)*(X3[0]^3)*X0*X1)+4*((X3[1]^3)*(X3[0]^3)*X1*X2)+(X3[1]^2)*(X3[0]^4)*(X0^2)+2*((X3[1]^2)*(X3[0]^4)*X0*X2)+4*((X3[1]^2)*(X3[0]^4)*(X1^2))+4*(X3[1]*(X3[0]^5)*X0*X1)+(X3[0]^6)*(X0^2))
  
  p[1] = ((X3[1]^7)*(X1^2)+2*((X3[1]^6)*X3[0]*X0*X1)-1*((X3[1]^6)*X3[0]*X1*X2)+(X3[1]^5)*(X3[0]^2)*(X0^2)-1*((X3[1]^5)*(X3[0]^2)*X0*X2)+(X3[1]^5)*(X3[0]^2)*(X1^2)+3*((X3[1]^4)*(X3[0]^3)*X0*X1)-2*((X3[1]^4)*(X3[0]^3)*X1*X2)+2*((X3[1]^3)*(X3[0]^4)*(X0^2))-2*((X3[1]^3)*(X3[0]^4)*X0*X2)-1*((X3[1]^3)*(X3[0]^4)*(X1^2))-1*((X3[1]^2)*(X3[0]^5)*X1*X2)+X3[1]*(X3[0]^6)*(X0^2)-1*(X3[1]*(X3[0]^6)*X0*X2)-1*(X3[1]*(X3[0]^6)*(X1^2))-1*((X3[0]^7)*X0*X1))/((X3[1]^6)*(X2^2)+4*((X3[1]^5)*X3[0]*X1*X2)+2*((X3[1]^4)*(X3[0]^2)*X0*X2)+4*((X3[1]^4)*(X3[0]^2)*(X1^2))+(X3[1]^4)*(X3[0]^2)*(X2^2)+4*((X3[1]^3)*(X3[0]^3)*X0*X1)+4*((X3[1]^3)*(X3[0]^3)*X1*X2)+(X3[1]^2)*(X3[0]^4)*(X0^2)+2*((X3[1]^2)*(X3[0]^4)*X0*X2)+4*((X3[1]^2)*(X3[0]^4)*(X1^2))+4*(X3[1]*(X3[0]^5)*X0*X1)+(X3[0]^6)*(X0^2))

Step 7: alpha = (r[i]*r[i]) / (p[i]*A[i][j]*p[j])
  alpha = 

Step 8: r[i] = r - alpha * A[i][j] * p[j]
  r[0] = 0

  r[1] = 0

Step 9: x[i] = x + alpha*p[i]
  x[0] =

  x[1] = 

END