why matlab Make a mistake in calculating sin and cos while teta in wiling to zero ?

3 views (last 30 days)
dear all
as it is obvious base on the Taylor series or l'hopital's rule , answer the below equation is 1/2
Taylor series :
(1-cos(teta))/(teta^2)= 1/2! - (teata^2)/4!
and while teta wiling to zero others terms become zero and the answer is 1/2
l'hopital's rule
(1-cos(teta))/(teta^2)= sin(teta)/(2*teta) =(1/2)*(cos(teta))
and again while teta wiling to zero answer is 1/2
but matlab act different and it gives right answer when you assume that teta is 0.01 or 0.001 but when you assume that teta is 0.00001 matlab give zero answer to this equation
is it my fault or something wrong with matlab and i really appreciated if some one can help me and if let me know if there is code for my problem

Accepted Answer

Walter Roberson
Walter Roberson on 26 Jan 2018
By teta = .0000001 then cos(teta) becomes 1 to the limit of double precision, considering round-off error. Then 1-cos(teta) is 0, so you get 0 as the result.
This is a limitation in numeric processing with finite precision. If double precision had (for example) 80 bits instead of 53, then the problem would still occur, just with a smaller x.
Walter Roberson
Walter Roberson on 26 Jan 2018
The Symbolic toolbox is the main supported tool for extended precision. You could also look at https://www.mathworks.com/matlabcentral/fileexchange/6446-multiple-precision-toolbox-for-matlab (read the comments on the submission for some corrections.)

Sign in to comment.

More Answers (1)

Torsten on 26 Jan 2018
Edited: Torsten on 26 Jan 2018
Use symbolic maths instead of numerical:
Best wishes


Community Treasure Hunt

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

Start Hunting!

Translated by