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(b==real(b)));
end

n=1000;
m;
vec=zeros(m);
for