sturmcount.m

来自「麻省理工ComputationalMethod」· M 代码 · 共 10 行

function count=sturmcount(A,x)

n=size(A,1);
p=[0,1,A(1,1)-x];
for k=2:n
  p=[p,(A(k,k)-x)*p(end)-A(k,k-1)^2*p(end-1)];
end

count=sum(diff([1,sign(p(3:end))])~=0);