PetscInitialize({'-pc_type','ml','-ksp_monitor'petsc-3.5.4 2015-05-23}); v = PetscViewer('/Users/barrysmith/Datafiles/Matrices/medium',Petsc.FILE_MODE_READ); A = PetscMat; A.Load(v); v.Destroy; ksp = PetscKSP; ksp.SetOperators(A); ksp.SetFromOptions; b = PetscVec(1:181); x = PetscVec(1:181); ksp.Solve(b,x); ksp.Destroy; A.Destroy; PetscFinalize;