% A sample run of ml(k)bicgstab for the solution of A x = b with
% preconditioner M.

clc
clear

N = 100;  
k = 10;
A = randn(N);
M = randn(N);
b = randn(N,1);

x = zeros(N,1);
if k == 1
    Q = b-A*x;
else
    Q = [b-A*x,  randn(N,k-1)];
end
max_it = 10*N;
tol = 10^(-7);

[x,err,iter,flag] = mlbicgstab(A,x,b,Q,M,max_it,tol);

err
iter
flag