function [g]=g_fun(H)
%calculate the gap of the metrix H
[m]=size(H,1);
zero_length=zeros(1,m);
for i=1:m
    none_zero(i)=zero_fun(H(i,:));
    zero_length(i)=i+none_zero(i);
end;
g=m-min(zero_length);