lcd.m

来自「一本很简单」· M 代码 · 共 15 行

function d = LCD(m)

% d = LCD(m) returns the least common denominator of the entries of m.

[num,denom] = rat(m, 1e-10);
[m, n] = size(denom);
d = 1;
for i = 1:m
  for j = 1:n
    g = gcd(d, denom(i,j));
    d = d*denom(i,j)/g;
  end;
end;
return;