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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Erwin
1 / 1 / 0
Регистрация: 12.03.2011
Сообщений: 8
#1

Численное интегрирование(ошибка при расчете) - C++

30.04.2011, 14:05. Просмотров 512. Ответов 0
Метки нет (Все метки)

Прога считает определенные интеграл на заданном отрезке.
Пользователь вводит начало и конец отрезка (я ввожу 0.1 и 0.6 соответственно в связи с особенностью интегрирования функции), ну и понеслась.
Ошибка возникает в 46 строчке при расчете I.
Трейсил прогу, выяснил, что значение y[1][0] почему то равно -9.2559631349317831e+061, хотя при заполнении массива y[5][5] значение присваивается вполне нормально...

Подскажите, где я накосячил? скрин

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
#include <iostream>
#include <conio.h>
#include <math.h>
using namespace std;
void main()
{
    setlocale(LC_ALL,"Russian");
    int i, j, k;
    double xn, xk, x[5], y[5][5], h, I, iz;
 
    cout << "Введите начало отрезка (нижний предел): ";
    cin >> xn; cout << endl;
    cout << "Введите конец отрезка (верхний предел): ";
    cin >> xk; cout << endl;
    
    h = (xk - xn)/4;
 
    cout << "Шаг = " <<h <<endl <<"\n";
 
    for (i=0; i<5; i++)
    {   
        x[i] = xn + i * h; 
    }
        
    for (i=0; i<5; i++)
    {   y[i][0] = asin(sqrt(x[i]));
        
        cout << "x["<<i<<"] = " << x[i] << cout<< "\ty["<<i<< "] = " << y[i][0] << endl;
    }
    cout << "\n\t\tТАБЛИЦА РАЗНОСТЕЙ" <<endl <<"\n";
    cout << "\t1\t\t2\t\t3\t\t4" <<"\n";
    k = 0;
    for(i=0; i<5; i++)
    {   for(j=0; j<5-k; j++)
        {y[j][i+1] = y[j+1][i] - y[j][i];}
        k++;    
    }
    k = 0;
    for(i = 0; i<5; i++)
    {   for(j=1; j<5-k; j++)
        {cout << y[i][j]<<"\t" ;}
        cout<<endl;
        k++;
    }
 
    I=h*((y[0][0]+y[4][0])/2.0+y[1][0]+y[2][0]+y[3][0]);
    cout << "\nПриближенное значение = " <<I;
 
    iz=(asin(sqrt(x[4]))*x[4]+1/2.0*((x[4]*sqrt(-x[4]*(x[4]-1)))/(x[4]*x[4]-x[4]*(x[4]-1)))+1/2.0*atan(sqrt((-x[4]*(x[4]-1)))/x[4]))-
        (asin(sqrt(x[0]))*x[0]+1/2.0*((x[0]*sqrt(-x[0]*(x[0]-1)))/(x[0]*x[0]-x[0]*(x[0]-1)))+1/2.0*atan(sqrt((-x[0]*(x[0]-1)))/x[0]));
    cout << "\n\nИстинное значение = " <<iz;
    cout << "\n\nПогрешность = " << iz-I;
    _getch();
    }
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
30.04.2011, 14:05     Численное интегрирование(ошибка при расчете)
Посмотрите здесь:

ошибка при расчете! найти ее - C++
#include &lt;iostream&gt; #include &lt;math.h&gt; using namespace std; int main() { double x; cin &gt;&gt; x; double sum = ((x-2)/(x-1)); ...

Численное интегрирование системы дифференциальных уравнений методом Рунге - Кутта - C++
Доброго времени суток. Пытаюсь разработать функцию для численного интегрирования систем дифференциальных уравнений методом Рунге-Кутта. При...

Вывод nan при расчете - C++
При расчете codeblocks выводит nan как сделать так, чтобы он считал это ошибкой? #include &lt;iostream&gt; #include &lt;math.h&gt; using...

Переполнение при расчете функции - C++
рассчитываю функцию, получается переполнение, что делается в таких случаях? // raschet # include &lt;iostream.h&gt; # include...

Интегрирование, ошибка - C++
При вводе данных в консоль выскакивает ошибка: 'Invalid floating operation' и выделяет часть кода: if (error==0)...

Зачем в целом числе писать вещественную часть при расчете триганометрических функций в библиотеке math.h - C++
Доброго времени суток. В чем собственно петрушка. При расчете функции наткнулся на ошибку вычислений. В результате было обнаружено...

Простое интегрирование по времени - C++
while (t&lt;=interval) { v=vp+((1-pa/p)*g-(9/2)*na*vp/(r*r*p))*dt; z=zp+vp*dt; ...

Интегрирование дифференциального уравнения - C++
Задание: Подготовить программу интегрирования дифференциального уравнения у’=cos(x+1)*tg(x)+y на отрезке 0&lt;=x&lt;=1 при у(0)=0. Интерфейс...

Численное дифференцирование основанное на методе Лагранжа - C++
Подскажите, пожалуйста, как реализовать программу которая осуществляет численное дифференцирование методом Лагранж. На вход программе...

Интегрирование методом левых прямоугольников - C++
Задание: &quot;В прикладной программе проводится функция f(x)=cos(x), пределы интегрирования А и В. Решение провести методом левых...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru