n=input('Inserisci la dimensione della matrice ') ;
A=zeros(n) ;
B=A ;

risp=input('Vuoi inserire gli elementi della matrice? (s/n) ','s') ;

if risp=='s'
    disp('Inserisci A riga per riga')
    for i=1:n
        A(i,:)=input('Valori di A \n') ;
    end
else
    disp('A(i,j)=(1-i/2)/sqrt((j+2))')
    for i=1:n
        for j=i:n
           A(i,j)=(1-i/2)/sqrt((j+2)) ;
        end
    end
end
B(:,1:2)=A(:,1:2) ;
B(:,3:n)=2*ones(n,n-2) ;

if n<=7
    disp(A)
    disp(B)
end
s=sum(diag(B)) ;
disp(s)
