Форум программистов [Powered by Invision Power Board]

Версия для печати темы
Нажмите сюда для просмотра этой темы в оригинальном формате

Форум программистов > Центр помощи > [C++] Вычислить сумму членов бесконечного ряда с заданно

Автор: Daremez 15.1.2012, 18:51

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

Автор: Fynivx 16.1.2012, 05:00

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

Код


#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

Цитата(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 и т.д.

Автор: bsa 16.1.2012, 14:21

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

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

Автор: Daremez 16.1.2012, 15:50

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

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

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

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

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

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

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

Автор: feodorv 16.1.2012, 19:18

Цитата(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;
}

???