Hallo,
Erf[x>6]=1.
Leider spinnt MMA manchmal, wenn für sehr große Zahlen die Erf berechnet
werden soll:
Plot[Erf[n], {n, -10, 100000}, Frame -> True, PlotRange -> All]
Dann hilft
Errorfunktion[(z_)?NumericQ] := If[Abs[z] > 6, Sign[z], Erf[z]];
Errorfunktion[(z0_)?NumericQ, (z1_)?NumericQ] := Errorfunktion[z1] -
Errorfunktion[z0];
Plot[{Erf[n], 1.025*Errorfunktion[n]}, {n, -6, 6}, Frame -> True,
PlotRange -> All]
Plot[{Errorfunktion[n]}, {n, -6, 60000}, Frame -> True, PlotRange -> All]
Diese "eigene" Errorfunktion würde ich gerne im Rahmen einer
Unprotect[Erf] Protect[Erf] Definition in die Erf einbauen, was aber in
einer unendliche Iteration endet:
Unprotect[Erf]
Erf[(z_)?NumericQ] := If[Abs[z] > 6, Sign[z], Erf[z]];
Protect[Erf]
Erf[8]
In[45]:=
N[Erf[2]]
$IterationLimit::"itlim":"Iteration limit of \!\(4096\) exceeded.
\!\(\*ButtonBox[\"Mehr\[Ellipsis]\", ButtonStyle->\"RefGuideLinkText\",
ButtonFrame->None, ButtonData:>\"$IterationLimit::itlim\"]\)"
Hold[Erf[2.]]
Was mach ich Falsch?
Gruß
Peter
|