In this paper, we study conjugate gradient (CG) and precondition conjugate gradients methods (PCG), such as bi-conjugate gradient method (Bi-CG), bi-conjugate gradient stabilized method (Bi-CGSTAB), composite step bi-conjugate gradient method (CSBCG), composite step conjugate gradient squared method (CSCGS) and also the quasi-minimal residual methods (QMR) for solving linear systems. Then a new software package is provided for the mentioned methods. In addition, using a category of the family L^2 (a,b) is used to displaying integral operators as Toeplitz and sparse systems. These categories include B-spline functions and wavelet functions. We have compared the order of operations and the time limit using each of the methods mentioned. Numerical examples are also provided.
