Pow - некорректность возведения в степень

@RINGARD · Регистрация: 27.10.2012

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

Написал программу решения СМО, однако при выполнении натрнулся на некорректность возведения в степень, что с помощю pow, что самописного возведения, но выяснил если pow проводить не над вычисляемыми значениями, то результат корректен, однако вывод самих вычесляемых переменных показывает, что они вполне нормальные.

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
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
  //---------------------------------------------------------------------------
 
#include <iostream>
 using namespace std;
#include <vcl.h>
#include <math.h>
#include <stdlib.h>
#include <stdio.h>
#include <iomanip>
#include <vcl.h>
#include <string.h>
#include <math.h>
#include <sysmac.h>
#pragma hdrstop
 
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
double m, n, T, aL, A, ju, k, plas=0, g=0;
double d=1;
//---------------------------------------------------------------------------
double fact(double a)
{
        double r=1;
        for(double i = 1; i <= a; ++i) {r = r * i;}
        return r;
}
//---------------------------------------------------------------------------
float step(double x, double y)
{
        double r;
        r = exp(y * log(x));
        return r;
}
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::Button1Click(TObject *Sender)
{
        m=StrToFloat(Edit1->Text);
        n=StrToFloat(Edit2->Text);
        T=StrToFloat(Edit3->Text);
        aL=StrToFloat(Edit4->Text);
        ju=d/T;
        A=aL/ju;
        for(k=0; k<m; k++){
                if(k<n)
                        g=(((fact(m))*(pow(A,k))) / ((fact(k))* (fact(m-k))));
                else
                        g=(((fact(m))*(pow(A,k))) / ((fact(n))*(pow(n,k-n))*(fact(m-k))));
                plas=plas+g;
 
        }
        /*
        double g=0.4;
        double e=10; */
        Memo1->Lines->Add(d/plas);
       /* Memo1->Lines->Add(k);
        Memo1->Lines->Add(A);
        Memo1->Lines->Add(pow(A,k-1));
        Memo1->Lines->Add(pow(g,e));  */
        double pp=plas;
        plas=d/plas;
        for(k=1; k<=m; k++){
                if(k<n){
                        g=((((fact(m))*(pow(A,k))) / ((fact(k))* (fact(m-k))))*pp);
                        Memo1->Lines->Add(g);
                        plas=plas+g;
                        }
                if(k>=n)
                        {
                        g=((((fact(m))*(pow(A,k))) / ((fact(n))*(pow(n,k-n))*(fact(m-k))))*pp);
                        Memo1->Lines->Add(g);
                        plas=plas+g;
                        }
        }
        Memo1->Lines->Add(plas);
 
}
//---------------------------------------------------------------------------

AntonChik · 25.04.2013, 06:34

что в вашем понимании вычисляемые и не вычисляемые значения, нормальные и ненормальные?
такой код выдает корректные значения?

C++
1
2
3
4
5
for(k=1; k<=10; k++)
{
g=pow(A,k);
Memo1->Lines->Add(g);
}

@BRcr · 25.04.2013, 07:51

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
        for(k=1; k<=m; k++){ 
                if(k<n){
                        g=((((fact(m))*(pow(A,k))) / ((fact(k))* (fact(m-k))))*pp);
                        Memo1->Lines->Add(g);
                        plas=plas+g;
                        }
                if(k>=n)
                        {
                        g=((((fact(m))*(pow(A,k))) / ((fact(n))*(pow(n,k-n))*(fact(m-k))))*pp);
                        Memo1->Lines->Add(g);
                        plas=plas+g;
                        }
                if ( ::errno == ERANGE ) Memo1->Lines->Add( "ошибка диапазона" ); 
                if ( ::errno == EDOM ) Memo1->Lines->Add( "некорректные входные значения" );
        }

Еще неплохо бы глянуть вывод такого кода.

Новые блоги и статьи Все статьи Все блоги /
делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .
Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .