VPF::[C++] Вычислить сумму членов бесконечного ряда с заданно

Daremez

Дата 15.1.2012, 18:51 (ссылка)

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

Загрузка ...

Новичок

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

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

Помогите решить - Си практически незнаю. Может кто подкинет код?

Присоединённый файл ( Кол-во скачиваний: 30 )

__________.PNG 11,31 Kb

Fynivx

Дата 16.1.2012, 05:00 (ссылка)

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

Загрузка ...

Шустрый

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

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

Не знаю, правильно ли, но вот самый простой, первый, пришедший в голову, вариант:

Код


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

unsigned long factorial ( int x ) {
    unsigned long result = 1;
    for ( int i = 1; i <= x; i++) result *= i;
    return result;
}

float calc ( float e, float x ) {
    float fn = e;
    float result = 0;
    for ( int n = 0; fn >= e; n++) {
        fn = pow( x, 2 * n + 1 ) / factorial( 2 * n + 1 );
        result+= fn * pow( 0-1, n );
    }
    return result;
}

int main ( int argc, char** argv ) {
    float e;
    float x;
    printf( "Введите e:\n" );
    scanf( "%f", e );
    printf( "Введите x:\n" );
    scanf( "%f", x );
    printf( "Ответ: %f\n", calc( e, x ) );
}

А вообще здесь для таких вещей специальный раздел есть.
Только вот проблема - факториал семнадцати (девятая итерация) уже не влезет в 64 бита, а бигинт использовать - уж точно не этого уровня задачка.

feodorv

Дата 16.1.2012, 07:00 (ссылка)

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

Загрузка ...

Эксперт

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

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

Цитата(Fynivx @ 16.1.2012, 05:00

)

pow( x, 2 * n + 1 )

Возведение в степень требуется выполнять путём умножения ;) Да и факториал не стОит так уж в лоб вычислять...

В задаче легко можно вычислить следующий член ряда, зная предыдущий:

Код


double cur;
double x2;
unsigned int num;

void getNext( void )
{
  num = num + 1;
  cur = -cur;
  cur = cur * x2;
  cur = cur / ((double) (2*num) * (double) (2*num+1));
}

Ну и инициализация и проверка на эпсилон:

Код


double getSum( double x, double e)
{
  double sum = x;

  num = 0;
  cur = x;
  x2 = x * x;
  
  while( fabs(cur) >= e )
  {
     getNext();
     sum = sum + cur;
  }

  return sum;
}

Итого:

Код


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

double cur;
double x2;
unsigned int num;

void getNext( void )
{
  num = num + 1;
  cur = -cur;
  cur = cur * x2;
  cur = cur / ((double) (2*num) * (double) (2*num+1));
}

double getSum( double x, double e)
{
  double sum = x;

  if( e <= 0 ) e = 1.e-6;
  num = 0;
  cur = x;
  x2 = x * x;

  while( fabs(cur) >= e )
  {
     getNext();
     sum = sum + cur;
  }

  return sum;
}

int main( void )
{
  float x = 0.5, e = 1.e-6;

  printf( "Введите число:\n" );
  scanf( "%f", &x);

  printf( "Введите точность:\n" );
  scanf( "%f", &e);

  printf( "sum = %e, sin = %e\n", getSum( x, e), sin(x));
}

Цитата(Daremez @ 15.1.2012, 18:51

)

Си практически незнаю

Гм. Стоит в него углубиться, иначе так и не понятно будет, что это за printf/scanf, 1.e-6, %f/%e, почему везде double, а "x" и "e" почему-то float и т.д.

Это сообщение отредактировал(а) feodorv - 16.1.2012, 07:04

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

Напильник, велосипед, грабли и костыли - основные инструменты программиста...

bsa

Дата 16.1.2012, 14:21 (ссылка)

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

Загрузка ...

Эксперт

Профиль
Группа: Модератор
Сообщений: 9185
Регистрация: 6.4.2006
Где: Москва, Россия

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

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

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

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

Правильно заданный вопрос - половина ответа
Отзывы о работе модератора.

Daremez

Дата 16.1.2012, 15:50 (ссылка)

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

Загрузка ...

Новичок

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

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

вообще то так и планировал углубиться но сессия пришла раньше.

feodorv, твой код "итого" у меня запускается но пишет ошибку машинный код завершился 0х0000013 (0х0)
Что делать?

Добавлено через 54 секунды

Цитата(bsa @ 16.1.2012, 14:21)

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

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

bsa хорошо буду знать.

feodorv

Дата 16.1.2012, 19:18 (ссылка)

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

Загрузка ...

Эксперт

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

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

Цитата(Daremez @ 16.1.2012, 15:50

)

твой код "итого" у меня запускается но пишет ошибку машинный код завершился 0х0000013 (0х0)

ФигОво

А чем и на чём компилишь?

Добавлено @ 19:20
А если main() переписать так:

Код


int main( void )
{
  double x = 0.5, e = 1.e-6;
  printf( "sum = %e, sin = %e\n", getSum( x, e), sin(x));
  return 0;
}

???

Это сообщение отредактировал(а) feodorv - 16.1.2012, 19:23

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

Напильник, велосипед, грабли и костыли - основные инструменты программиста...

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