VPF::[Pascal] Интеграл методом треугольников

andru2012

Дата 27.3.2014, 20:33 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 4
Регистрация: 4.3.2012

Репутация: нет
Всего: нет

У меня есть программа которая считает интеграл --- f(x)=5x^2-x+2

Код


program fourrr;

var a, b : real; { границы отрезка}
       h : real; { приращение аргумента }
       s : real; { приближенное значение интеграла }
       n : real; {количество интервалов }
       x : real; { аргумент }
       y : real; { значение функции в начале интервала }
       i : integer;
begin
    writeln('  f(x)=5x^2-x+2  ');

    write('Нижняя граница a: '); read(a);
    write('Верхняя граница b: '); read(b);
    write('приращение аргумента h: '); read(h);

     n := (b - a) / h + 1;
     x := a;
     s := 0;

    i := 1;
    while i <= n do
      begin
        y := x * x + 2;  { значение функции в начале интервала }
        s := s + (y * h);
        x := x + h;
      end;

   writeln('значение интеграла: ', s); 
end.

Мне надо переписать ее на вот такой интеграл --- f(x)=(p^5 * x * e^x) / ( (p + e^x)^6 ) * dx
где p меняется от -10 до 10 В программе использовать процедуры или функции
Огромное спасибо за помощь!

v4v4n

Дата 4.4.2014, 07:52 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 45
Регистрация: 26.9.2008

Репутация: 2
Всего: 4

andru2012, в твоем примере не метод трапеций, не функция f(x)=5x^2-x+2 и, вообще, пример не работает ))

Код


program fourrr;
var a, b : real; { границы отрезка}
       h : real; { приращение аргумента }
       s : real; { приближенное значение интеграла }
       x : real; { аргумент }
       p : real; 
  function f(p,x:real):real;
    begin
        f:=(exp(5*ln(p))*x*exp(x))/(exp(6*ln(p+exp(x))));
    end;

begin
    writeln('  f(x)=(p^5 * x * e^x) / ( (p + e^x)^6 )  ');
    write('Параметр p: '); readln(p);
    write('Нижняя граница a: '); readln(a);
    write('Верхняя граница b: '); readln(b);
    write('приращение аргумента h: '); readln(h);
    s:=(f(p,a)+f(p,b))/2;
    x:=a+h;
    while x < b do
      begin
        s := s + f(p,x);
        x := x + h;
      end;
    s:=s*h;
    writeln('значение интеграла: ', s); 
end.

Это сообщение отредактировал(а) v4v4n - 4.4.2014, 08:17

Poseidon

Дата 4.4.2014, 12:01 (ссылка)

(нет голосов)

Загрузка ...

Delphi developer

Профиль
Группа: Комодератор
Сообщений: 5273
Регистрация: 4.2.2005
Где: Гомель, Беларусь

Репутация: 5
Всего: 133

Для домашних заданий, контрольных и курсовых существует Центр помощи.

Тема перенесена!

--------------------

Если хочешь, что бы что-то работало - используй написанное,
если хочешь что-то понять - пиши сам...

1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »