VPF::[C++] Сумма бесконечного ряда - Форум программистов

Artiom18

Дата 16.12.2008, 21:00 (ссылка)

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

Загрузка ...

Новичок

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

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

Помогите, а?

Задачка найти сумму бесконечного ряда:

y=1+1/1!+1/2!+...+1/n!
с точностью эпсилон=0,001.

Artiom18

Дата 17.12.2008, 17:57 (ссылка)

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

Загрузка ...

Новичок

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

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

Вот что получилось:

Код


#include <stdio.h>
#include <math.h>
#pragma hdrstop
#include <conio.h>
#include <condefs.h>
#pragma argsused
#include <iostream.h>
int main()
{
    float eps=0.001;
    int n=0;
    long int f=1;
    float y=0;
    float y1=0;
    float s=1;
    long int r;
m: n=n+1;
   f=f*n;
   y=1/f;
   y1=1/(f*n);
   s=s+y;
   r=fabs(y-y1);
   if ( r>=eps );
   goto m;
   cout<<"Summa="<<s<<endl;
getch();
return 0;
}

Не работает. :dash
Помогите.

Это сообщение отредактировал(а) Artiom18 - 17.12.2008, 17:59

Artiom18

Дата 20.12.2008, 22:47 (ссылка)

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

Загрузка ...

Новичок

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

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

Так тоже не работает smile

Код


#include <condefs.h>
#pragma hdrstop
#include <iostream.h>
#include <conio.h>
#include <math.h>
#pragma argsused
int main()
{
 int n=1;
 int f=1;
 float eps=0.001;
 float yp=0;
 float s=1;
 float r=0;
 float ypr=0;
m: ypr=yp+1/f;
   yp=ypr+1/(f*(n+1));
r=fabs(ypr-yp);
n=n+2;
f=f*n;
s=s+yp;
if(r>=eps) goto m;
cout<<"summa riada = " << s << endl;
getch();
        return 0;
}

IKM2007

Дата 20.12.2008, 23:03 (ссылка)

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

Загрузка ...

Зима близко

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

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

Код


#include <iostream>
int fac(int x)
{
int f = 1, i;
for( i = 2; i <= x; i++)
f *= i;
return f;
}
void main()
{
double y=0;
int n;
std::cin>>n;

for(int i = 0; i <= n; i++)
y+=(double)1/fac(i);

std::cout<<y<<std::endl;
}

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

"К чёрту обстоятельства, я создаю возможности."

Брюс Ли

volvo877

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

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

Загрузка ...

Эксперт

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

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

IKM2007, неправда, надо считать

Цитата(Artiom18 @ 16.12.2008, 20:00

)

с точностью эпсилон

, а не определенное количество элементов...

Код


#include <iostream>
#include <cmath>

using namespace std;

int main()
{
    double next = 1, s = 0;
    const double eps = 0.001;

    int i = 1;
    while(fabs(next) > eps) {
        s += next;
        next /= ++i;
    }
    cout << s << endl;
    return 0;
}

Artiom18

Дата 21.12.2008, 13:01 (ссылка)

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

Загрузка ...

Новичок

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

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

volvo877, большое спасибо.
Пришлось добавить в Borlande немного библиотек:

Код


#include <vcl\condefs.h>
#include <stdio.h>
#include <iostream.h>
#include <conio.h>
#pragma hdrstop

И в конце:

Код


getch();

И еще, ряд должен быть равен e, т.е. 2.71.
А программа не считает факториал 0!=1.
Насколько я понял. Т.е. нужно добавить еще к сумме +1.
Еще раз всем спасибо.

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