clear all
close all
clc

codigo_tipo = input('Tipo: 1-SIMPLES; 2-COMPOSTO; 3-DUPLO (Default = 1)  => ');
if isempty(codigo_tipo)
    codigo_tipo = 1;
end
switch codigo_tipo
    case 1
        tipo='SIMPLES';
    case 2
        tipo='COMPOSTO';
    case 3
        tipo='DUPLO';
end
str = sprintf('\nTipo selecionado: %d - %s\n',codigo_tipo,tipo);
disp(str);
switch codigo_tipo
    case 1 % SIMPLES
        Ns=29;
        Np=29;
        Nr=Ns+2*Np;
        AsAc=2*(Ns+Np)/Ns;
        ArAc=2*(Ns+Np)/Nr;
        AsAr=-Nr/Ns;
        fprintf('As/Ac = %7.3f        Ac/As = %7.3f\n',AsAc,1/AsAc)
        fprintf('Ar/Ac = %7.3f        Ac/Ar = %7.3f\n',ArAc,1/ArAc)
        fprintf('As/Ar = %7.3f        Ar/As = %7.3f\n',AsAr,1/AsAr)
    case 2 % COMPOSTO
        Np1=47;
        Np2=25;
        Ns=31;
        Nr=102;
        AsAc=(Np1*Nr+Np2*Ns)/(Np2*Ns);
        ArAc=(Np1*Nr+Np2*Ns)/(Np1*Nr);
        AsAr=-Np1*Nr/(Np2*Ns);
        fprintf('As/Ac = %7.3f        Ac/As = %7.3f\n',AsAc,1/AsAc)
        fprintf('Ar/Ac = %7.3f        Ac/Ar = %7.3f\n',ArAc,1/ArAc)
        fprintf('As/Ar = %7.3f        Ar/As = %7.3f\n',AsAr,1/AsAr)
    case 3 % DUPLO
        Np1=13;
        Np2=31;
        Ns1=30;
        Ns2=12;
%         Np1=53;
%         Np2=67;
%         Ns1=42;
%         Ns2=22;
        As1Ac=(Np2*Ns1-Np1*Ns2)/(Np2*Ns1);
        As2Ac=-(Np2*Ns1-Np1*Ns2)/(Np1*Ns2);
        As1As2=Np1*Ns2/(Np2*Ns1);
        soma=(Np1+Ns1)-(Np2+Ns2);
        fprintf('As1/Ac  = %7.3f        Ac/As1  = %7.3f\n',As1Ac,1/As1Ac)
        fprintf('As2/Ac  = %7.3f        Ac/As2  = %7.3f\n',As2Ac,1/As2Ac)
        fprintf('As1/As2 = %7.3f        As2/As1 = %7.3f\n',As1As2,1/As1As2)
        fprintf('(Np1+Ns1)-(Np2+Ns2) = %8.3f\n',soma)
end






