Разработка алгоритма вычисления параметрического интеграла методом Симпсона

[killbore]

Кто разбирается, помогите, плиз. Очень нужно. Только начинаю изучать программирование и с математикой такой не совсем лады, а тут бац:

[p51x]

Метод Симпсона описан много где, даже с исходниками. Математики тут совсем нет, кроме значков...

[Stilet]

http://www.programmersforum.ru/showthread.php?t=56969
Так пойдет?

И
http://www.delphiworld.narod.ru/_all_articles_.html
Раздел Математика - [80]

Ну и от меня на всякий случай:

MODULE Симпсон;
IMPORT StdLog,Math,In;
VAR
 par:REAL;

PROCEDURE Fn(x:REAL):REAL;
BEGIN
 RETURN Math.Tan(x);
END Fn;

PROCEDURE Интеграл(a,b:REAL;n:INTEGER):REAL;
VAR x,dx,Res:REAL;i:INTEGER;
BEGIN
 dx:=(b-a)/n; Res:=0;
 x:=a;
 In.Open;
 (*FOR i:=1 TO n DO*)
 REPEAT 
  In.Real(x);
  StdLog.Real(x);
  Res:=Res+dx*(Fn(x)+4*Fn(x+dx/2)+Fn(x+dx))/6;
  (*x:=x+dx; *)
 UNTIL ~In.Done; 
 RETURN Res;
END Интеграл;

PROCEDURE Интергировать*;
VAR n:INTEGER;
  a,b,Eps,Prev,Res:REAL;
BEGIN
 a:=7;b:=9;Eps:=0.001;
 StdLog.Clear;
  Res:=Интеграл(a,b,4);
  n:=4;
  REPEAT
   Prev:=Res;
   n:=n+n;
  UNTIL (Res=0)OR(ABS((Res-Prev)/Res)<Eps); 
  StdLog.Ln; 
  StdLog.String('Результат');
  StdLog.Tab;
  StdLog.RealForm(Res,5,0,0,' ');
END Интергировать;


END Симпсон.

Результат:

10.0 20.0 30.0 40.0 50.0 50.0
Результат -0.13279

Чесно говоря не проверял на вшивость.

[killbore]

Thanx, Stilet.
P.S. Блин, на VB6 было б супер..