GMRES.m

function [x]=GMRES(A,b,m,x0)
  r0=b-A*x0;
  eps=10^-6;
  I=eye(m+1);
  [V,H]=mArnoldi_GS(A,m,r0);
  Vm=V;
  [Q,R]=qr(H,0);
  [Rl,Rc]=size(R);
   Rt=zeros(Rl,Rc+1);
   Rt(1:Rl,1:Rc)=R;
   b=norm(r0)*(Q'*I(:,1));
   y=Rt\b;
   x=x0+Vm*y;