Нелинейные уравнения

[prikolist]

Здравствуйте программисты! Дали решить 2 нелинейных уравнения,фото ниже. Нужно отделить графически или аналитически корни уравнений(найти отрезок [a,b] на котором находится только один корень) и уточнить по одному корню для каждого уравнения с точностью E=0.001 c помощью:
-метода половинного деления;
-метода хорд;
-метода касательной;
-комбинированного метода;
-метода итераций.

Я попробывал методом половинного деления, проверьте пожалуйста, просто, мне нужно чтоб выдавало результат 2-х интегралов, а выдаёт одного, и не находит корня.
Вот исходник:

#include <vcl>
#include <iostream>
#include <cmath>
using namespace std;
const double epsilon = 1e-2;

double f(double x)
{
    return 9.25*pow(x,4)-3*x*x+6.25*x;
}

int main()
{
    double a, b, c;
    a = 0;
    b = 2;
    while (b - a > epsilon){
        c = (a + b) / 2;
        if(f(b) * f(c) < 0)
            a = c;
        else
            b = c;
    }
    cout << (a + b) / 2 << endl;
    system("pause");
    return 0;
}

Вот 2 уравнения,которые нужно решить:

[Sazary]

Вот тут:

while (b - a > epsilon){
        c = (a + b) / 2;
        if(f(a) * f(c) < 0)
            a = c;
        else
            b = c;
    }

Сообщение от prikolist

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

В каком смысле? У вас программа для нахождения одного корня одного уравнения (правого). Если нужно еще и первое, то добавьте еще одну функцию и прокрутите цикл еще раз (а лучше сделать функцию, которая все это будет делать и передавать ей две другие функции по ссылке).

и не находит корня.

Что значит не находит?

[prikolist]

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

#include <vcl>
#include <iostream>
#include <cmath>
using namespace std;
const double epsilon = 1e-2;

double f1(double x);
double f2(double x);
void integral(double (*f) (double) );

int main()
{
cout<<"   --------Metod polovinogo delenia---------";
integral(&f1);
integral(&f2);
system("pause");
return 0;
}

double f1(double x)
{
    return 9.25*pow(x,4)-3*x*x+6.25*x;
}

double f2(double x)
{
return 2*x-log(x)-7;
}

void integral(double (*f) (double) )
{
    double a, b, c;
    a = 0;
    b = 2;
    while (b - a > epsilon){
        c = (a + b) / 2;
        if(f(b) * f(c) < 0)
            a = c;
        else
            b = c;
    }
    cout << (a + b) / 2 << endl;
    system("pause");
  
}

[Sazary]

Во-первых, вы не исправили то, на что я указал.
Во-вторых, вот что у меня выводит:

0.00390625
1.00391

Первый похож на ответ. А для второго, видимо, нужно выбрать другие начальные границы.

[prikolist]

А как для другого изменить начальные, границы, в каком месте,ведь эти границы и изменятся для обоих уравнений. И ещё,у меня вывело вот так результат,странный какой-то,какая-то ошибка log: SING error, и корня целого не вижу. Вот фото:

[Sazary]

Измените функцию вот так:

void integral(double (*f) (double), double start, double end )
{
    double a, b, c;
    a = start;
    b = end;
    while (fabs(b - a) >= epsilon)
    {
        c = (a + b) / 2;
        if(f(a) * f(c) < 0)
            a = c;
        else
            b = c;
    }
    cout << (a + b) / 2 << endl;
    cout<<f((a+b)/2)<<endl;
}

И при вызове посылайте границы:

integral(&f1,0,2);

[prikolist]

Попробывал вот так, но останавливается курсор возле первой ссылки. Вот как написал:

#include <vcl>
#include <iostream>
#include <cmath>
using namespace std;
const double epsilon = 1e-2;

double f1(double x);
double f2(double x);
void integral(double (*f) (double) );

int main()
{
cout<<"   --------Metod polovinogo delenia---------"<<endl<<endl;
integral(&f1,0,1);
integral(&f2);
system("pause");

}

double f1(double x )
{
    return 9.25*pow(x,4)-3*x*x+6.25*x;
}

double f2(double x)
{
return 2*x-log(x)-7;
}

void integral(double (*f) (double), double start, double end )
{
    double a, b, c;
    a = start;
    b = end;
    while (fabs(b - a) >= epsilon)
    {
        c = (a + b) / 2;
        if(f(a) * f(c) < 0)
            a = c;
        else
            b = c;
    }
    cout << (a + b) / 2 << endl;
    cout<<f((a+b)/2)<<endl;
}

[Sazary]

Ну а прототип функции-то, вообще-то, тоже поменять нужно.
И сделайте точность побольше. Например, 1e-4.

[prikolist]

Прототип прописал, но в этом месте,останавливается уже на &f2.

integral(&f1,0,2);
integral(&f2);
system("pause");

Если прописать integral(&f2,0,2);
выдаст такой результат:
(Но нужно только 2 цифры чтоб было в результате, 1 корень,другой не корень)

[Sazary]

Попробуйте что-нибудь такое:

integral(&f2,0.001,8);

И у вас в функции f2 ошибка - десятичный логарифм не log, а log10.