numerisches differenzieren

[Dmo11]

Hallo Leute....

Hab diesen Quellcode zum numerischen differenzieren hier im Forum gefunden.

Delphi-Quellcode:

function foo_(x: Double; degree: integer): Double;
begin
if degree = 0
  then Result := foo(x)
  else Result := (foo_(x + PRECISION, degree - 1) - foo_(x - PRECISION, degree - 1)) / (2*PRECISION);
end;

Kann mir bitte jemand erklären wie ich diese funktion für

Delphi-Quellcode:

u:=u1*sin(2*pi/m*s1*i)

anwenden kann.
Ich bräuchte davon nur lediglich die 1.Ableitung.

mfg dmo

[gammatester]

Zitat:

Zitat von Dmo11

Kann mir bitte jemand erklären wie ich diese funktion für

Delphi-Quellcode:

u:=u1*sin(2*pi/m*s1*i)

anwenden kann.
Ich bräuchte davon nur lediglich die 1.Ableitung.

Hier brauchst Du gar keine numerische Ableitung. Die Ableitung von sin(ax) ist a*cos(ax), also sollte Dein Code ungefähr so aussehen (unter der Verausstzung, daß i eine Indexvariabe ist, die dem x entspricht)

Delphi-Quellcode:

a  := 2*pi/m*s1;
u  := u1*sin(a*i);
du := u1*a*cos(a*i);

Gruß Gammatester

[Khabarakh]

Zitat:

Zitat von Dmo11

Hab diesen Quellcode zum numerischen differenzieren hier im Forum gefunden.

Wäre keine schlechte Idee, dann auch den Original-Link mit anzugeben.
"foo" ist die zu differenzierende Funktion, du musst also deine mathematische Funktion in eine neue Funktion

Delphi-Quellcode:

function foo(x : Double) : Double;

stecken. Da nicht wirklich erkenntlich ist, was nun die Funktionsvariable ist, überlasse ich das Weitere dir ;) .

[edit]
Roter Kasten :gruebel: ?
Gut, das wäre natürlich eine Alternative :mrgreen: .
[/edit]