Hallo,
als Anlage übersende ich Euch ein kleines Notebook. Darin ist die Definition
für eine Funktion enthalten, die aus einigen Erf und einer Gausschen
Glockenkurven bzw. einer inversen Gausschen Glockenkurven zusammengesetzt
ist:
\!\(\((\[ExponentialE]\^\(1\/4\ k\ \((h - T)\)\^2\ \((\(-1\) + T\/h)\)\^m\)\
\
\((Erf[1\/4\ a\ \@k\ \((\(-1\) + T\/h)\)\^\(m/2\)\ \((l -
2\ r\ Cos[\[Alpha]])\)] +
Erf[1\/4\ a\ \@k\ \((\(-1\) + T\/h)\)\^\(m/2\)\ \((l +
2\ r\ Cos[\[Alpha]])\)])\)\ \((Erf[
1\/4\ a\ \@k\ \((\(-1\) + T\/h)\)\^\(m/2\)\ \((b -
2\ r\ Sin[\[Alpha]])\)] +
Erf[1\/4\ a\ \@k\ \((\(-1\) + T\/h)\)\^\(m/2\)\ \((b +
2\ r\ Sin[\[Alpha]])\)])\))\)/\((4\ \((\(-1\) + \
\[ExponentialE]\^\(1\/4\ k\ \((h - T)\)\^2\ \((\(-1\) + T\/h)\)\^m\))\))\)\)
Ich will diese nun numerisch über h integrieren.
Numerisch, weil leider ein kein analytisches Integral nach h existiert.
Numerisch auch, weil eine Serienbildung nicht trivial ist (es entstehen
immer unbestimmte Ausdrücke, die ich auch über die üblichen Tricks, wie
Limit, nicht umgehen konnte).
NIntegrate hat aber auch öfters Probleme.
Den Grund kann man erkennen, wenn man die Funktion über h aufträgt (siehe
Notebook). Bei sehr kleinen h fängt die Funktion an zu flattern.
Hat jemand einen Tipp, wie man das umschiffen kann?
Gruß, Peter Klamser
Flattern.zip
Description: Zip compressed data
