Sobre:
Cita:
|
Empezado por roman
¡Eeeeh! ¡No vale! Yo quería ver la implementación de esta función!

|
Puedes encontrar la implementación si búscas información sobre cálculo numérico, en general pueden haber varias implementaciones según el método que decidas usar, aqui te va una de las implementaciones:
Código:
function dfx(fx:tfx;x:real;h:real):real; {funcion que deriva fx respecto a x}
var ii,ji,ki:integer;
p,s:real;
xi:array[0..3]of real;
const
n=3; {tercer grado}
begin
Result:=0;
xi[0]:=x-2*h; xi[1]:=x-h; xi[2]:=x; xi[3]:=x+h;
for ii:=0 to n do
begin
s:=0;
for ki:=0 to n do if ki<>ii then
begin
p:=1;
for ji:=0 to n do if (ji<>ki) and (ji<>ii) then p:=p*(x-xi[ji]);
s:=s+p;
end;
p:=1;
for ji:=0 to n do if ji<>ii then p:=p*(xi[ii]-xi[ji]);
Result:=Result+s*fx(xi[ii])/p
end;
end;
Servido
