Alternating harmonic series

Alternating harmonic series

1 - 1/2 + 1/3 - 1/4 + ...,

n = 2;

p = 6*n+1;

t = 1:1:4*n;

s = isprime(p);

f = fractional(n);

plot(t,f,'--');

title("Alternating harmonic series");

legend("H(n)");

f = sum(f);

format rational

[num, dem] = rat(f);

fprintf("The sum of the series: %d/%d \n",num,dem);

q = mod(num,p);

if((q==0)&&(s==1))

fprintf("The fraction is a multiple of the prime number %d.",p);

end

function [S] = fractional(n)

S = [];

for i = 1:4*n

S(end+1) = ((-1).^(i+1))./i;

end

end