DMUG-Archiv 2001

Frühere   Chronologischer Index   Spätere
Vorherige   Thematischer Index   Nächste

NIntegrate

Hallo,

leider habe ich eine Funktion, die sich nicht symbolisch nach dh integrieren
lässt.

\!\(\[Integral]\_0\%h\((\ \[ExponentialE]\^\(\(1.0914477214063762`\ \((h - \
Z)\)\^2\)\/\((\(-1\) + Z\/h)\)\^1.0216118780371475`\)\ \
\((Erf[\(0.5223618768168233`\ a\ \((l - 2\ r\ Cos[\[Alpha]])\)\)\/\((\(-1\)
+ \
Z\/h)\)\^0.5108059390185737`] +
                Erf[\(0.5223618768168233`\ a\ \((l + 2\ r\
Cos[\[Alpha]])\)\)\
\/\((\(-1\) + Z\/h)\)\^0.5108059390185737`])\)\ \
\((Erf[\(0.5223618768168233`\ a\ \((b - 2\ r\ Sin[\[Alpha]])\)\)\/\((\(-1\)
+ \
Z\/h)\)\^0.5108059390185737`] +
                Erf[\(0.5223618768168233`\ a\ \((b + 2\ r\
Sin[\[Alpha]])\)\)\
\/\((\(-1\) + Z\/h)\)\^0.5108059390185737`])\))\)/\((4\ \((\(-1\) + \
\[ExponentialE]\^\(\(1.0914477214063762`\ \((h - Z)\)\^2\)\/\((\(-1\) + \
Z\/h)\)\^1.0216118780371475`\))\))\) \[DifferentialD]h\)

Nicht ganz so schlimm, es gibt ja NIntegrate.

Leider läuft das nicht ganz so problemlos.

Immer wieder gibt es Fehlermeldungen, dass die Funktion oszillieren würde.

Stimmt nicht ganz. Es gibt nur für sehr kleine h numerische Instabilitäten.
Die Funktion fängt an zu flattern. Wenn ich das umschiffe bekomme ich
trotzdem Probleme, weil z.B. NIntegrate Funktionswerte außerhalb der
Integrationsgrenzen berechnet.

Hat jemand eine Idee, wie man z.B. NIntegrate entlocken kann, mit welchen
Parametern er im Augenblick die Funktion aufruft?

Wie verwendet man z.B. Trace etc.?

Gruß, Peter Klamser


Frühere   Chronologischer Index   Spätere
Vorherige   Thematischer Index   Nächste

DMUG DMUG-Archiv, http://www.mathematica.ch/archiv.html