Python. I have a list of codes need to be debuged. The topic is using the BiCGSTAB to solve linear equation(numerical method problem). See photo one is the question. The second photo is the idea.

Here is my codes but I am not sure where is wrong.

from math import * import numpy as np import matplotlib.pyplot as plt

A = np.array([[1.0,6.0,3.0],[23.0,1.0,80.0],[7.0,8.0,91.0]]) b = np.array([[4.0],[1.0],[6.0]])

kmax = 0

 def bicgstab(x0, b = b, A = A, eps, kmax):  n = b.shape[0]  errtol = eps*np.linalg.norm(b)  error = []  x = x0    rho = np.zeros(kmax+1)    r = b - np.dot(A,x)  rhat0 = r  rhat0T = np.array(zip(*rhat0))   #return x, rho, r  return np.dot(rhat0T,r)

do different states or regions of interest. Note for this project you will also need to use a quadrature rule. Solving linear systems. Gaussian elimination was discussed as a method for solving linear systems. There are other methods that may be used for solving linear system for example iterative methods may also be considered. Iterative methods are an alternative to Gaussian elimination that focus on matrix multiplication to approximately solve linear systems. Two groups could do this project. One group researching the GMRES algorithm for solving a linear system. A second group could discus BiCGSTAB as an algorithm for solving linear systems. Implementing a GMRES or BiCGSTAB solver. Bunnies and Wolves (or the Princess and Frog). Model the population of a predator interacting with a prey population. There are several Ordinary differential equatio that can be use to do this using Runge-Kutta methods. A this modeling should look up the Lotka-Volterra model and see me for more details. ns group interested in doing do different states or regions of interest. Note for this project you will also need to use a quadrature rule. Solving linear systems. Gaussian elimination was discussed as a method for solving linear systems. There are other methods that may be used for solving linear system for example iterative methods may also be considered. Iterative methods are an alternative to Gaussian elimination that focus on matrix multiplication to approximately solve linear systems. Two groups could do this project. One group researching the GMRES algorithm for solving a linear system. A second group could discus BiCGSTAB as an algorithm for solving linear systems. Implementing a GMRES or BiCGSTAB solver. Bunnies and Wolves (or the Princess and Frog). Model the population of a predator interacting with a prey population. There are several Ordinary differential equatio that can be use to do this using Runge-Kutta methods. A this modeling should look up the Lotka-Volterra model and see me for more details. ns group interested in doing