BiCGSTAB Fortran solver using Biconjugate Gradient Stabilized (BiCGStab) method in order to solve large systems of equations