sturm.m

来自「various matlab programs to slove various」· M 代码 · 共 8 行

function [p]=sturm(dd,bb,x)
n=length(dd); 
p(1)=1;  
p(2)=d(1)-x;
for i=2:n
   p(i+1)=(dd(i)-x)*p(i)-bb(i-1)^2*p(i-1); 
end
return