Problem mit FindRoot

[Michael Gamer <Michael.Gamer@XXXXXXX.de>]

Hallo zusammen, 

  zu Demonstrationszwecken habe ich die Berechnung des Schnittes zweier Kreise via numerischer Integration dargestellt. Dazu dient mir die Funktion A, definiert durch:

A[r_, R_, d_] := Module[{h, s = (R^2 - r^2 - d^2)/(2 d)},
  h = Abs[Sqrt[r^2 - s^2]];
  2 NIntegrate[ d + Sqrt[r^2 - x^2] - Sqrt[R^2 - x^2], {x, 0, h}]

A berechnet einfach die Schnittfläche der Kreise.

  ]

Um nun den Punkt zu bestimmen, an dem aus dem kleinen Kreis ein bestimmter Wert ausgeschnitten wird wollte ich diese Fläche etwa mit

FindRoot[A[r, R, x] == 0.14, {x, R}]

das geht nicht (FindRoot hat Attribut HoldAll, dann klappt das mit der Integrationsgrenze nicht). Wie kann ich es denn anstellen, eine derartige Gleichung zu lösen. Natürlich könnte man schnell ein Bisektionsverfahren oder so programmieren, ich glaube aber, daß das auch mit FindRoot geht, bekomme es aber im Moment nicht hin. Bin dankbar für einen Hinweis.

Viele Grüße

m.g.

Michael Gamer

Zur Römerbrücke 19

63456 Hanau 

Telefon: 06181 663363

Fax: 06181 663363

Mail: mg@XXXXXXX.de

home:  www.michaelgamer.de

photos:  www.mg-photo.de