Ряды в С++

@marik2012 · Регистрация: 05.05.2012

Author24 — интернет-сервис помощи студентам

Найти сумму ряда с точностью E=10^-4, общий член которого An=(n!)/3(n^n)

Оригинал

Знайти суму ряду з точністю E=10^-4, загальний член якого An=(n!)/3(n^n)

Добавлено через 14 минут
вот формула $https://www.cyberforum.ru/cgi-bin/latex.cgi?n!\div 3{n}^{n}$

@zss · 05.05.2012, 19:11

Что-то я не могу найти красивого рекуррентного соотношения A(n+1)=f(An),
т.к. (n+1)^(n+1) надо раскрывать через бином Ньютона.
Предлагаю такой вариант (конечно не оптимальный):

C++

int pown(int n)
{
    int t=1;
    for(int i=1;i<=n;i++)
       t*=n;
    return t;
}
 
double sum(double eps)
{
    double an,s;
    int n;
    an=1./3.;
    s=an;
    n=2;
    while(an>eps)
    {
        an=an*n*pown(n-1)/pown(n);
        s+=an;
        n++;
    }
    return s;
}

@marik2012 · 05.05.2012, 19:22 **[ТС]**

Желательно через while бы

@zss · 05.05.2012, 19:28

А я Вам через что сделал (16 строка)?

@gGrn-7DA · 05.05.2012, 20:04

(n+1)^(n+1)/n^n в пределе имеет e*(n+1).
http://ru.wikipedia.org/wiki/Замечательные_пределы
Ну и предел A(n+1)/A(n)=1/e

Но эти формулы точны только при больших n...
Для наглядности. Погрешность при n=1000 достаточно большая
http://www.google.ru/search?q=e-(1001/1000)^1000

@marik2012 · 05.05.2012, 21:05 **[ТС]**

хотелось бы увидеть реализацию программы при маленьких n.

Jupiter · 05.05.2012, 22:24

Сообщение от marik2012

хотелось бы

Комментарий модератора

Хотелось бы от Вас увидить здравый смысл:
1. Форум с доменом .ru по-любому только русскоязычный.
2. Тема и так находится в разделе С++. Зачем ещё и к названию лепить язык!?
Следующая ваша тема с заданием на украинском языке просто будет удалена!

@gGrn-7DA · 06.05.2012, 10:43

Сообщение от marik2012

реализацию программы при маленьких n.

выше пост с кодом и есть реализация.

@marik2012 · 14.05.2012, 21:50 **[ТС]**

пишет что ошибка в 4 строке, может быть кто-то предложит другую реализацию?

@zss · 14.05.2012, 22:03

Сообщение от marik2012

пишет что ошибка в 4 строке, может быть кто-то предложит другую реализацию?

Компилируйте компилятором С++, а не С
(если используете MS VC, то расширение файла должно быть .cpp)

@marik2012 · 14.05.2012, 22:09 **[ТС]**

использую dev.cpp

@gGrn-7DA · 14.05.2012, 22:52

замени код функции.

C

int pown(int n)
{
    int t=1;
    int i;
    for(i=1;i<=n;i++)
       t*=n;
    return t;
}

еще есть вариант с использованием флага -std=c99 при компиляции

@marik2012 · 19.05.2012, 08:17 **[ТС]**

C++

#include <iostream>
#include <cmath>
using namespace std;
long factorial(int n);
int main()
{
double eps=1.0e-4; //точность расчета
    int n=1;
    double a;
    double s=0;
 
do
{
    a=factorial(n)/3*pow(n,n);
    s+=a;
    n++;
}while(a>=eps);
cout<<"s="<<s<<endl;
 
    return 0;
}
 
long factorial(int n)
{
    long f=1;
    for(int i=1; i<=n;i++)
        f*=i;
    return f;
}

подскажите правильно ли написал, если нет то исправьте ошибки.

Добавлено через 7 часов 1 минуту
но мне что-то не нравиться то что выводит на екран.

@zss · 19.05.2012, 15:26

я написал рекуррентное вычисление члена ряда
a=a*n*pown(n-1)/pown(n);
чем оно Вам не понравилось?
А Ваша формула
a=factorial(n)/3*pow(n,n);
неверна (pow(n,n) - попало в числитель).

@marik2012 · 19.05.2012, 17:06 **[ТС]**

в таком виде программа не работает.

C++

#include<iostream>
using namespace std;
 
int main()
{
int pown(int n)
 
 
{
    int t=1;
    for(int i=1;i<=n;i++)
       t*=n;
    return t;
}
 
double sum(double eps)
{
    double an,s;
    int n;
    an=1./3.;
    s=an;
    n=2;
    while(an>eps)
    {
        an=an*n*pown(n-1)/pown(n);
        s+=an;
        n++;
    }
    return s;
}
}

@zss · 19.05.2012, 17:30

Функции не вставляют внутрь других функций

C++

#include <iostream>
using namespace std;
int pown(int n)
{
    int t=1;
    for(int i=1;i<=n;i++)
       t*=n;
    return t;
}
 
double sum(double eps)
{
    double an,s;
    int n;
    an=1./3.;
    s=an;
    n=2;
    while(an>eps)
    {
        an=an*n*pown(n-1)/pown(n);
        s+=an;
        n++;
    }
    return s;
}
int main()
{
   cout<<sum(1e-6);
}

@marik2012 · 20.05.2012, 11:43 **[ТС]**

C:\Users\Admin\Desktop\1.cpp(29) : warning C4508: 'main' : function should return a value; 'void' return type assumed

1.obj - 0 error(s), 1 warning(s)
с++ 6.0 выдает такую ошибку.

@gGrn-7DA · 20.05.2012, 11:51

в конец main добавь

C++

return 0;

@marik2012 · 20.05.2012, 12:04 **[ТС]**

выводит такой результат 3,50491e-005 это точность какая?

@zss · 20.05.2012, 12:34

замените функцию pow на такую

C++

double pown(int n)
{
   double t=1;
    for(int i=1;i<=n;i++)
       t*=n;
    return t;
}

Была ошибка - переполнение разрядной сетки.

@marik2012 0 / 0 / 0 Регистрация: 05.05.2012 Сообщений: 41
		1
	Ряды в С++ 05.05.2012, 18:45. Показов 13204. Ответов 25 Метки нет (Все метки) Найти сумму ряда с точностью E=10^-4, общий член которого An=(n!)/3(n^n) Оригинал Знайти суму ряду з точністю E=10^-4, загальний член якого An=(n!)/3(n^n) Добавлено через 14 минут вот формула $https://www.cyberforum.ru/cgi-bin/latex.cgi?n!\div 3{n}^{n}$ 0

@marik2012 0 / 0 / 0 Регистрация: 05.05.2012 Сообщений: 41
	05.05.2012, 19:22 [ТС]	3
	Желательно через while бы 0

@zss Модератор 13507 / 10757 / 6412 Регистрация: 18.12.2011 Сообщений: 28,712
	05.05.2012, 19:28	4
	А я Вам через что сделал (16 строка)? 0

@gGrn-7DA мну довольно <(-__-)l 217 / 206 / 15 Регистрация: 17.01.2010 Сообщений: 2,462
	05.05.2012, 20:04	5
	(n+1)^(n+1)/n^n в пределе имеет e*(n+1). http://ru.wikipedia.org/wiki/Замечательные_пределы Ну и предел A(n+1)/A(n)=1/e Но эти формулы точны только при больших n... Для наглядности. Погрешность при n=1000 достаточно большая http://www.google.ru/search?q=e-(1001/1000)^1000 0

@marik2012 0 / 0 / 0 Регистрация: 05.05.2012 Сообщений: 41
	05.05.2012, 21:05 [ТС]	6
	хотелось бы увидеть реализацию программы при маленьких n. 0

@gGrn-7DA мну довольно <(-__-)l 217 / 206 / 15 Регистрация: 17.01.2010 Сообщений: 2,462
	06.05.2012, 10:43	8
	Сообщение от marik2012 реализацию программы при маленьких n. выше пост с кодом и есть реализация. 0

@marik2012 0 / 0 / 0 Регистрация: 05.05.2012 Сообщений: 41
	14.05.2012, 21:50 [ТС]	9
	пишет что ошибка в 4 строке, может быть кто-то предложит другую реализацию? 0

@zss Модератор 13507 / 10757 / 6412 Регистрация: 18.12.2011 Сообщений: 28,712
	14.05.2012, 22:03	10
	Сообщение от marik2012 пишет что ошибка в 4 строке, может быть кто-то предложит другую реализацию? Компилируйте компилятором С++, а не С (если используете MS VC, то расширение файла должно быть .cpp) 0

@marik2012 0 / 0 / 0 Регистрация: 05.05.2012 Сообщений: 41
	14.05.2012, 22:09 [ТС]	11
	использую dev.cpp 0

@zss Модератор 13507 / 10757 / 6412 Регистрация: 18.12.2011 Сообщений: 28,712
	19.05.2012, 15:26	14
	я написал рекуррентное вычисление члена ряда a=anpown(n-1)/pown(n); чем оно Вам не понравилось? А Ваша формула a=factorial(n)/3*pow(n,n); неверна (pow(n,n) - попало в числитель). 0

@marik2012 0 / 0 / 0 Регистрация: 05.05.2012 Сообщений: 41
	20.05.2012, 11:43 [ТС]	17
	C:\Users\Admin\Desktop\1.cpp(29) : warning C4508: 'main' : function should return a value; 'void' return type assumed 1.obj - 0 error(s), 1 warning(s) с++ 6.0 выдает такую ошибку. 0

@marik2012 0 / 0 / 0 Регистрация: 05.05.2012 Сообщений: 41
	20.05.2012, 12:04 [ТС]	19
	выводит такой результат 3,50491e-005 это точность какая? 0