sturm.m

various matlab programs to slove various problems

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