VPF::Вычисление опред. интеграла разными методами

LastOne

Дата 8.5.2010, 18:25 (ссылка)

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

Загрузка ...

Новичок

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

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

Очень прошу помочь....задали курсовик: "Вычислить определенный интеграл методами правых и левых прямоугольников,трапеций и методом Симпсона",времени мало осталось...а я еще новичек в С++. :(

вот интеграл : sin(ln(x))dx/x в границах от 1 до e и с точностью eps=0,0001

Заранее СПАСИБО!!!

toxx

Дата 8.5.2010, 19:21 (ссылка)

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

Загрузка ...

Опытный

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

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

тут посмотри.
если писать сам не хочешь тут лучше создавать тему тут

LastOne

Дата 8.5.2010, 19:35 (ссылка)

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

Загрузка ...

Новичок

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

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

ну я написал большую часть,но мне кажется,что считает не правильно...и не хватает одного метода.
Вот код,если интересно

Код


#include <stdio.h>
#include <stdlib.h>
#include <conio.h>
#include <math.h>

int main()
{
const int n=100;
double a=1, b=exp(x);
double h,x;
double F[n],I1,I2,tmp;
int i;
h=(b-a)/n;

for(i=0,x=a; i<n; i++,x+=h)
 F[i] = sin(log(x))/x;

I1 = 0;
for(i=1;i<n-1;i++)
 I1 += F[i];
I1 += F[0]/2 + F[n-1]/2;
I1 *= h;
//-------------
I2 = 0;
for(i=1; i<n; i+=2)
 I2 += F[i];
I2 *= 4;
tmp = 0;
for(i=2; i<n-1; i+=2)
 tmp += F[i];
I2 += tmp*2 + F[0] + F[n-1];
I2 *= h/3;

printf("\n integral po formule trapeciy raven = %f \n",I1);
printf("\n integral po formule Simpsona raven %f \n",I2);

getch();
return 0;
}

toxx

Дата 8.5.2010, 19:42 (ссылка)

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

Загрузка ...

Опытный

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

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

Цитата

но мне кажется,что считает не правильно...

почему так кажется?

а вы считали вручную?

LastOne

Дата 8.5.2010, 19:52 (ссылка)

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

Загрузка ...

Новичок

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

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

Нет.Просто ответы сильно отличаются друг от друга.

toxx

Дата 8.5.2010, 20:00 (ссылка)

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

Загрузка ...

Опытный

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

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

определение переменной выше чем вычисление

Код


double a=1, b=exp(x);
double h,x;

у вас х не инициализировано тут

Код


double a=1, b=exp(x);

должно быть так

Код


double h,x=1.0;
double a=1.0, b=exp(x);

Алгоритмы у вас вроде правильные, поэтому ошибки наверно тока в этих строчках

Это сообщение отредактировал(а) toxx - 8.5.2010, 20:09

LastOne

Дата 10.5.2010, 10:42 (ссылка)

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

Загрузка ...

Новичок

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

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

Да,точно!!! smile

Спасибо за помощь!

Только вот с методом прямоугольников вобще разобраться не могу,делал как по ссылке на алгоритмы,но ничего не получилось....
Не подскажете как его написать?

Всё,разобрался!Уже не надо smile

Это сообщение отредактировал(а) LastOne - 10.5.2010, 21:29

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