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

Корни функции - C++

Восстановить пароль Регистрация
 
Sun4
0 / 0 / 0
Регистрация: 10.03.2011
Сообщений: 3
10.03.2011, 22:09     Корни функции #1
Нужно найти корни функции y(x)=x^4+6x^2-12x-8
Код почти готов, только он находит 1 корень.
Корни которые должны быть: x1=-0,523425 x2=1.78967

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
#include <iostream>
#include <cmath>
using namespace std;
double f( double x ) 
{
    return x*x*x*x+6*x*x-12*x-8;
}
double findRoot(double (*f)(double), double a=1.5, double b=2., double eps=1e-6) {
    double t;
   while (fabs(b-a)>=eps) 
{
      //1 t =( b*f(a)-f(b)*a)/(f(a)-f(b)); 
      t = a+(f(b)*(b-a))/(f(b)-f(a));  //2
      if (f(a)*f(t)<0) 
{
         b=t;
//p         cout << "b=" << b;
      }
else if (f(t)*f(b)<0)
{
         a=t;
//p         cout << "a=" << a;
      } else 
         return t;
   }
   return t;
}
int main() {
   double t = findRoot( f );
   cout << "x=" << t << " f(x)=" << f(t) << endl;
   return 0;
}
Помогите пожалуйста, буду благодарен
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
resource2008
 Аватар для resource2008
109 / 109 / 2
Регистрация: 11.03.2011
Сообщений: 421
12.03.2011, 00:59     Корни функции #2
а это что за метод?

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

попробуй (8-я строка):

findRoot(double (*f)(double), double a=-1, double b=2, double eps=1e-6)
Sun4
0 / 0 / 0
Регистрация: 10.03.2011
Сообщений: 3
12.03.2011, 16:27  [ТС]     Корни функции #3
Метод хорд.
Этот вариант не работает, он вообще показывает не правильный ответ.
resource2008
 Аватар для resource2008
109 / 109 / 2
Регистрация: 11.03.2011
Сообщений: 421
12.03.2011, 22:14     Корни функции #4
хах... точно. нашел.
тебе нужно в этом случае проводить вычисления в два этапа:
первый - диапазон от 0 до 2
второй - диапазон от -1 до 0

нашел у себя старую работу, писалась очень давно на BC еще под DOS
подправил, теперь и на VS пойдет


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
#include<stdio.h>
#include<conio.h>
#include<math.h>
 
 
float function(float x)
{
float h;
//h=tan(0.58*x+0.1)-x;
h=x*x*x*x+6*x*x-12*x-8;
if (h==0) printf("x=%g",h);
return h;
}
 
int main()
{
float e,a,b,c,c1;
unsigned int j;
c=c1=j=0;
printf("Enter correction level: ");
scanf("%f",&e);
printf("Enter a: ");
scanf("%f",&a);
printf("Enter b: ");
scanf("%f",&b);
if (function(b)*function(a)>0) {return -1;}
while (1)
{
j++;
c1 = c;
c=a-(function(a)/(function(b)-function(a)))*(b-a);
if (function(c)*function(a)>=0)
{if (function(c)*function(b)>=0) return -2;
else a=c;}
else b=c;
if (fabs(fabs(c1)-fabs(c))<e) break;
}
printf("x=%2.3f\n",c);
printf("Total: %i steps\n", j);
getch();
return 0;
}
Вложения
Тип файла: doc doc.doc (56.5 Кб, 21 просмотров)
resource2008
 Аватар для resource2008
109 / 109 / 2
Регистрация: 11.03.2011
Сообщений: 421
12.03.2011, 22:15     Корни функции #5
там в файлике чуть-чучть теории, плюс блок схема алгоритма, на всякий случай
Sun4
0 / 0 / 0
Регистрация: 10.03.2011
Сообщений: 3
13.03.2011, 12:35  [ТС]     Корни функции #6
Спасибо большое.
А что бы в один этап не получиться?
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
13.03.2011, 13:58     Корни функции
Еще ссылки по теме:

Найти корни уравнений ax^3+bx^2+cx=0 и y^3-ay^2-b=0, используя процедуру. Есть ли у этих уравнений совпадающие корни? C++
Возвращение значения в функции, ссылки, вычислить корни квадратного уравнения C++
Найти все корни функции f(x) на интервале [-1;1] методом дихотомии с погрешностью 0.02. f(x)= (e^x)-arccos(x) C++

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

Или воспользуйтесь поиском по форуму:
M128K145
Эксперт C++
 Аватар для M128K145
8272 / 3491 / 142
Регистрация: 03.07.2009
Сообщений: 10,707
13.03.2011, 13:58     Корни функции #7
Sun4, посмотрите Нахождение корней уравнения
Yandex
Объявления
13.03.2011, 13:58     Корни функции
Ответ Создать тему
Опции темы

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