You can use cumsum
.
You probably want to take the absolute value of the eigenvalues before summming:
P = 90; %// desired percentage
eigenValues = eigs(Matrix);
summation = cumsum(abs(eigenValues));
summation = summation/summation(end)*100; %// normalize to get percentage
k = find(summation>=P,1); %// first index for whith P is exceeded