Ошибка в программе

@Vaden · Регистрация: 25.04.2019

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

Написал программу, не могу понять в чем ошибка, нужна помощь

Задача:

Для х изменяющегося от a до b с шагом (b-a)/k, где (k=10), вычислить функцию f(x), используя ее разложение в степенной ряд в трех случаях:

а) для заданного n;

б) для заданной точности e (e=0.0001);

с) для «точного» значения (по аналитической формуле).

Код:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
#include <iostream>
#include <iomanip>
#include <cmath>
#include <cstdlib>
#include <stdio.h>
using namespace std;
double Se(double x, double e);
 
double f(double x){                     
    double pi = 3.1415926;
    double y  = ((((1.0+x*x)/2)*(atan (x) * 180.0 / pi))-(x/2));
    return y;
}
 
double S(double x, long n){  // ПЕРВЫЙ ЭЛЕМЕНТ СУММЫ РЯДА
    double pi = 3.1415926;
    double ak = (pow(x,3.0)/3.0);
 
    double sum = 0;
    for(long k = 0; k < n; k++)
    {
        sum += ak;
        ak  *= k*x*(((2.0*((-1.0)^(n+2.0))) *n*x^2.0 -((-1.0)^n+2.0) *x^2.0)/(2.0*n*((-1.0)^n+1.0) +3.0*((-1.0)^n+1.0))); //РЕКУРЕНТНАЯ Ф.
 
            }
    return sum;
}
 
 
 
int main(){
 
    double x, xn = 0.1, xk = 1.0, dx, e1, e2, y1, y2, y3;
    dx = ((xk-xn)/10.0);
    printf("+--------+--------------+---------------+------------+---------------+-----------------+\n");
    printf("|    x   |      y1      |      y2       |      y3    |       e1      |       e2        |\n");
    printf("+--------+--------------+---------------+------------+---------------+-----------------+\n");
    for(x = xn; x < xk; x += dx) {
 
 
        y1=S(x, 30.0);
        y2=Se(x, 0.0000001);
        y3=f(x);
        e1=fabs( (y3-y1)/y3  );
        e2=fabs( (y2-y1)/y2  );
        printf("| %6.2f |%12.5g  |%12.5g   |  %6.3f    |         %6.3f|      %6.3f     |\n",
               x,y1,y2,y3,e1,e2);}
printf("+--------+--------------+---------------+------------+---------------+-----------------+\n");
     system("pause");
    return 0;
}
double Se(double x, double e){   //СУММА РЯДА С ТОЧНОСТЬЮ ЭПСИЛОН
    double sum = 0;
    double pi = 3.1415926;
    double an  = (pow(x,3.0)/3.0);
    for(double n = 1.0; e < fabs(an); n = n + 1.0)
    {
        sum += an;
        an  *=(x^3.0*x)/3.0);
    }
    return sum;
}

Написано на языке С++ в codeblocks

@Yetty · 15.05.2019, 21:45

Vaden, сбросьте скрин функции и её разложения в ряд

@oleg-m1973 · 15.05.2019, 21:55

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

Написал программу, не могу понять в чем ошибка, нужна помощь

А оператор ^ это разве не bitwise XOR? Там, наверное, надо pow()

@Vaden · 15.05.2019, 21:57 **[ТС]**

1-ая строка - функция, 2- диапазон изменения аргумента, 3- значение n, 4- сумма

@Vaden · 15.05.2019, 21:59 **[ТС]**

код исправил

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
#include <iostream>
#include <iomanip>
#include <cmath>
#include <cstdlib>
#include <stdio.h>
using namespace std;
double Se(double x, double e);
 
double f(double x){                    
    double pi = 3.1415926;
    double y  = ((((1.0+x*x)/2)*(atan (x) * 180.0 / pi))-(x/2));
    return y;
}
 
double S(double x, long n){  // ПЕРВЫЙ ЭЛЕМЕНТ СУММЫ РЯДА
    double pi = 3.1415926;
    double ak = (pow(x,3.0)/3.0);
 
    double sum = 0;
    for(long k = 0; k < n; k++)
    {
        sum += ak;
        ak  *= k*x*(((2.0*(pow(-1.0,n+2.0))) *n*pow(x,2.0) -(pow(-1.0,(n+2.0))) *pow(x,2.0))/(2.0*n*(pow(-1.0,(n+1.0))) +3.0*(pow(-1.0,(n+1.0)))));//рекурентная ф.
            }
    return sum;
}
 
 
 
int main(){
 
    double x, xn = 0.1, xk = 1.0, dx, e1, e2, y1, y2, y3;
    dx = ((xk-xn)/10.0);
    printf("+--------+--------------+---------------+------------+---------------+-----------------+\n");
    printf("|    x   |      y1      |      y2       |      y3    |       e1      |       e2        |\n");
    printf("+--------+--------------+---------------+------------+---------------+-----------------+\n");
    for(x = xn; x < xk; x += dx) {
 
 
        y1=S(x, 30.0);
        y2=Se(x, 0.0000001);
        y3=f(x);
        e1=fabs( (y3-y1)/y3  );
        e2=fabs( (y2-y1)/y2  );
        printf("| %6.2f |%12.5g  |%12.5g   |  %6.3f    |         %6.3f|      %6.3f     |\n",
               x,y1,y2,y3,e1,e2);}
printf("+--------+--------------+---------------+------------+---------------+-----------------+\n");
     system("pause");
    return 0;
}
double Se(double x, double e){   //СУММА РЯДА С ТОЧНОСТЬЮ ЭПСИЛОН
    double sum = 0;
    double pi = 3.1415926;
    double an  = (pow(x,3.0)/3.0);
    for(double n = 1.0; e < fabs(an); n = n + 1.0)
    {
        sum += an;
        an  *=(pow(x,(3.0*x))/3.0);
    }
    return sum;
}

Только теперь проблема другая - результаты слишком сильно разнятся :

+--------+--------------+---------------+------------+---------------+-----------------+
| x | y1 | y2 | y3 | e1 | e2 |
+--------+--------------+---------------+------------+---------------+-----------------+
| 0.10 | 0.00033333 | 0.00040014 | 2.834 | 1.000| 0.167 |
| 0.19 | 0.0022863 | 0.0026259 | 5.478 | 1.000| 0.129 |
| 0.28 | 0.0073173 | 0.0082627 | 8.294 | 0.999| 0.114 |
| 0.37 | 0.016884 | 0.018983 | 11.357 | 0.999| 0.111 |
| 0.46 | 0.032445 | 0.036626 | 14.735 | 0.998| 0.114 |
| 0.55 | 0.055458 | 0.06333 | 18.488 | 0.997| 0.124 |
| 0.64 | 0.087381 | 0.10178 | 22.670 | 0.996| 0.141 |
| 0.73 | 0.12967 | 0.15573 | 27.326 | 0.995| 0.167 |
| 0.82 | 0.18379 | 0.23106 | 32.496 | 0.994| 0.205 |
| 0.91 | 0.25119 | 0.33838 | 38.211 | 0.993| 0.258 |
| 1.00 | 0.33333 | 0.5 | 44.500 | 0.993| 0.333 |
+--------+--------------+---------------+------------+---------------+-----------------+

@Vaden · 15.05.2019, 22:01 **[ТС]**

...

Verevkin · 15.05.2019, 22:02

^ - это не возведение в степень. Тут тебе не бейсик.

@Vaden · 15.05.2019, 22:04 **[ТС]**

Уже исправил, посмотрите выше.

@oleg-m1973 · 15.05.2019, 22:12

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

an **=(pow(x,(3.0*x))/3.0);

Судя по предыдущему коду здесь должно быть (pow(x,3.0) * x)/3.0); (о формуле молчу)

Добавлено через 3 минуты

Сообщение от oleg-m1973

Судя по предыдущему коду здесь должно быть (pow(x,3.0) * x)/3.0); (о формуле молчу)

Хотя, скорее ошибаюсь

@Vaden · 15.05.2019, 22:27 **[ТС]**

После исправления (pow(x,(3.0*x))/3.0); на (pow(x,3.0) * x)/3.0); результаты те же

@Yetty · 16.05.2019, 08:35

Vaden,

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
double Sum(double x, int n)
{
    int p;
    double a=-x, S=0.0;
    for(int i = 1; i <= n; i++)
    {       
        p=4*i*i-1;
        a*=-x*x;
        S+=a/p;
    }
    return S;
}
 
double Sum_e(double x, double e)
{   
    int i=2, p=3; 
    double a=x*x*x, S=0.0;
    while(fabs(a/p)>e)
    {        
        S+=a/p;        
        p=4*i*i-1;
        a*=-x*x;
        i++;
    }
    return S;
} 
 
double y(double x)
{   
    return (1.0+x*x)/2.0*atan(x)-x/2.0;    
} 
 
int main()
{   
    int n=30, k=10;
    double a=0.1, b=1.0, e=0.0001, h=(b-a)/k;
    
    cout << " x        S          Se         y\n";
    for (double x=a; x<=b+e/2.0; x+=h)    
    cout << setprecision(2) << fixed << x << "  " 
    << setprecision(7) << Sum(x,n) << "  " << Sum_e(x,e) << "  " << y(x) << "\n";    
    
system("pause");
return 0;
}

Новые блоги и статьи Все статьи Все блоги /
Новый CodeBlocs. Версия 25.03 palva 04.01.2026 Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .	Модель микоризы: классовый агентный подход anaschu 02.01.2026 Раньше это было два гриба и бактерия. Теперь три гриба, растение. И на уровне агентов добавится между грибами или бактериями взаимодействий. До того я пробовал подход через многомерные массивы,. . .	Учёным и волонтёрам проекта «Einstein@home» удалось обнаружить четыре гамма-лучевых пульсара в джете Млечного Пути Programma_Boinc 01.01.2026 Учёным и волонтёрам проекта «Einstein@home» удалось обнаружить четыре гамма-лучевых пульсара в джете Млечного Пути Сочетание глобально распределённой вычислительной мощности и инновационных. . .	Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .	Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .
Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Рецензия / Мнение/ Перевод Нашел на реддите интересную статью под названием The Thinkpad X220 Tablet is the best budget school laptop period . Ниже её машинный перевод. Thinkpad X220 Tablet —. . .	PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .

@Yetty 7438 / 5030 / 2892 Регистрация: 18.12.2017 Сообщений: 15,692
	15.05.2019, 21:45
	Vaden, сбросьте скрин функции и её разложения в ряд 0

@Vaden 0 / 0 / 0 Регистрация: 25.04.2019 Сообщений: 63
	15.05.2019, 21:57 [ТС]
	1-ая строка - функция, 2- диапазон изменения аргумента, 3- значение n, 4- сумма Миниатюры 0

@Vaden 0 / 0 / 0 Регистрация: 25.04.2019 Сообщений: 63
	15.05.2019, 22:01 [ТС]
	... Миниатюры 0

Verevkin Злостный нарушитель 10246 / 5673 / 1262 Регистрация: 12.03.2015 Сообщений: 26,284
	15.05.2019, 22:02
	^ - это не возведение в степень. Тут тебе не бейсик. 0

@Vaden 0 / 0 / 0 Регистрация: 25.04.2019 Сообщений: 63
	15.05.2019, 22:04 [ТС]
	Уже исправил, посмотрите выше. 0

@Vaden 0 / 0 / 0 Регистрация: 25.04.2019 Сообщений: 63
	15.05.2019, 22:27 [ТС]
	После исправления (pow(x,(3.0x))/3.0); на (pow(x,3.0) x)/3.0); результаты те же 0