eig_power.m
33bba04159d6master

eig_power.m
View Options

			function [ vec, val ] = eig_power(inputmatrix)
			[Ni, Nj] = size(inputmatrix);
			assert(Ni == Nj \|\| Ni == 0, 'Matrix must be square or non-zero')
			vec = zeros(Ni, 1);
			vec(1) = 1;
			val = ctranspose(vec) * inputmatrix * vec;
			oldval = 0;

			while abs(oldval - val) > eps
			vec = inputmatrix * vec;
			vec = vec / norm(vec);
			oldval = val;
			val = ctranspose(vec) * inputmatrix * vec;
			end
			end