New guy and programming problem

08282017, 10:55 AM
Post: #21




RE: New guy and programming problem  
08282017, 02:20 PM
Post: #22




RE: New guy and programming problem
So for completeness' sake, I'll add the response I left to you when you asked this question on my 42s programming video.
Sure, so what's happening is that if you enter a new value for R1 or R2 the stack is lifting and the old value is now in the Y portion of the stack, so ultimately your program is calculating the old R2 + 1/(new R2) and then 1/x. A couple of solutions: the most straightforward is probably to store the intermediate value, so after the first 1/X, STO 00 and then after the second 1/X, RCL 00. Or, perhaps a better practice is to create the variables first, then perform the operations: Code: LABEL "PAR2" Lastly, if you wanted to be really streamlined, you could do away with INPUT altogether and just enter your values into the stack before executing the program Code: LABEL "PAR2" 

08282017, 06:11 PM
Post: #23




RE: New guy and programming problem
(08282017 02:20 PM)Logan Wrote: So for completeness' sake, I'll add the response I left to you when you asked this question on my 42s programming video. Yes, I did read your reply, also lemme use the instance to thank you so much for your videos, I've watched every single one of them (calculators related) and its partially your fault that I'm so inclined to the 42S hahaha. Greetings 

08282017, 06:13 PM
(This post was last modified: 08282017 06:14 PM by Dieter.)
Post: #24




RE: New guy and programming problem
(08272017 11:12 PM)Paul Dale Wrote: Yes to your other suggestions for improving accuracy, I got in what I could  limited flash space and running out of energy justifying everything. Life is short and RAM is full. ;) (08272017 11:12 PM)Paul Dale Wrote: Apart from 1  sin(x) which doesn't have accuracy concerns that I can see Well, for 90±5E–7° even on the 34s sin(x) rounds to 1 (in SP mode). On a 10digit machine this already happens for 89,9995°...90,0005°. Even for 89°...91° three significant digits of 1 – sin(x) are lost. (08272017 11:12 PM)Paul Dale Wrote: ...and √(x²+y²) which is present already as complex ABS. Ah, yes, that's something I forgot. I even use complex ABS for this purpose here and there, both on the 34s and the 35s. Dieter 

08282017, 06:13 PM
Post: #25




RE: New guy and programming problem
Ha! I'll accept that blame!


08292017, 04:22 AM
Post: #26




RE: New guy and programming problem
(08282017 06:13 PM)Dieter Wrote: Well, for 90±5E–7° even on the 34s sin(x) rounds to 1 (in SP mode). Good point, I was wrong. For some reason I thought the flattest gradient was 1 not 0. Pauli 

09012018, 12:23 PM
Post: #27




RE: New guy and programming problem
(08282017 06:13 PM)Dieter Wrote: On a 10digit machine this already happens for 89,9995°...90,0005°. We can use this differencetoproduct formula to calculate e.g. \(1\sin(89.9995°)=\sin(90°)\sin(89.9995°)\): \(\sin \theta  \sin \varphi =2\sin \left({\frac {\theta  \varphi }{2}}\right)\cos \left({\frac {\theta + \varphi }{2}}\right)\) And then since \(\theta = \frac{\pi}{2}\) we can set \(\varepsilon = \frac {\theta  \varphi }{2}\) and use \(\cos \left({\tfrac {\pi }{2}}\varepsilon \right)=\sin \varepsilon\) and end up with: \(1  \sin \varphi =2\sin^2 \varepsilon\) Thus we set \(\theta=90\) and \(\varphi=89.9995\) and get: 90 ENTER 89.9995  2 ÷ SIN x² 2 × 3.807717748e11 (Calculated using an HP41.) Cheers Thomas 

09012018, 02:16 PM
Post: #28




RE: New guy and programming problem
(09012018 12:23 PM)Thomas Klemm Wrote:(08282017 06:13 PM)Dieter Wrote: On a 10digit machine this already happens for 89,9995°...90,0005°. Another way, rearrange 1  sin(x), so more accurate for x close to right angle. 1  sin(x) = (1  sin(x)) * (1 + sin(x)) / (1 + sin(x)) = cos(x)^2 / (1 + sin(x)) 1  sin(89.9995°) ~ cos(89.9995°)^2 / 2 ~ 3.807717747e11 

09012018, 03:18 PM
(This post was last modified: 09012018 05:23 PM by sasa.)
Post: #29




RE: New guy and programming problem
(09012018 12:23 PM)Thomas Klemm Wrote: We can use this differencetoproduct formula to calculate... Even result is good enough, there is still too much trouble to get accurate result to the last significant digit with arguments near 90 degree. The first problem is what is each calculator's internal mantissa length, the second exact representation for PI in internal format (internal conversion from degree to radian) and proper round off mechanism. And as well the used format (decimal or binary floating point) and proper algorithms to calculate trig functions. The Wolfram Alpha returns: 3.8077177473096929677620869980236729697355177044497495... × 10^11 

09012018, 04:28 PM
(This post was last modified: 09012018 04:35 PM by Thomas Klemm.)
Post: #30




RE: New guy and programming problem
(09012018 03:18 PM)sasa Wrote: Even result is good enough, there is still too much trouble to get accurate result to the last significant digit with arguments near 90 degree. For sin(0.00025°) WolframAlpha gives: 4.3633231299719786898764073194250596667396107572873548... × 10^6 The HP41C returns: 4.363323130e06 This result is correctly rounded to the last digit. Then we square the rounded number and get: 1.90385887367929969 × 10^11 The HP41C returns: 1.903858874e11 Again the result is rounded correctly. Last thing to do is to multiply the number by 2: 3.807717748 × 10^11 Which is exactly what the HP41C returns: 3.807717748e11 As we can see the HP41C rounds these intermediate results correctly. The problem with the cancelation when calculating \(1sin(x)\) is not related to the accuracy of calculating \(sin(x)\) for \(x\) close to \(90°\). Similar to Albert's trick we can multiply \(\sqrt{1+x^2}1\) by \(\frac{\sqrt{1+x^2}+1}{\sqrt{1+x^2}+1}=1\) and get: \(\begin{align*} \sqrt{1+x^2}1 &= \sqrt{1+x^2}1 \cdot \frac{\sqrt{1+x^2}+1}{\sqrt{1+x^2}+1} \\ &= \frac{(\sqrt{1+x^2}1)(\sqrt{1+x^2}+1)}{\sqrt{1+x^2}+1} \\ &= \frac{1+x^21}{\sqrt{1+x^2}+1} \\ &= \frac{x^2}{\sqrt{1+x^2}+1} \end{align*}\) And again cancelation for small \(x\) is avoided. 

09012018, 05:20 PM
Post: #31




RE: New guy and programming problem
(09012018 03:18 PM)sasa Wrote: Even result is good enough, there is still too much trouble to get accurate result You might have the impression that Wolfram answer is easily calculated. It is *not* You don't feel the difficulty, because modern computers are just too fast. Or, in the old days, where you actually get *charge* for computing time If you don't believe me, try to reproduce Wolfram answer by hand. 

09012018, 05:33 PM
(This post was last modified: 09022018 08:31 AM by sasa.)
Post: #32




RE: New guy and programming problem
(09012018 05:20 PM)Albert Chan Wrote: You might have the impression that Wolfram answer is easily calculated. It is *not* Yes, I have had to go myself through all these headache with my own library for MCUs... Simply, there is too many things to balance and decide proper path, if resources are quite limited. 

09012018, 07:28 PM
(This post was last modified: 09012018 07:33 PM by Vtile.)
Post: #33




RE: New guy and programming problem
(08272017 05:12 PM)Gerson W. Barbosa Wrote: Is there another calculator that offers this as a builtin function? Yes and it is even a standard factory build: Casio FX61F. What comes to the OPs original question.. Dieters two last program snippets are the way it should be done with the RPN or RPL machine (atleast the 2ns last should also work without issues with impedances). Then there is some ways to do it "better" to avoid some insignificant rounding errors which are irrelevant most all practical levels. 

09022018, 06:14 PM
(This post was last modified: 09022018 06:40 PM by Albert Chan.)
Post: #34




RE: New guy and programming problem
Third way to calculate accurate 1  sin(x), is to remove sine function, replaced with cosine.
We can then pick whatever double angle formula is suitable: cos(2x) = cos(x)²  sin(x)² = 2 cos(x)²  1 = 1  2 sin(x)² 1  sin(89.9995°) = 1  cos(0.0005°) = 1  (1  2 sin(0.00025°)²) = 2 sin(0.00025°)² = 3.807717747e11 Note: above same as Klemm differencetoproduct formula (post 27), but probably easier to remember. BTW, you can build multiple angles cosine formula with cosines ... Neat cos(n x) = 2 cos(x) cos((n1) x)  cos((n2) x) 

09022018, 07:33 PM
Post: #35




RE: New guy and programming problem
Fun fact: These trigonometric identities were used for multiplication and division before logarithms were known.
Cf. Prosthaphaeresis Quote:Notice how similar the above algorithm is to the process for multiplying using logarithms, which follows these steps: scale down, take logarithms, add, take inverse logarithm, scale up. It's no surprise that the originators of logarithms had used prosthaphaeresis. Indeed the two are closely related mathematically. In modern terms, prosthaphaeresis can be viewed as relying on the logarithm of complex numbers, in particular on Euler's formula: I find it interesting that this technique was developed using complex numbers under the hood while they weren't discovered yet. Cheers Thomas 

09022018, 09:21 PM
(This post was last modified: 09022018 09:22 PM by Albert Chan.)
Post: #36




RE: New guy and programming problem
I never knew trig function were used as logarithm !
Example: arccos/cos vs ln/exp 0.123 * 0.456 = cos(82.9347°) * cos(62.8707°) = 1/2 (cos(82.9347° + 62.8707°) + cos(82.9347°  62.8707°)) = 1/2 (cos(145.8054°) + cos(20.0640°)) = 0.05608822 ~ 0.056088 0.123 * 0.456 = exp(ln(0.123) + ln(0.456)) = exp(2.09557 + 0.78526) = exp(2.88083) = 0.05608819 ~ 0.056088 But, how can trig function used for division ? 

09022018, 09:33 PM
Post: #37




RE: New guy and programming problem  
09022018, 09:40 PM
Post: #38




RE: New guy and programming problem  
09032018, 10:32 PM
Post: #39




RE: New guy and programming problem  
« Next Oldest  Next Newest »

User(s) browsing this thread: