Re: NIntegrate mit einer compelierten Funktion

[rolfm <rolfm@XXXXXXX.com>]

Subsidencefunction$internal$xy$ohne$if$compiled = Compile[{{a, _Real}, {b,
_Real}, {l, _Real}, {M, _Real}, {x, _Real}, {y, _Real}, 
     {H, _Real}}, ((-(1/4))*(Erf[((l - 2*y)*a)/(1.4142135623730951*(-1 +
H/M)^(1/2))] + 
       Erf[((l + 2*y)*a)/(1.4142135623730951*(-1 + H/M)^(1/2))])*(Erf[((b +
2*x)*a)/(1.4142135623730951*(-1 + H/M)^(1/2))] + 
       Erf[((b - 2*x)*a)/(1.4142135623730951*(-1 + H/M)^(1/2))]))/(-1 +
E^(2*M*(-H + M))), CompilationTarget -> "C"]; 
suf[(z__)?NumberQ] := Subsidencefunction$internal$xy$ohne$if$compiled[z]; 
NIntegrate[suf[1, 2, 3, 4, x, 6, 7], {x, 0, 10}, AccuracyGoal -> 8]

funktioniert.

Man sollte meinen dass man so eine Hilfsfunktion (suf) nicht braucht, aber
in dem Fall wohl doch.

Gruss,

Rolf Mertig
GluonVision GmbH
Berlin