Compute the infinite sum of pi/4 up to 5 correct decimals.

1 view (last 30 days)
I want to approximate pi/4 by using the infinite sum pi/4=((-1)^(n))/(2*n+1) from 0 to infinity. Here is my MATLAB code:
if true
tol = 10^(-5) %Up to five correct decimals
d=1; %Arbitrarily chosen, only condition is that d>tol at the start
sn=0; %Initial starting value
n=0;
while d>tol
sn = ((-1)^(n))/(2*n+1); This is the approximation for pi/4.
d=abs(pi/4-sn); Value that decides whether the loop continues.
n=sn; I'm not sure but I want n to increase by one each
time to prevent
circular calculations.
end
Why do I get Inf+ Nani?
  2 Comments
dpb
dpb on 19 Sep 2013
Why do I get Inf+ Nani?
Because of
n=sn;
Try each step directly at the command line and see what happens...
You need to set a value for the total of the summation of each term initially to zero and then accumulate the terms into that variable -- sn is ok for the name but you don't accumulate a sum of the terms in n, you replaced in with the individual term going forward.
To get the subsequent terms you need to increment n, that is correct --
n=n+1;
where the n=sn; line is instead; the implementation of accumulating the summation is left for you to consider how to do that a little more...

Sign in to comment.

Accepted Answer

Azzi Abdelmalek
Azzi Abdelmalek on 19 Sep 2013
Edited: Azzi Abdelmalek on 19 Sep 2013
tol = 10^(-5)
d=1;
sn=0; %
n=0;
while d>tol
sn = sn+(-1)^n/(2*n+1);
n=n+1;
d=abs(pi/4-sn);
end
disp(sn)
  2 Comments
Jan
Jan on 19 Sep 2013
Edited: Jan on 19 Sep 2013
This is obviously a homework question. Solving it does not allow the OP to find the solution by his own.
Azzi, please do not post solutions of homework questions, because this reduces the reputation and efficiency of the forum. As a teacher on a university you should think of your colleagues, who do not want to get solutions created by others than their students.

Sign in to comment.

More Answers (0)

Categories

Find more on Loops and Conditional Statements in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!