线性方程组求解、预处理迭代法、最小残差法fortran程序