Разработка программы для вычисления интеграла по формуле трапеций с тремя десятичными знаками / Turbo Pascal

[Настена95]

Здравствуйте, напишите пожалуйста программу для вычисления интеграла по формуле трапеций с тремя десятичными знаками : ∫dx/2x+3, при n= 5; a=1 b=9

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

[Настена95]

ну пожалуйста не игнорируйте

[type_Oleg]

Это - изобретение велосипеда. Поэтому никому неохота, наверное все устали.
Посмотрите например сюда - http://programmersforum.ru/showthrea...=261769&page=2
Пост #12 - только саму функцию заменить, одну строку. Ну и интервалы a,b

Насчет именно 3-х десятичных знаков - не знаю..

Впрочем, у вас функция простая, можно взять вторую производную и вычислить погрешность ( остаточный член) .Точнее, наоборот, по заданной погрешности ( 0,0001 вроде у вас )определить, какой должен быть шаг h, и сколько должно быть n=(b-a)/h

[Mad_Cat]

Насчет именно 3-х десятичных знаков - не знаю..

точность 0.001 всего лишь

n берем какое-нить по умолчанию и каждый раз увеличиваем в 2 раза
как только abs(Integral(n)-Integral(2*n))<0.0001 ==> расчет окончен

const
a = 1;
b = 9;

var n : integer;
    Int:real;
function f(x: real): real;
begin
f := 1/2*x+3;
end;

function IntegrTrap(x1,x2:real;n:Integer):real;
var i:Integer;
    dx,s:real;
begin
  s:=0;
  dx:=(x2-x1)/n;
  for i:=1 to n do
     s:=s+dx*(f(x1+i*dx)+f(x1+(i-1)*dx))/2;
  IntegrTrap:=s;
end;

begin
 n:=5;
 while abs(IntegrTrap(a,b,n)-IntegrTrap(a,b,2*n))>0.0001 do
 begin
 n:=n*2;
 end;
 Writeln('число интервалов разбиения = ',n);
 WriteLn('приближенное значение = ',IntegrTrap(a,b,n):6:3);
 ReadLn;
end.

число интервалов разбиения = 5
приближенное значение = 44.000

если же функция будет 1/(2*x+3), то ответ будет иным

число интервалов разбиения =20
приближенное значение = 0.722

upd

Три знака должны быть точными.

точными так точными, но ничего по сути не изменилось, кроме числа итераций во 2 случае

число интервалов разбиения = 80
приближенное значение = 0.722

[type_Oleg]

Mad_Cat, нет, там в прикрепленном файле в конце приведена формула остаточного члена. Не зря же ее привели.
И наверное не 0,001 , а 0,0001 или хотя бы 0,0005. Три знака должны быть точными.

[Настена95]

Огромнейшее всем спасибо, простите,что отняла ваш личное время)))