Delphi курсовой (м. Ньютона с ф-ей Бесселя)!

[TokkyDog]

Большая просьба ко всем. Мне нужно сделать курсовой в delphi по решению СНАУ методом Ньютона и секущих, т.е. сравнить их. Загвоздочка вот в чем: нужно использовать ф-ю Бесселя нулевого порядка, а м.Ньютона нужна еще и её производная, а я не знаю где её взять. Все почти готово вот только с производными не знаю. Плиз помогите кто нить. Буду очень благодарен)))

[Lionel Messi]

Насколько я понял, услуги по созданию курсовых проектов платные, и даром никто не сделает

[Пепел Феникса]

а почему мы должны делать за вас нахаляву?вы проворонили учебу...не учились...так что тут извиняйте.
если почти все готово, то выложите ваше, и помогут(но не сделают за вас).

надо отличать помогите от сделайте за меня.

[TokkyDog]

Ну уж извините не так выразился. Мне нужна помощь.
Вот для метода Вегстейна, но здесь по ходу считает не Бесселя ф-ю?
Не могли бы вы помочь с нахождением Z для метода секущих, плиз)
Так же и м. Ньютона в отдельном юните, но там тоже есть f1 и f2, а мне для бесселя нужно

unit Vegstein;

interface
uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, ExtCtrls, StdCtrls, ComCtrls,Math,unit1,bessel;
Type Metod2=class
function f_c(x:extended):extended;
Procedure schet(a:extended;b:extended;e:exten ded);


end;
implementation

Function Metod2.f_c(x:extended):extended;
begin
Result:=x-4*cos(x);
end;



Procedure Metod2.schet(a:extended;b:extended; e:extended);
var Time1,Time2,TotalTime:TDateTime;
x,x1,z,fc,d:extended;
n:integer;
begin
x1:=(a+b)/2;
x:=x1*0.1;
n:=0;
fc:=f_c(x1);
z:=x1+((x1-x)*(fc-x1))/(1-(fc-x1));
repeat

inc(n);
x:=x1;
x1:=fc;
fc:=f_c(z);
z:=x1+(((x1-x)*(fc-x1))/(1-(fc-x1)));

until abs(fc-z)<e;

Form1.Memo2.Lines.add('Корень: '+floattostr(z));
Form1.Memo2.Lines.add('Количество итераций: '+floattostr(n));
Form1.Memo2.Lines.Add('Время выполнения: '+floattostrf(Totaltime,fffixed,8,1 0));

end;
end.

[TokkyDog]

Здесь есть и бесселя юнит, но по-моему считает ф-ю x-4*cos(x);
unit MN;
interface
uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, ExtCtrls, StdCtrls, ComCtrls,Math,unit1,bessel;
Type Metod1=class
function f1(x:extended):extended;
function f2(x:extended):extended;
procedure schet(a:extended;b:extended;eps:ext ended);

end;
implementation
Function Metod1.f1(x:extended):extended;
begin
Result:=x-4*cos(x); // Основная функция
end;

Function Metod1.f2(x:extended):extended;
begin
Result:=1+4*sin(x); //Производная функции
end;

procedure Metod1.schet(a:extended;b:extended; eps:extended);
var
n:integer;
x,t:extended;
Time1,Time2,TimeTotal:TDateTime;
begin
n:=0;
time1:=now;
x:=b;
repeat
inc(n);
t:=f1(x)/f2(x);
x:=x-t;
until abs(t)<=eps;

Form1.Memo1.Lines.Add('Корень: '+FloatToStr(x));
Form1.Memo1.Lines.Add('Количество итераций: '+FloatToStr(n));
Time2:=now;
TimeTotal:=Time2-Time1;
Form1.Memo1.Lines.Add('Время выполнения: '+floattostrf(TimeTotal,fffixed,8,1 0));
end;
end.