Turbo Pascal, метод Симпсона

[venera-8]

Turbo Pascal, метод Симпсона

Очень нужна помощь.. В общем суть- вычислить интеграл методом симпсона.. Я написала,всё считает,но, оказываетсЯ,еще нужен эпсилон- точность,которую надо забивать при расчете( напр. 0,1 ;0,01 и т.д.). Я не могу сообразить как вставить эпсилон сюда.. Я учусь только программировать, подошла к преподу- а он даже объяснить ничего не может((надеюсь на вашу помощь.... уже долго сижу с этой проблемой и никак не могу понять... вот мой код:

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,n: integer;  {n-chislo tochek na interv}
     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;



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

 

[Vladimir_S]

Посмотрите эту тему:
http://www.tehnari.ru/f43/t48865/
Там программа по расчету численного интеграла тремя методами с заданной точностью. Суть: вот Вы задали число разбиений промежутка, сосчитали интеграл, запомнили результат. Потом удвоили число разбиений, снова сосчитали и сравнили результаты. Если расхождение не превышает заданной точности, то всё, задача решена, а если превышает - то снова мельчите шаг вдвое и т.д.

И еще - во избежание неприятностей, НЕ ДУБЛИРУЙТЕ ТОПИКИ! У нас это не допускается. Дубликат удален.

 

[venera-8]

объясню что хочет препод..(я не знаю как реализовать,уже как только не делала). В общем вот у меня например при разбиении 100 и эпсилоне 0,1 выдает ответ -6,0325642 (это пример,ответ другой)
при 0,01 должен быть такой чтоб совпадала первая цифра после запятой, при 0,001 - две цифры после запятой и т.д. -6,0325642 потом 6,0345664 потом 6,0326456 и т.д. Я уже полгода ничего не могу сделать..и он объяснить не может... я делала с эпсилоном, но ему не понравилось, циферки там не менялись.. только в конце при больших эпсилон.. На самом деле нужна помощь...

 

[Vladimir_S]

Девушка, мне кажется, что предложенный мною метод как раз и решает именно эту задачу. Обратите внимание - там каждая из процедур заканчивается строкой
Until ABS(Sum2-Sum1)<0.001;
Это и означает, что при измельчении шага вдвое изменяются только значащие цифры, меньшие порядка, указанного в этом самом эпсилон (там фиксированное значение 0.001, но можно, естественно, подставить любое: 0.1, 0.01, 0.0001 и т.п.). Чего же Вам еще?