Compute_Matrix.m
6e883b19735amaster

Compute_Matrix.m
View Options

			function [MOI_deg] = Compute_Matrix(A,M,deg)

			MOI_deg = 0;

			% First term
			for k = 1:deg
			MOI_deg = MOI_deg + Compute_ck(k) * (MA^k+A^kM);
			end

			% Second term
			for k1 = 1:deg-1
			for k2 = 1:deg-k1
			MOI_deg = MOI_deg - (Compute_ck(k1)Compute_ck(k2)) A^k1 * M * A^k2;
			end
			end
			end