VPF::Получить число e - Форум программистов

Kipter

Дата 3.10.2007, 14:39 (ссылка)

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

Загрузка ...

Шустрый

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

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

Народ у нас есть число e = 2.7182818284590

Я пытаюсь написать программу чтобы его получить.....
по формуле 1/fact(0) + 1/fact(1) ... + 1/fact(n)

Код


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

long double fact(int n)
{
    int i;
    int p=1;

    for(i=1;i<=n;i++)
        p*=i;

    return p;
}

int main()
{
    long double e=0;
    long int i = 0;

    while(i < 15)
    {
        e += (1.0 / fact(i));

        printf("%0.13f \n", e);

        i++;
    }

    printf("\n%0.13f \n", e);
}

Собственно проблемы с точностью, немогу получить именно 2.7182818284590

Это сообщение отредактировал(а) Kipter - 3.10.2007, 14:40

archimed7592

Дата 3.10.2007, 14:44 (ссылка)

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

Загрузка ...

Архимед

Профиль
Группа: Завсегдатай
Сообщений: 2531
Регистрация: 12.6.2004
Где: Moscow

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

Цитата(Kipter @ 3.10.2007, 14:39

)

Собственно проблемы с точностью

Угу...

Код


long double fact(int n)
{
    long double p=1;
    for(int i = 1; i <= n; ++i)
        p*=static_cast< double >(i);

    return p;
}

Добавлено через 43 секунды
Ибо fact(15) ~ 1'307'674'368'000.

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

If you have an apple and I have an apple and we exchange apples then you and I will still each have one apple. But if you have an idea and I have an idea and we exchange these ideas, then each of us will have two ideas.
© George Bernard Shaw

UnrealMan

Дата 3.10.2007, 16:18 (ссылка)

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

Загрузка ...

Опытный

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

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

Цитата(Kipter @ 3.10.2007, 15:39

)

по формуле 1/fact(0) + 1/fact(1) ... + 1/fact(n)

З-з-з-з-з-з-з-зачем?
exp(1) - вот те и число e.

archimed7592

Дата 3.10.2007, 17:04 (ссылка)

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

Загрузка ...

Архимед

Профиль
Группа: Завсегдатай
Сообщений: 2531
Регистрация: 12.6.2004
Где: Moscow

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

Цитата(UnrealMan @ 3.10.2007, 16:18

)

exp(1) - вот те и число e.

Так не интересно smile

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

MAKCim

Дата 4.10.2007, 08:45 (ссылка)

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

Загрузка ...

Воін дZэна

Профиль
Группа: Экс. модератор
Сообщений: 5644
Регистрация: 10.12.2005
Где: Менск, РБ

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

Kipter,
так никто не считает
нужно использовать реккурентную формулу

Код


static inline void calculate(double * result, int count) {
    asm {
        fninit
        fld1
        fld1
        fld1
        mov ecx, count
        fdiv
        fadd st2, st0
        fld
        faddp st2, st0
        dec ecx
        db 75h, F5h
        fxch st2, st0
        fst qword ptr [result]
    }
}

Это сообщение отредактировал(а) MAKCim - 4.10.2007, 09:05

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

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