Cos (v. 1.0.1)

function [] = calculo_cos_b(x, Epsilon)
% paulogp
% calculo de y = cos(x)
% erro < E
% funcao: calculo_cos_b(angulo, k)

    disp('-------- INICIO --------');

    syms n;
    syms Rk;
    syms aux;
    syms cos_x;

    % numeros longos
    format long

    % contador a zero
    n = 0;
    k = 0;

    % valor inicial
    Rk = 20;
    TruncError = Epsilon / 2;

    % variavel para somatorio
    cos_x = 0;

    % tabela
    disp(' Iteracao  Valor');

    % ciclo
    while Rk >= TruncError
        % formula de Rk
        Rk = (-1).^(k+1) * ( x.^(2.*k+2) / factorial(2.*(k + 1)) );

        % cos (x)
        cos_x = cos_x + ((-1).^n * (x^(2*n) / factorial(2*n)));
        aux = [n cos_x];
        disp(aux);

        n = n + 1;
        k = n;
    end

    % saida de valores
    disp('Erro de Truncatura =');
    disp(TruncError);

    disp('Rk = ');
    disp(Rk);

    disp('--------- FIM ---------');

    % apagar variaveis
    clear;
end