вычислить интеграл методом симпсона, не тот ответ Опции

[Ven]

Надо вычислить интеграл методом симпсона

пример взяла с этого форума, всё работает,вот только в ответ выводит 2, а в он-лайн подсчете интеграла на стороннем сайте ответ= -6.. где ошибка?))помогите разобраться)

program Simpson;
uses crt;
{integriruemaya funkciya}
function F(x:Real):Real;


begin

F:=exp(2*sin(x)) * sin(5 + 2* cos(x) );
end;

var  a,b,h,x: real;
     i,l,n: integer;
     integ: real;
begin
 clrscr;
 {interval integrirovaniya}

 a:=0;
 b:=2*pi;

 {osnovnaya programma}

 write ('Vvedite kolichestvo razbieniy intervala (chetnoe) n= '); readln(n);
 if (n mod 2) > 0 then
 begin
 n:=n+1;
 writeln('Chislo n bilo vvedeno nechetnoe i bilo zameneno na n=',n);
 end;
 h:=(b-a)/n;
 integ:=F(a)+F(b)+4*F(a+h);
 for i:=1 to (n div 2 ) -1 do
 begin

 x:=a+2*h*i;
 integ:=integ + 2*F(x) + 4*F(x+h);
 end;
 integ:=h*integ/3;

 l:=round(-ln(integ)/ln(10));

 writeln('Integral= ',integ:l:l);
 readln;
 end.

там в конце строчка integ:=h*integ/3;

к чему она? так было в примере решения..

Сообщение отредактировано: Ven - 13.10.2011 1:31

[Krjuger]

Ну для начала ты должна посчитать точное значение, затем ты сравниваеш точное значение с значением,полученным при n=2,если не удовлетворяет твоей точности,то увеличиваеш n на 2 и пересчитываеш симпсона и так до тех пор пока не получишь ответ,но тут есть маленькая тонкость,как и у любого другого метода,есть функции для которых метод расходится и заданную точность невозможно достигнуть,в таком случае нужно поставить какое нибудь естественно ограничение,например n=1000,и если ответ за это время не был получен,то вывести соответствующее сообщение.