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

Численное решение уравнения - C++

Восстановить пароль Регистрация
 
akaelwins
1 / 1 / 2
Регистрация: 11.04.2012
Сообщений: 17
20.05.2012, 13:27     Численное решение уравнения #1
Добрый день, имеется следующее уравнение:
http://www.cyberforum.ru/cgi-bin/latex.cgi?sin(t) = b*{e}^{v*t}
b и v - константы, v > 0
Решаю его численно методом Ньютона, и в принципе, всё хорошо, 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
#include <math.h>
#include <iostream>
#include <iomanip>
#include <algorithm>
using namespace std;
 
double f(double t, double b, double v)
 {
  return sin(t) - b*exp(v*t);
 }
 
double dfdt(double t, double b, double v)
 {
  return cos(t) - b*v*exp(v*t);
 }
 
int main()
{
double b = 0.5;
double v = 0.2;
double t[] = {0, 0};
const int count=100;
double fa=f(t[0], b, v);
double e=10e-15;
for(int i=0; i<count && fabs(fa) > e ; i++)
 {
  
 t[0] = t[0] - f(t[0], b, v) / dfdt(t[0], b, v);
 fa= f(t[0], b, v);
 
 cout<<"t["<<i<<"]= "<< setprecision(10) << t[0] <<"  f("<< t[0] <<")="
                     << setprecision(10) << fa << endl;
 }
cout<<"result: t1 = "<< t[0] <<endl;
system("pause");
return 0;
}
Добавлено через 30 минут
Вопрос снимается, сам догнал)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
20.05.2012, 13:27     Численное решение уравнения
Посмотрите здесь:

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

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

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

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