clearglobal

// Definição das variáveis
c = input ('Entre com o número de termos da série que se deseja visualizar:');
n=1
A = zeros (c)
S = zeros (c)
s = 0
P = zeros (c)
p = 1
Q = zeros (c)
q = 1
k = input ('Entre com um número entre -1 e 1:');

for i=1:c
    P(i) = 1/((2*i)-1)
end

for i=1:c
    Q(i) = 1/(2*i)
end

// Definição da série mágica

for n=1:c
    if s<=k then
        A(n)=P(p)
        p=p+1
    else
        A(n)=-Q(q)
        q=q+1
    end
    s=s+A(n)
    S(n)=s
end

mprintf ('Os primeiros %i termos da sequência que define a série que converge para %g são',c,k);
disp(A)
mprintf ('As primeiras %i somas parciais da série que converge pra %g são',c,k);
disp(S)
mprintf ('A soma parcial s(%i) que aproxima %g é dada por',c,k);
disp(s)
