(* ::Package:: *)

(************************************************************************)
(* This file was generated automatically by the Mathematica front end.  *)
(* It contains Initialization cells from a Notebook file, which         *)
(* typically will have the same name as this file except ending in      *)
(* ".nb" instead of ".m".                                               *)
(*                                                                      *)
(* This file is intended to be loaded into the Mathematica kernel using *)
(* the package loading commands Get or Needs.  Doing so is equivalent   *)
(* to using the Evaluate Initialization Cells menu command in the front *)
(* end.                                                                 *)
(*                                                                      *)
(* DO NOT EDIT THIS FILE.  This entire file is regenerated              *)
(* automatically each time the parent Notebook file is saved in the     *)
(* Mathematica front end.  Any changes you make to this file will be    *)
(* overwritten.                                                         *)
(************************************************************************)



Arnoldi[A_,b_,nMax_]:=Module[
{m=Length[A],Q,H,n,j,h,v},
Q=ConstantArray[0,{m,nMax+1}];
H=ConstantArray[0,{nMax+1,nMax}];
Q[[All, 1]] = b/Norm[b];
Do[
v = A.Q[[All, n]];
Do[
H[[j,n]] = Q[[All, j]].v;
v = v - H[[j,n]] Q[[All, j]],
{j, 1, n}
];
H[[n+1,n]] = Norm[v];
Q[[All, n+ 1]] = v/H[[n+1,n]],
{n, 1, nMax}
];
{Q,H}
]



