Форум программистов, компьютерный форум CyberForum.ru

Помогите Исправить Ошибку - C++

Восстановить пароль Регистрация
 
Алексей_123
0 / 0 / 0
Регистрация: 22.09.2011
Сообщений: 51
09.11.2011, 16:14     Помогите Исправить Ошибку #1
Помогите пожалуйста исправить ошибку в проге:
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
87
#include <iostream>
#include <cmath>
using namespace std;
double f1(double);
double f2(double);
double f3(double);
double B(double,double, double, double (*f)(double));
double E(double,double, double, double (*f)(double));
 
    int main(){
    double a0, a1, eps;
    double (*f_array[])(double ,double, double, double (*)(double))={B,E}; 
    double (* integ [])(double)={f1, f2, f3}; 
    int a, u;
        cout<<"Input eps "<<endl;
        cin>>eps;
        cout<<"Input a0, a1 for function "<<endl;
        cin>>a0>>a1;
        cout<<"Input the number of methods from 2, 5 "<<endl; 
        cin>>a;
        cout<<"Input the number of integral from 1 to 3 "<<endl;
        cin>>u;
        printf("%.10f",(f_array[a-1])(eps, a0, a1, integ[u-1]));
        }
 
 
double B(double eps,double a, double b, double (*f)(double))
{
    int n=4; 
    double h=(b-a)/4,k1, k0, l;
        for(int i=1; i<=n; i++)
            {
                l=a+i*h; k1+=h*f(l);
            }
                do{
                    k0=k1;
                    k1=0;
                    n*=2;
                    h=(b-a)/n;
                    for(int i=1; i<=n; i++)
                    {
                        l=a+i*h;
                        k1+=h*f(l);
                    }
            } 
                while (abs(k1-k0)>eps);
                return k1;
}
 
double E(double eps,double a, double b, double (*f)(double))
{
    int n=4;
    double h=(b-a)/4, k1, k0, l, m;
        for(int i=1; i<=n; i++)
        {
            l=a+(i-1)*h; k1+=h*f(l);
        }
            do
            {
            k0=k1;
            k1=0;
            n*=2;
            h=(b-a)/n;
                for(int i=1; i<=n; i++)
                {
                    l=a+(i-1)*h;
                    m=a+i*h; 
                    k1+=(m-l)/6*(f(l)+4*f((l+m)/2)+f(m));}
                } 
            while (abs(k1-k0)>eps);
            return k1;
}
 
double f1(double x)
{
    return (x*sqrt(1+x));
}
 
double f2(double x)
{
    return (asin(sqrt(x))/(sqrt(x*(1-x))));
}
 
double f3(double x)
{
    return (pow(2,x)/(1+pow(4,x)));
}

Смысл программы в том, что она вычисляет значение интеграла с заданной точностью eps.
Мне нужно высчитать его двумя методами:по формуле правых прямоугольников(в проге метод 2) и по формуле Симпсона (параболических трапеций)(в проге метод 5).
Если выбирать формулу правых прямоугольников, то прога считает все 3 интеграла, а на формулу симпсона выскакивает непонятная ошибка.
a0:2, a1:7(для первого интеграла), a0:0.2, a1:0.3(для второго интеграла), a0: -13, a1: -2(для второго интеграла).
Приближенные значения: 54.551508(для первого интеграла), 0.12101312(для второго интеграла), 0.35325339(для второго интеграла).

Добавлено через 2 часа 0 минут
eps приблизительно 10^5

Добавлено через 1 час 20 минут
мб вывод данных неправильный?

Добавлено через 18 часов 37 минут
ну гляньте пжл
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
09.11.2011, 16:14     Помогите Исправить Ошибку
Посмотрите здесь:

помогите исправить ошибку C++
Помогите исправить ошибку C++
C++ Помогите исправить ошибку
Помогите исправить ошибку C++
C++ Помогите исправить ошибку.....
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Jupiter
Каратель
Эксперт C++
6542 / 3962 / 226
Регистрация: 26.03.2010
Сообщений: 9,273
Записей в блоге: 1
Завершенные тесты: 2
09.11.2011, 16:19     Помогите Исправить Ошибку #2
Цитата Сообщение от Алексей_123 Посмотреть сообщение
а на формулу симпсона выскакивает непонятная ошибка.
какая именно ошибка?
Алексей_123
0 / 0 / 0
Регистрация: 22.09.2011
Сообщений: 51
09.11.2011, 19:42  [ТС]     Помогите Исправить Ошибку #3
ошибка выскакивает не при компиляции, а уже при самом выполнении программы

Добавлено через 30 секунд
сейчас не могу назвать ее
silent_1991
10.11.2011, 00:12     Помогите Исправить Ошибку
  #4

Не по теме:

Цитата Сообщение от Алексей_123 Посмотреть сообщение
сейчас не могу назвать ее
За вами следят?

Yandex
Объявления
10.11.2011, 00:12     Помогите Исправить Ошибку
Ответ Создать тему
Опции темы

Текущее время: 10:24. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru