solution.m

gmres.m为迭代算法GMRES的实现

clear all;
close all;
%function [x, error, iter, flag] = gmres( A, x, b, M, restrt, max_it, tol )
A=[1 1 1;1 2 1;2 1 1];
x=[1 0 0 ];
b=[1 3 5];