%akolouthia Sturm gia ton upologismo twn idiotimwn enos mitrwou
function solution=SturmS(A)
[n,n] = size(A);
ff(1,:)=[1 0 0 0 0]; 
ff(2,:)=[A(1,1) -1 0 0 0];
for i=3:n+1
h=[A(i-1,i-1) -1];
ff(i, 1) = h(1)*ff(i-1,1)-A(i-1, i-2)^2*ff(i - 2, 1);
for z=2:n+1
ff(i, z) = h(1)*ff(i-1, z)+h(2)*ff(i-1, z-1)-A(i - 1, i - 2)^2*ff(i - 2, z); 
end; 
end;
for i=1:n+1
y(i)=ff(n+1,n+2-i);
end;
eigval=roots(y)
end