Autor |
Rundungsfehler in der Praxis |
|
Wally
Senior  Dabei seit: 02.11.2004 Mitteilungen: 9811
Wohnort: Dortmund, Old Europe
 | Themenstart: 2019-12-11
|
Als ich an Squires zweitem Rätsel ein wenig herumprobierte, habe ich eine seltsame Entdeckung gemacht- ein Integral einer positiven Funktion wird negativ ;)
\sourceon Maple
evalf(int(21*x^21/(x^2+4*x+7), x = 0 .. 1));
-8.
\sourceoff
Wolfram Alpha gibt folgendes witziges Ergebnis: https://www.matheplanet.de/matheplanet/nuke/html/uploads/b/6234_mp.png
So ein ganz klein bischen ist natürliche Intelligenz noch besser als künstliche - tröstlich!
Wally
|
Profil
|
Caban
Senior  Dabei seit: 06.09.2018 Mitteilungen: 3264
Wohnort: Brennpunkt einer Parabel
 | Beitrag No.1, eingetragen 2019-12-11
|
Hallo
Da kenne ich auch ein Beispiel: Ein Taschrechner hat als Nullstelle von x*exp(x) -527 und viele andere solche Werte ausgegeben. Der Schüler hat das aber ggeglaubt.
Gruß Caban
|
Profil
|
markusv
Senior  Dabei seit: 24.01.2017 Mitteilungen: 339
Wohnort: Leipzig
 | Beitrag No.2, eingetragen 2019-12-11
|
Der "praktische" Bezug dazu fehlt mir zwar, weshalb ich den Titel nicht verstehe, aber Wolfram Alpha gibt durch Klicken auf "more digits" 0,0813 als Ergebnis zur Schau. Hätte zumindest meine natürliche Intelligenz jetzt nicht so schnell ermittelt.
|
Profil
|
Kitaktus
Senior  Dabei seit: 11.09.2008 Mitteilungen: 7276
Wohnort: Niedersachsen
 | Beitrag No.3, eingetragen 2019-12-11
|
Aber warum wird da im Ergebnis überhaupt ein Faktor 105 abgespaltet? Wenn da 8.*10-2 stünde, wäre das ja verständlich.
Wahrscheinlich sind die einzelnen Summanden betragsmäßig sehr groß, löschen sich aber im Wesentlichen gegenseitig aus.
|
Profil
|
wessi90
Senior  Dabei seit: 16.09.2011 Mitteilungen: 2127
 | Beitrag No.4, eingetragen 2019-12-11
|
Mathematica gibt mir nach ca. 2.2s Rechenzeit das Ergebnis 0.0813247 aus, warum Maple hier Murks veranstaltet weiß ich aber auch nicht.
|
Profil
|
Diophant
Senior  Dabei seit: 18.01.2019 Mitteilungen: 11121
Wohnort: Rosenfeld, BW
 | Beitrag No.5, eingetragen 2019-12-11
|
\(\begingroup\)\(\newcommand{\ba}{\begin{aligned}}
\newcommand{\ea}{\end{aligned}}
\newcommand{\bpm}{\begin{pmatrix}}
\newcommand{\epm}{\end{pmatrix}}
\newcommand{\bc}{\begin{cases}}
\newcommand{\ec}{\end{cases}}
\newcommand{\on}{\operatorname}\)
Hallo zusammen,
dann ist Mathcad Prime doch nicht so schlecht wie sein Ruf: ich bekomme den Wert \(0.0813242161\) unmittelbar, ohne irgendeine Wartezeit. Mit Maxima 5.41 bekomme ich noch mehr Stellen, ebenfalls unmittelbar ohne spürbare Wartezeit.
@wessi90:
Warum weicht das Ergebnis von Mathematica hier ab der 7. Nachkommastelle ab, oder ist das ein Tippfehler?
Gruß, Diophant \(\endgroup\)
|
Profil
|
wessi90
Senior  Dabei seit: 16.09.2011 Mitteilungen: 2127
 | Beitrag No.6, eingetragen 2019-12-11
|
Nein, das war kein Tippfehler. Ich habe für die erste Eingabe
\sourceon Mathematica
N[Integrate[]]
\sourceoff
benutzt, was lange dauert, da er zunächst exakt integriert und dann den numerischen Wert angibt. Mit NIntegrate, also numerischer Integration, kommt der korrekte Wert sofort ohne merkbare Wartezeit.
|
Profil
|
ThomasRichard
Senior  Dabei seit: 08.04.2010 Mitteilungen: 476
Wohnort: Aachen
 | Beitrag No.7, eingetragen 2019-12-11
|
Das ist kein Murks, sondern ein ganz normales Phänomen in jedem CA-System: der exakte Ausdruck, den Wally mit int ermittelt und dann direkt an evalf weitergereicht hat, erfordert eine etwas höhere als die Default-Genauigkeit von 10 Stellen. Einfach Digits:=15 setzen, dann liefert Maple eine gute Näherung - und benutzt noch die Hardware-Arithmetik, also ohne wesentliche Geschwindigkeitseinbuße. Ab Digits:=16 werden (beliebig genaue) Software Floats verwendet.
Wenn man nur am numerischen Ergebnis interessiert ist, sollte man direkte numerische Integration verwenden, also entweder int mit der Option 'numeric' oder die "inerte" Form Int in evalf einwickeln. Dann überspringt Maple die exakte Integration (die ja eh nicht immer möglich ist).
|
Profil
|
Diophant
Senior  Dabei seit: 18.01.2019 Mitteilungen: 11121
Wohnort: Rosenfeld, BW
 | Beitrag No.8, eingetragen 2019-12-11
|
Hallo ThomasRichard,
Dankeschön für die Erläuterung:
\quoteon(2019-12-11 15:13 - ThomasRichard in Beitrag No. 7)
Wenn man nur am numerischen Ergebnis interessiert ist, sollte man direkte numerische Integration verwenden, also entweder int mit der Option 'numeric' oder die "inerte" Form Int in evalf einwickeln. Dann überspringt Maple die exakte Integration (die ja eh nicht immer möglich ist).
\quoteoff
Ich verwende zwar kein Maple, sondern Mathcad Prime, und dort scheint das ähnlich zu sein, so dass ich an dieser Stelle jetzt auch wieder etwas gelernt habe: auch in Mathcad weichen die Werte zwischen der numerischen Integration und der exakten mit anschließender Auswertung voneinander ab, wenn auch ein klein wenig geringer als in Mathematica (an der gleichen Stelle, aber nur etwa halb so weit).
Ein Werkzeug zur Messung der Rechenzeit gibt es meines Wissens in Mathcad nicht, jedoch dauerte die exakte Integration bei mir dann auch so ca. 2s.
Gruß, Diophant
|
Profil
|
Caban
Senior  Dabei seit: 06.09.2018 Mitteilungen: 3264
Wohnort: Brennpunkt einer Parabel
 | Beitrag No.9, eingetragen 2019-12-11
|
Hallo
Mein GTR liefert in schneller Zwit, das Ergebnis von Diophant.
Gruß Caban
|
Profil
|