function [z,theta]=compito(n,alpha)
beta=2+alpha;  % 1
z=zeros(n-1,1);
z(n-1)=beta;
z(n-2)=beta*z(n-1)-1;  %2
for i=1:n-3
    z(n-(i+2))=beta*z(n-(i+1))-z(n-i);  %  2 op per ogni iterazione
end
theta=beta*z(1)-z(2);  % 2 op