Поиск:

Ответ в темуСоздание новой темы Создание опроса
> не могу запустить программу! 
:(
    Опции темы
Dmi3ev
Дата 10.3.2009, 18:27 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


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

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



Цитата

Dmi3ev, а чем ты компилируешь, а то у меня возникает ошибка(в C++ Builder 2007): Floating point overflow...

ms vs 2008
в какой строке ошибка?
Цитата

Floating point overflow...

сделай так, чтобы нули не получались и все будет ок... скорее всего поэтому ошибка...

Это сообщение отредактировал(а) Dmi3ev - 10.3.2009, 18:40


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

PM MAIL   Вверх
molodzo
Дата 10.3.2009, 19:30 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

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



Пробую откомпилировать в Visual Studio 2005. Создал progect(по умолчанию), код вставил в name.h... А вот что нужно выбрать при: executable for debug session?
PM MAIL   Вверх
molodzo
Дата 10.3.2009, 22:48 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

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



Ура! Удалось запустить! Только вот странность: при любом значении step, программа выдает одно и тоже значение (k = 17)... Кто-нибудь знает в чем ошибка?
PM MAIL   Вверх
molodzo
Дата 10.3.2009, 23:38 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

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



Исправил все значения на double... прога стала выдавать значения k = 43 при любом x!
Если вместо float поставить long double, то k уже будет равняться 50!
Также не получается нормально вывести left(пишет left = 1.#IND), а вот у right нормальное число!
Кто-нибудь знает почему такое происходит и как это исправить?

Код

#include "stdafx.h"
#include <math.h>
#include <stdio.h>
#include <iostream>
#include <conio.h>
using namespace std;

double fact(double val)
{
    if (val==0)
        return 1;
    else
        return (val*fact(val-1));
}
int _tmain()
{
   double x;    
   double k = 0;
   double left = 0;
   double right = 0;
   double step = pow(0.1, 15);
    do
    {
    cout << "Vvedite x: ";
    cin >> x;
    }
    while ((x>=0.25) || (x<=-0.25));
    do
    { 
        k++;
        left = (fact(4 * k + 1) / pow((long double)fact(2 * k), 2)) * pow(x, 2 * k);    
        right = 0.5 *  (pow((1/(1 - 4 * x)), 1.5)) + pow((1/(1 + 4 * x)), 1.5) ;        
    }
    while (fabs(right - left) >= step);
    cout << "k = " << k << endl;
    cout << "right = " << right << endl;
    cout << "left = " << left << endl;
    system("pause");
    return 0;
}

PM MAIL   Вверх
Dmi3ev
Дата 10.3.2009, 23:55 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


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

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



Цитата

Кто-нибудь знает в чем ошибка? 

короче, molodzo, программа у тебя правильно работает, осталось только предусмотреть те случае, когда могут быть 0 (в отрицалтельной степени) и еще немного въехать в то, что компьютер некоторые числа воспринимает одинаково, хотя они не равны...

Добавлено через 12 минут и 23 секунды
molodzo, я понял, где косяк, завтра выложу полное решение...


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

PM MAIL   Вверх
Anikmar
Дата 11.3.2009, 00:10 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
****


Профиль
Группа: Завсегдатай
Сообщений: 2513
Регистрация: 26.11.2006
Где: Санкт-Петербург

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



molodzo, А вы уверены, что уравнение сходится при любом х?

Мне лень было кодить и я загнал вашу формулу банально в excel, и напоролся на х (простите, не помню какой), который в принципе не сходился. Может быть добавить проверку, что при таком-то х, к просто не существует?
PM MAIL ICQ   Вверх
Dmi3ev
Дата 13.3.2009, 10:53 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


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

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



Проблема в том, что мы  считаем (ты говоришь, что к-бывает 43):
((2*43)!)^2=(86!)^2=(1*2*3*...*86)*(1*2*3*...*86)
(4*43+1)!=173! - это число точно не почитает, отсюда и ошибки... нужно пользоваться чем-то другим для хранения таких больших чисел... попробуй просто написать прогу:
Код

for(int i=0; i<1000; i++)
  cout<<fact(i)<<endl;

не пробовал! но на вскидку, дальше 100 точно не должен уйти...


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

PM MAIL   Вверх
Anikmar
Дата 13.3.2009, 12:36 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
****


Профиль
Группа: Завсегдатай
Сообщений: 2513
Регистрация: 26.11.2006
Где: Санкт-Петербург

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



Выражение надо упростить и обойтись без факториала. Использовать несколько другой метод.
Посмотрим на формулу повнимательнее и раскроем ее
При к=1
Дробь получится 5!/(2!*2!)
В развернутом виде (1*2*3*4*5)/(1*2*1*2)

При к=2
9!/(4!*4!)
В развернутом виде (1*2*3*4*5*6*7*8*9)/(1*2*3*4*1*2*3*4)

При к=3
13!/(6!*6!)
В развернутом виде (1*2*3*4*5*6*7*8*9*10*11*12*13)/(1*2*3*4*5*6*1*2*3*4*5*6)
Т.е. сразу видно какую часть можно сократить, а подсчет уже производить не втупую вычислением числителя и знаменателя, а последовательным умножением 7 * 8/2 * 9/3 * 10/4 * 11/5 * 12/6 * 13
Никакого переполнения у нас не будет.

И т.д. для любого к
PM MAIL ICQ   Вверх
Dmi3ev
Дата 13.3.2009, 12:59 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


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

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



Anikmar, сильно не вдумывался, но тут, конечно, надо что-то преобразовывать и считать не в лоб, потому как в лоб не выходит стандартными средствами... или нужно уже брать длинные числа и шпарить...


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

PM MAIL   Вверх
Dmi3ev
Дата 14.3.2009, 00:02 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


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

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



что вроде...
только я брал порядок 10^(-5), думаю, ты легко это поправишь...
Код

#include "stdafx.h"
#include <math.h>
#include <stdio.h>
#include <iostream>
using namespace std;

double myright(double x)
{
    return( 0.5 * ( pow(1 - 4 * x, 1.5) + pow( (1 + 4 * x), 1.5) ) );
}

double myleft(double x, int k)
{
    int u1=2*k;//граница 1
    int u2=4*k;//граница 2
    double res=1;//результат
    for (int i=1, j=u1+1; i<=u1, j<=u2; i++, j++)
        res*=(float)j/i;
    res*=4*k+1;
    res*=pow(x, 2 * k);
    return res;
}

int _tmain()
{
    double xb=-0.25;//начало диапазона х
    double xe=0.25;//конец диапазона х, т.е. х из (-0,25; 0,25) 
    
    int k=1;//к, которое и нужно найти...
    double min=1;//текущий минимум разницы левой и правой части, для проверки
    double step = pow(0.1, 5);//шаг
    double eps=pow(0.1, 5);//погрешность
    for (;;)
    {
        double i=xb+step;//текущий х
        while (i<xe)
        {
            cout<<"x="<<i<<endl;
            double l=myleft(i, k);//левая часть
            cout<<"left="<<l<<endl;
            double r=myright(i);//правая часть
            cout<<"right="<<r<<endl;
            double res=fabs(l-r);//разница по модулю
            cout<<"res="<<res<<endl;
            if (res<min)
                min=res;
            cout<<"min="<<min<<endl;
            if ( res<=eps)
            {
                cout<<k<<endl;
                system("pause");
                return 0;
            }
            i+=step;

        }
        k++;
    }
    cout<<k<<endl;
    system("pause");
    return 0;
}


Это сообщение отредактировал(а) Dmi3ev - 14.3.2009, 12:03


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

PM MAIL   Вверх
molodzo
Дата 15.3.2009, 12:54 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

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



Хм... странная программа... А она у тебя работает? Какие числа выдает?
У меня она считает как будто до бесконечности! А если вручную ввести i, то тоже непонятно когда остановится... Я поставил if k >1500 break; и прога у меня посчитала до этого значения за 3 сек!!! странно... 
P.S. Спасибо за участие!

PM MAIL   Вверх
Dmi3ev
Дата 15.3.2009, 13:01 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


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

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



у меня она выдает последние такие:

Добавлено через 11 секунд
ща досчитает... smile

Добавлено через 3 минуты и 46 секунд
x=-0.21262
left=1.28753
right=1.28753
res=4.8107e-006
min=4.8107e-006
2
собственно 2  и есть ответ...

Добавлено через 9 минут и 50 секунд
а ты паузу сделай и посмотри min уменьшается?
в смысле в консоли (я нажимаю правую кнопу мыши или pause и смотрю)


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

PM MAIL   Вверх
molodzo
Дата 15.3.2009, 16:19 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

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



Да! У меня тоже наконец посчитала  smile ! 
А вот почему, если я ввожу x сам(допустим x = 0.21), у меня прога выдает другие значения! И left при этом равняется 1.#IND ?
PM MAIL   Вверх
Dmi3ev
Дата 15.3.2009, 16:46 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


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

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



Цитата

А вот почему, если я ввожу x сам(допустим x = 0.21),

как ты ввод организуешь? я же не вижу, что ты там творишь?


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

PM MAIL   Вверх
molodzo
Дата 15.3.2009, 19:48 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

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



Вот так:

Код

int _tmain()
{
    double xb=-0.25;//начало диапазона х
    double xe=0.25;//конец диапазона х, т.е. х из (-0,25; 0,25) 
    
    int k=1;//к, которое и нужно найти...
    double min=1;//текущий минимум разницы левой и правой части, для проверки
    double eps=pow(0.1, 5);//погрешность

    double i;
    cout << "Vvedite X: ";  
    cin >> i;
    
    for (;;)
    {
            double l=myleft(i, k);//левая часть
            cout<<"left="<<l<<endl;
            double r=myright(i);//правая часть
            cout<<"right="<<r<<endl;
            double res=fabs(l-r);//разница по модулю
            cout<<"res="<<res<<endl;
            if (res<min)
                min=res;
            cout<<"min="<<min<<endl;
            if ( res<=eps)
            {
                cout<<k<<endl;
                system("pause");
                return 0;
            }
        k++;
    }
    cout<<k<<endl;
    system("pause");
    return 0;
}

PM MAIL   Вверх
Страницы: (4) Все 1 2 [3] 4 
Ответ в темуСоздание новой темы Создание опроса
Правила форума "С++ Builder"
Rrader

Запрещается!

1. Публиковать ссылки на вскрытые компоненты

2. Обсуждать взлом компонентов и делиться вскрытыми компонентами

  • Литературу по С++ Builder обсуждаем здесь
  • Действия модераторов можно обсудить здесь
  • С просьбами о написании курсовой, реферата и т.п. обращаться сюда
  • Настоятельно рекомендуем заглянуть в DRKB (Delphi Russian Knowledge Base) - крупнейший в рунете сборник материалов по Дельфи


Если Вам понравилась атмосфера форума, заходите к нам чаще! С уважением, Rrader.

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


 




[ Время генерации скрипта: 0.1013 ]   [ Использовано запросов: 21 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.