Вычисление определнных интегралов, Помогите плиз, не знаю метод =(( Опции

[sintanial]

Численное интегрирование
Спасибо огромнейшее за то что выложили так много способов вычисления определенных интегралов, но всё же в универе задали вычислить интеграл методом правых прямоугольников, а здесь я нашел только центральных прямоугольников, и просто метод прямоугольников =). Я просто болел когда была эта лекция и поэтому не знаю чем отличается метод прямоугольников от метода правых прямоугольников и левых прямоугольников =), а интеграл завтра сдавать( поэтому лекцию взять не успею =(( ). Прошу помощи, за ранее огромное спасибо

[Джери]

у мня метод левых прямоугольников и вот что я пробнул набрать
-------------

Program integ;
uses crt;
var 
a,b,eps,k,sum,sum1:real;
n:integer;
function f(x:real):real;
 begin
 f:=x/cos(x);
 if x=pi/2+pi*k then
 writeln('wrong');
 end;

function calc(n:real):real;
var
h,sqr,x:real;
begin
 h:=abs(b-a)/n;
 x:=a;
 sqr:=0;
  while x<b do
   begin
     x:=x+h;
    sqr:=sqr+f(x);

  end;
   sqr:=abs(sqr*h);
   calc:=sqr;
end;

begin
TextColor(2);
 writeln('Zadanie f(x)=x/cos(x)');
 write('Vvedite levuiu granitsu a=');    {left_border}
 readln(a);
 write('Vvedite pravuiu granitsu b=');   {right_border}
 readln(b);
 write('Vvedite tochnost eps=');           {tochnost}
 readln(eps);
 sum:=0;

 if {(a=pi/2+pi*k) or} (a=b) then
    writeln('error')
 else
  begin
   n:=1;
   sum1:=calc(n);
   repeat
   sum:=sum1;
   n:=n*2;
   sum1:=calc(n);
   until abs(sum-sum1)<=eps;

   writeln(' Result= ' ,sum1);
   writeln(' n= ',n);
  end;

readln;
end.

---------
но оно как я понял:
1 работает ток в определенных значениях раз;
2 не учитываетса то что в значении х=pi/2+pi*k(real) функция y=x/cos(x) проскакивает ( например х=0 0/cos0=0)
Проблема в том что я не наю как исправить ни 1 ни 2 проблему

Сообщение отредактировано: Lapp - 16.01.2010 0:16