function [m,b] = num_real_eigenvalues(A)
%Computes the eigenvalues of a matrix and returns the number
%of real ones
b=eig(A);
m=length(find(abs(imag(b))<1e-6));
end


