Finding an approximation for cos(x) using for loops

Question

Laura 2013-11-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/106482-finding-an-approximation-for-cos-x-using-for-loops

评论： Weeraphol Egsaphung 2017-9-23

采纳的回答： Umair Nadeem

在 MATLAB Online 中打开

Hi there,

I recently got help on creating a script to calculate the value of cos(x) using for loops. The question asks:

"Write a script which will calculate the value of cos(x) (ask the user for the value of x and n):

cos(x)=( from k=1 –> n )∑(-1) ̂(k-1)(x ̂(2(k-1))/(2(k-1))!)

" I had real problems getting the answer and so asked a tutor to help me with it. She managed to get the code working but I failed to understand how it actually works. (her explanation was not very clear)

So, I have sat staring at it for ages and it still doesn't make sense to me.

Here it is:

% Ask user for values x, n
x=input('Please enter a value for x: '); %cos angle
x_rad=x*pi/180;
n=input('Please enter a value for n: '); %nth term
%factorial
fact=1;
sum=0;
for k=1:n
  initial=-1;
  for a=0:(k-1)
    initial=(-1)*initial;
  end
  numerator=1;
  for b=1:2*(k-1)
    numerator=numerator*x_rad;
  end
  for t=1:2*(k-1)
    fact=fact*t;
  end
  total=(initial*numerator)/fact;
  sum=sum+total;
end
fprintf('cos(%d)=%6.6f \n',x,sum);

Can anyone explain how this code is working in the for loop? I am so confused.

Also, any advice on a simpler way to do this?

Thanks so much, Laura

2 个评论
显示无隐藏无

Walter Roberson 2013-11-17

Which "for" loop ?

Laura 2013-11-17

All of them. I don't understand why we need nested for loops in the first place and why the end values for each one is different. (e.g. for b=1:2*(k-1))

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Umair Nadeem 2013-11-17

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/106482-finding-an-approximation-for-cos-x-using-for-loops#answer_115501

编辑：Umair Nadeem 2013-11-17

在 MATLAB Online 中打开

The code has been developed using a more complex approach. You can simply implement the equation you gave above. Here is how to do it.

% Ask user for values x, n
x=input('Please enter a value for x: '); %cos angle
x_rad=x*pi/180;
n=input('Please enter a value for n: '); %nth term
% Initialize resultant variable
sum=0;
for k=1:n
  initial = (-1)^(k-1);
  numerator = x_rad^(2*(k-1));
        denominator = factorial(2*(k-1));
  total=(initial*numerator)/denominator;
  sum=sum+total;
end
fprintf('cos(%d)=%6.6f \n',x,sum);