Simulasi Barisan dan Deret Menggunakan Matlab

clc;clear all;
disp('Program   : Simulasi Barisan Deret');
disp('Programer : Syaharuddin');
disp('Pilihan:');
disp('1. Barisan Deret Aritmatika');
disp('2. Barisan Deret Geometri');
p=input('Pilihan = ');
disp('==============================================');
switch p
    case 1
        a1=input('Suku ke-1 = ');
        a2=input('Suku ke-2 = ');
        n=input('n = ');
        b=a2-a1;
        disp(['Beda = ',num2str(b)]);
        if n<1
            disp('Tidak Terdefinisi');
        else
        t=[a1:b:n*b];
        Un=a1+(n-1)*b;
        Sn=(n/2)*(a1+Un);
        disp(['Barisan: ',num2str(t)]);
        disp(['Jadi, Un = ',num2str(a1),' + (',num2str(n),'-1) x ',num2str(b),' = ',num2str(Un)]);
        disp(['Jadi, Sn = (',num2str(n),'/2) x (',num2str(a1),' + ',num2str(Un),') = ',num2str(Sn)]);
        end
    case 2
        a1=input('Suku ke-1 = ');
        a2=input('Suku ke-2 = ');
        n=input('n = ');
        r=a2/a1;
        disp(['Rasio = ',num2str(r)]);
        t=[a1:r:n*r];
        Un=a1*(r^(n-1));
        disp(['Barisan: ',num2str(t)]);
        disp(['Jadi, Un = ',num2str(a1),' x (',num2str(r),'^(',num2str(n),'-1) = ',num2str(Un)]);
        if r>1
        Sn=(a1*((r^n)-1))/(r-1);
        disp(['Jadi, Sn = [',num2str(a1),' x (',num2str(r),'^',num2str(n),'-1)]/',num2str(r),'-1 = ',num2str(Sn)]);
        else
          r<1
        Sn=(a1*(1-(r^n)))/(1-r);
        disp(['Jadi, Sn = [',num2str(a1),' x (1-',num2str(r),'^',num2str(n),')]/1-',num2str(r),' = ',num2str(Sn)]);
        end
end