Error in a code which computes sqrt(a)
显示 更早的评论
One of the methods to compute sqer (a) , a>0 is
X(n+1) = (a + (X(n)*X(n-1))/(X(n)+X(n-1)), n = 1, 2, …, with X0=1 and X1=a (That is, it is known that lim n-> infin of Xn = sqrt(a)
Write a function [sqa, nitr] = mySqrt(a) which implements this calculation. The function should use a while-loop, terminate when the difference between Xn+1 and Xn becomes smaller than eps(10*a), and output Xn+1 in sqa and the value of n at which the while-loopwas terminated in nitr. Test your function for a = 103041
I have written this but it does not work
function [sqa, nitr] = mySqrt (a)
%[sqa, nitr] = mySqrt (a)
% computes square root of a
% sqa = a;
sqaprev = a;
nitr = 0;
X(n+1) = (a + (X(n)*X(n-1))/(X(n)+X(n-1))); %find the second term
sqa= X(n+1)
while abs (sqaprev-sqa) >= eps (10*a)
sqaprev = sqa;
sqa = (1/2) *(sqaprev+ (a/sqaprev));
nitr = nitr + 1;
end %while
end
4 个评论
Please give more information. Just saying it 'does not work' isn't much use to anyone wanting to help.
For a start though, na does not exist so that line will error immediately, irrespective of what the rest of the code may do.
This line also does not make sense, for numerous reasons:
X(n+1) = (a + (X(n)*X(n-1))/(X(n)+X(n-1)));
X also does not exist and you can't simply past an equation into Matlab and expect it to behave as code. I'm pretty sure the calculation of x(n+1) is what the while loop is supposed to be for.
Anastasia Kyriakou
2020-2-25
Anastasia Kyriakou
2020-2-25
Adam
2020-2-25
When you have a sequence in which the next term is defined in terms of the previous one you have to initialise the first term somehow, otherwise there is no previous term for the first case, to get the whole sequence going. In your case your next term is based on two previous terms, which is why you are given the first two terms of the sequence in the question so that you can then loop round increasing n and referrinig to the previous two terms, which will exist if you insert them as e.g.
X = [ 0 a ];
回答（1 个）
Samatha Aleti
2020-2-28
You may need to initialize the variables "sqaprev", "sqa" as follows:
function [sqa, nitr] = mySqrt (a)
sqaprev = 0;
sqa= a;
nitr = 0;
while abs (sqaprev-sqa) >= eps (10*a)
sqaprev = sqa;
sqa = (1/2) *(sqaprev+ (a/sqaprev));
nitr = nitr + 1;
end
end
类别
在 帮助中心 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息
2020-2-25
2020-2-28
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
